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1. Introduction 


Until not long ago continuum mechanics meant to most people the theories 
of inviscid and linearly viscous fluids and of linearly elastic solids. However, 
the behavior of only few real materials can be described adequately by these 
classical theories. Experimental scientists, who had to deal with real materials, 
developed a science of non-classical materials called rheology. But they did not 
succeed in fitting their experimental results into a general mathematical frame- 
work. Most of the rheological theories are either one-dimensional, and hence 
appropriate at best for particular experimental situations, or are confined 
to infinitesimal deformations, in which case they are only of limited use, because 
large deformations occur easily in the materials these theories are intended to 
describe. 

In the last few years a number of mathematicians have succeeded in devising 
general three-dimensional theories which are valid for large deformations*. 
Guided directly by physical experience or by the one-dimensional theories of 
the rheologists, these authors have proposed various special hypotheses intended 
to characterize the mechanical behavior of large classes of materials in arbitrary 
motions. Using the principle of objectivity of material properties (cf. Section 11 
of this paper) and usually also the assumption of material isotropy ** they then 
derived the general constitutive equation (also called rheological equation of state 
or stress-strain relation) compatible with the original hypothesis. In many cases 
they also solved particular problems in order to make possible a comparison 
with experimental results ***. 

This paper does not start with a special hypothesis. It starts (Chapter II) 
with a general principle, called the principle of determinism for the stress, which 
is implicit in all physical experience and applies to any material whatsoever as 
long as only its mechanical behavior is considered and as long as there are no 
constraints. From this principle and the principle of objectivity the most general 
constitutive equation for all materials is derived. 

Constitutive equations characterize material properties. But material proper- 
ties are of a local nature and may change from particle to particle in a body. 
It turned out to be necessary, therefore, to develop concepts which describe the 
local behavior of a motion at a particular particle. This is done in Chapter I. 
Some of the concepts introduced there are similar to those of the theory of 
“jets” introduced by EHRESMANN [2] into modern differential geometry. The 
treatment of strain and related notions given in Chapter I is, I believe, more 
concise and direct than previous treatments. It is based on the well known 
polar decomposition theorem for linear transformations. 

With only one exception **** the ideal materials previously considered in the 
literature are special cases of what is called a simple material in Chapter III of 


* A review of the work done up to 1953 was given by TRUESDELL [12]. For 
more recent researches we refer to [7], [4], [5], [8], [9], [138]. 

Sas These two things were not clearly distinguished for some time. A clarifica- 
tion was given by the author in [5], where the term “principle of isotropy of space’ 
was used for the principle of objectivity. 

*** Cf. [3], [5], [10], [12], [13], (14). 
*x*x* TRUESDELL’S theory of Maxwellian fluids (cf. [12], Chapter V D). 
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the present paper*. There is reason to believe that most real materials are 
adequately covered by this category. A detailed discussion of the theory of 
these simple materials is given. It includes new and very general definitions of 
isotropic and anisotropic solids and of fluids in terms of the special invariances 
of the corresponding constitutive equations. Fluids are always isotropic in the 
sense in which the term is used in this paper. I believe that physically observed 
anisotropies in fluids are not really intrinsic in these fluids but depend on the 
past history of the substance. 

In Chapter IV it is shown how two large classes of ideal materials fit into 
the general framework. In the isotropic case, these two classes were introduced 
by RIviin & ERICKSEN [8] and Cotrer & RIVLIN [1]. The general non-isotropic 
case is discussed here, too. 

The theory of simple materials presented in Chapter III generalizes the 
memory theory of linear visco-elasticity. The theory of materials of the rate 
type (Section 24) generalizes the differential operator theory of linear visco- 
elasticity. The materials of the differential type (Section 23) correspond to the 
special case when only the stress and not the stress rates occur in the basic 
equation. 

Notation and terminology. Physicists customarily use symbols to denote 
physical quantities. They do not use symbols for the functions that relate these 
quantities, except sometimes the symbol /, which then stands for any function. 
Pure mathematicians, on the other hand, distinguish carefully between functions 
and their values, and they use different symbols for different functions. Both 
methods have advantages and disadvantages. I have found it necessary here to 
follow the pure mathematicians. Thus, in this paper, / means something entirely 
different in nature from /’(t). The first is a function, and the second a value of 
that function and hence a constant. 

The terms “function” and ‘“‘functional’’ are used here as synonyms for 
“mapping”. A mapping is an operation that assigns to each element of one set, 
called the domain of the mapping, a value in another set. A mapping that has 
the same values as a given mapping « on a subset of the domain of « but is not 
defined outside this subset is called the restriction of « to the subset in question. 
Sometimes, if no confusion can arise, we use the same symbol « for such a 
restriction. 

If « and # are two suitable mappings, then «of denotes their composition 
defined by («0 8) (X) =a (8(X)). The inverse of a one-to-one mapping « is denoted 


by a. If F and G are two sets, then ¥ XY means the set of all ordered pairs 
of elements of F and Y. 

In the case of linear transformations /,G of a vector-space into itself we 
omit the symbol o and write simply FG for their composition. The inverse of a 
linear transformation G is denoted by G+, its transpose by G’. The identity 
transformation is denoted by 7. A linear transformation Q is orthogonal if and 
only if QQ7=I, i.e, Q?=Q. The term tensor is used as a synonym for linear 
transformation. Tensors of order higher than two do not occur in this paper. 


* More general materials are those “‘of order m”, discussed by the author in 
[7]. They include TRUESDELL’s theory of Maxwellian fluids. Simple materials are 
those of order one. 
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Points and vectors of a Euclidean space are denoted by boldface letters, 
except when they are values of mappings. If # and y are two points, then y—ax 
denotes the vector leading from a to y. If # is a point and w a vector, then 
a-+u denotes the point y uniquely determined by u=y—«x. 


I. Local kinematics 


2. Basic concepts 


Precise axiomatic definitions of the basic concepts used in this paper are 
given in [6]. Here is a brief summary: 

A body & is a smooth manifold of elements X, Z, ..., called particles. The 
configurations y,%,... of B are the elements of a set of one-to-one mappings 
of Z into a three-dimensional Euclidean point space &. The vector space as- 
sociated with & will be denoted by VW. The mass-distribution m of Z is a measure 
defined on all Borel subsets of #. If g is a configuration of Z, then there is a cor- 
responding mass density o, such that 
(2.1) m(€) = J e,(X)dV 

p(@) 
for all Borel subsets @ of Z. A motion is a one-parameter family {#,}, —co<t<oo, 
of configurations. The parameter Tt is called the t#me. A motion will often simply 
be denoted by #, which is then regarded as a point-valued function of two 
variables, a particle and a time. The mass density corresponding to the con- 
figuration #, is denoted by @(t). 

A body, its configurations, its mass-distribution, and its possible motions are 
subject to the axioms given in [6]. These axioms insure that the customary 
notions and operations are meaningful. However, in some of the considerations 
of this paper, continuity and differentiability conditions stronger than those 
implied by the axioms of [6] have to be assumed. It will be clear from the 
context when this is the case, and we shall not mention these conditions explicitly. 


3. Deformations and linear transformations 


A smooth homeomorphism 6 which maps a neighborhood of the null-vector O 
of the vector space Y onto another such neighborhood and which maps O into 
itself, 


(3.1) 6(O) =O, 


shall be called a local homeomorphism of ¥. We define an equivalence relation 
among all local homeomorphisms by local identity: 6 ~oifand only if 0(P) == 6(P) 
for all P in some neighborhood of O, however small. The resulting equivalence 
classes A will be called deformations. The equivalence class of 606 depends only 
on the equivalence classes A and A of 6 and 6, and it will therefore be denoted 
by MoA. With the law of composition thus defined, the deformations form a 
group denoted by &. 
A linear transformation G, i.e., a mapping of Y onto itself such that 


(3.2) G(R. +P) =G(B) +G(B), G(aP) =aG(P), 
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is a local homeomorphism provided that it is regular, 7.e., provided that GP =O 
is possible only for P =O. Two different linear transformations are never equi- 
valent in the sense defined above. Hence a regular linear transformation defines 
a unique deformation and may thus be regarded as a special deformation. The 
regular linear transformations thus form a subgroup of the group of deformations 
QJ, called the linear group and denoted by ¥Y. 

The gradient at P =O of a local homeomorphism 6 is a local property of 6 
and hence depends only on the equivalence class A to which 6 belongs. We can, 
therefore, define the gradient VA of a deformation A by 


(3.3) VA=VO6(O) if d5EA. 


The gradient VA is a regular linear transformation and hence itself a deforma- 
tion. Since, by the chain rule, 


(3.4) V (404) =(VA) (VA), 


it follows that the gradient operation V is an endomorphism of the group J. 
Its kernel is the group / of all null-deformations, 1.e., all deformations whose 
gradient is the identity 7. The quotient group Y/.V consists of the equivalence 
classes of deformations defined by the equivalence relation 


(3.5) A~A ifandonlyif VA=VA. 


The deformations in each class have the same gradient and any representing 
local homeomorphisms differ only by small terms of order one. The relation (3.5) 
may be read: A and A are equal up to a small term of order one. Each class 
has a unique representative which is a linear transformation. We have 


(3.6) A~VA, 


z.e. any deformation differs from its gradient only by a small term of order one. 
In the case when A =G is a linear transformation, we have 


(3.7) VG=G. 


The quotient group Y/.V is isomorphic to the linear group &%. 

A deformation A will be called isochoric if there is a volume preserving local 
homeomorphism 6 in the class 4. The isochoric deformations form a subgroup 
F% of J. The null-deformations are isochoric, and / is a subgroup of ¥. A linear 
transformation G is isochoric if and only if 


(3.8) | det G] =1. 


The isochoric linear transformations form a subgroup W of Y, called the wnt- 
modular group. They are also called unimodular transformations. The group of 
orthogonal transformations is a subgroup of %; it will be denoted by @. The 
space of all symmetric linear transformations will be denoted by ¥. The set of 
all positive definite and symmetric linear transformations will be denoted by -¥,. 
It is a subset of the linear group &. 
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4. Local configurations 
Consider a neighborhood W(X) of a particle X in a body &, 7.e., a part of 
B containing X. Let y be a smooth homeomorphism of (X) into the vector 
space Y mapping X itself into the null-vector, 


(4.1) p(X) =O. 


We define an equivalence relation among all such homeomorphisms by the con- 
dition that p~y if and only if y(Z) =(Z) for all Z in some neighborhood of X. 
The resulting equivalence classes ® will be called the local configurations of X. 
The set of all local configurations of X will be denoted by @x. 

Let 6 be a local homeomorphism of ¥. The equivalence class of the mapping 
doy depends only on the equivalence class ® of y and on the deformation A of 6. 
Therefore, it is meaningful to speak about the local configuration 4o@ of X, 
which is the composite of the local configuration ® of X and the deformation 4. 
If ® and @ are two local configurations of X, we define 


Ss ey S 
(4.2) A=@o@ if P=AoP@. 
This A will be called the deformation from the configuration ® to the configuration D. 


Let » be a configuration of the body 4. Then y=g—q(X) maps Z into 
the vector space ¥ and has the property (4.1). Its equivalence class will be 
denoted by ®(m, X), and it will be called the localization at X of mw. Assume 
that a motion #? of Z& is given. We then use the notation 


(4.3) O(t) = O(0,, X), 


and we call © the localization at X of the motion % We note that @ is a func- 
tion of a real variable whose values are in @,. Any sufficiently smooth function 
of this type will be called a local motion of X. 


Let @ be a local motion and @ a local configuration of X. We then write 


(4.4) Ome Ome. 


and call 2 the deformation function of X in the local motion @ relative to the local 
reference configuration ®. Q is a function of a real variable with values in ZY. 
Any sufficiently smooth function of this type will be called a deformation function. 


In the special case when ® = @(t) in (4.4), where ¢ is a particular time, we use 
the notation 


(4.5) Q,(z) = (2) 0 Oi 


and call 2, the deformation function of X in @ relative to the time ¢. 


Let y be a configuration of the body Z and 0, the corresponding mass density. 


The value g, (y(X)) depends only on the localization @(y, X) of p at X. We can, 
therefore, define 


(4.6) 02 =0(9(X)), G=Dy,X), 


and call og the mass density of X in the local configuration @, 
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5. Gradients 
Assume that two local configurations ® and D of a particle X are given. 
The regular linear transformation 
Be i —1 
(5.1) z» G=G(®,0)—VA,-A=Do® 


will then be called the deformation gradient from @® to OD. 


In the case when a local motion @ and a local reference configuration ® are 
given, we use the notation 


(5.2) Faj\=V 2G) =G4(O(r),®); 

633) B®) =P (7) =6 (00), 00). 

Assume that @® is another local reference configuration, and let 

(5.4) F(t) =V Q(t) =G (O(n), ). 

Then, by the chain rule for gradients, we have 

(5.5) F(t) =F(2)G, 

where G=GG@, @) is defined by (5.1). We note the following formulas: 
(5.6) G(®, 0) =G(G, 6), 

(5.7) F(t) =E(t)F(), Bl) =I. 


We define an equivalence relation among all local configurations BE Gy, by 
(5.8) @~@® ifandonlyif G(O,0) =I. 


The corresponding equivalence classes M will be called the configuration gradients 
at X. The class of local configurations equivalent to a given DE Gy is denoted 
by V@, so that 


(5.9) M=V®©® _  ifandonlyif @cM. 


The equivalence class of the local configuration 4 o ®, where A is a deformation, 
depends only on the equivalence class of ® and on the gradient VA of A. There- 
fore, it is meaningful to speak about the configuration gradient GM, which is 
the product of the configuration gradient M and the regular linear transforma- 
tion G. This product is characterized by the property 


(5.10) V(Ao®) = (VA) (V®). 

If M and M are two configuration gradients, we define 
(5.44) G=MM if M=GM. 

We have x uw ¥ : 
(5.12) G(®,®) =MM41, M=V0%, M=D2®, 


which shows that G(®, @®) depends only on the gradients of ® and @. In par- 
ticular, the deformation gradient function F in (5.2) depends only on V@, and 
we may call 


(5.13) Lt VO (a) ao 
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the deformation gradient at time 7 of the local motion O relative to the con- 
figuration gradient M as a reference. 

If @ and @ are two local configuration at X, then the corresponding densities 
0» and gz are related by 


(5.14) 06= | det G(®, P)| 09. 


It follows from (5.8) that og = e@il @ and @ have the same configuration gradient. 
Hence we can define 


(5.15) ou=ep if MVE, 


and we call o,, the mass density for the configuration gradient M at X. Equation 
(5.12) then implies ; 

(5.16) 0¢m =| det G| oy. 

It follows that 

(5217) On = Oy ad On yan G MM is unimodular. 


6. Rotation and strain tensors 


Of fundamental importance for the local kinematics of continuous media is 
the well known polar decomposition theorem: 

Let F be any regular linear transformation. Then there are unique decom- 
positions 


(6.1) F=RU.- P= 


in which R is orthogonal and U and V are symmetric and positive definite, 7.e., 
REO, U.VEY,. The following relations are valid: 


(6.2) U2=FTF, V? = FFT, 
(6.3) V=RUR?, V2=RU?R‘. 


. This theorem, when applied to a deformation gradient F =G(®, @), gives 
rise to the following terminology: We call R the rotation tensor, U the right strain 


tensor, and V the left strain tensor of the deformation from @ to @. The squares 


(6.4) C=U*=F'F, B=V2?=FFT=RCRt 


will be called the vight and the left Cauchy-Green tensors of the deformation from 
P to @. 


The strain tensors U and V describe adequately what is meant physically by 
strain because their eigenvalues are the principal extension ratios of the deforma- 
tion from ® to ®. However, it is often of advantage to use the Cauchy-Green 
tensors rather than the strain tensors as measures of strain because their com- 
ponents are rational functions of the components of the deformation gradient F, 
while the components of U and V are complicated irrational functions of the 
components of F. Both C and B are, of course, symmetric and positive definite. 


Their eigenvalues are the squares of the principal extension ratios of the deforma- 
tion. 
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In the special case when F is the tensor function defined by (Sr 2)Be Cala G} 
and B will also be tensor functions of a real variable. If F is replaced by the 
tensor function /; defined in (5.3), we use the analogous notation R,, U,, V,, C;, 
and 5,. These functions have the same smoothness properties as F or &, respec- 
tively. We note that 


(6.5) R(t) =R,() =U) =) =C,() = B,() =. 


7. Histories 


Let ¥ be an arbitrary set. The class of all functions with values in ¥ whose 
domain is the negative real axis will then be denoted by ¥*. 


Let « be any function of a real variable with values in ¥. We then define 
aC F* by 
(7.4) a (s\)=a{t-ts) "for ~s=0, 
and we call it the history up to time t of the function «. 


For functions «*¢ ¥* we generally use the notation 


(7.2) Agree (0) 
Then 
(7.3) a = a'(0) = a(t). 


Let @ be a local motion of a particle X, so that O(t) €@,. The corresponding 
history O6'€ @¥, defined by 


(7.4) G(s) =O(¢+s) for sso, 
will then be called the kinematical history up to t of X. If, in addition, a local 


=il 
reference configuration @ is given and if £2 =@o @ is the corresponding deforma- 
tion function, we call Q’¢€ Y*, defined by 


(7:5) Ore Ole) tors 0; 


the deformation history up to ¢ of X. 


8. Rate of strain and spin 
Consider a local motion @ of a particle X. Let & be the corresponding deforma- 
tion gradient function defined by (5.3) and F' the history of & up to time ¢ 
defined according to (7.1). 
We then call* 


(8.1) Ei) =k) =F (0) 

the velocity gradient at time t of the local motion O@. Similarly * 
(n) (n) 

(8.2) E(t) = BG) =H), Hp OeL Oyu 


* A superposed dot denotes the operation of differentiation and a superposed (7) 
its nth iteration. 
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is called the n acceleration gradient at time t of O. We have, by (6.5) and (8.1), 
(8.3) Buel y . aeete 

If @ is the localization at X of a global motion #, then E(t) (Z,(¢)) actually 
coincides with the gradient with respect to position of the velocity (n™ accelera- 
tion) at «=%,(X) and at time ?. 

The polar decomposition 
(8.4) F(t) = R,(t) U;,(t) 
defines the rotation tensor R,(t) and the right strain tensor U,(r) of the deforma- 


tion from the configuration at time ¢ to the configuration at time t. The histories 
R! and U/ up to time ¢ of R, and U, are defined according to (7.1). The tensor 


(8.5) Wit) = R(t) = Ki(0) 
is called the sfim at time ¢, and the tensor 


(8.6) D(t) = U,(t) = U;(0) 
the vate of strain at time ¢. Similarly, we define the n™ spin W,, and the n™ rate 
of strain D,, by 


(nm) (n) 


(8.7) W, (i) =R,(f) = Ri(0), 
(8.8) D, () = U,() = Uf). 


Replacing U, in (8.8) by the right Cauchy-Green tensor C,= U7, we get a tensor 
(n) (n) 

(8.9) A,,() =C,() =G(0), 

which we shall call the 2‘ Rivlin-Evicksen tensor*. By (6.5), (8.5), and (8.6) 

we have 


(8.10) R= Dp = 4521, (ARS DISD 


One could give definitions similar to (8.8) and (8.9) by replacing the right 
strain and Cauchy-Green tensors by their left counterparts. But the tensors 
defined in this way are of little interest. Moreover, for ~=1 one would get 
nothing new, because 


(8.11) D(t) = U,(t) =V,@), 


as is not hard to prove. 


9. Rational expressions for the rates 
It is possible to find explicit expressions for W,, D,,, and A,, as polynomials 


n? 


in the acceleration gradients FE, and their transpositions E7, k=1,2,...,. 
Differentiating C,H" F, [cf. (6.4)] » times, by the product rule we obtain 
(n) n (k) (n—k) 


C(t) = > (2) B7(0) E (2). 


k=0 


x Riviin & ERICKSEN [8] recognized the importance of these tensors and used 
them extensively. 
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Substituting t=¢ yields, according to (8.2) and (8.9), 


9.4 Vil AIA ee 
( ) n ai) k *“n—k 
Since Ey=TI by (8.3), this may be written in the form 
n—1 : 
(9.2) ApH, e, (Qieaee 
k=l 
Differentiating C{) =U; (t) and then substituting t=¢ we get 
n—1 
(9.3) A, =2D,+ >) (3) Di Due, 
k=1 
and hence 
1 n—1 By | 
(9.4) Dy =n — Di (f) Pa Pup 


This is a recursion formula which can be used to find explicit expressions for 
the D,, as polynomials in the A,, k=1,2,...,”. Hence, after substitution of 
(9.2), one would also obtain explicit expressions for the D, as polynomials in 
the FE, and ET, k=1,2,...,n. 

Differentiating /(t) =R,(t) U(r) and then putting +t =¢, we find 


n—1 
(9.5) Bae Dae, a Dede, 
and hence me 
n—1 
(9.6) EG RES DES DY (7) VO ee 
fie 


This is again a recursive formula which permits us to express W, as a polynomial 
in the E, and D,. Since the D, are polynomials in the E, and Ef, we can also 
find expressions for the W, as polanomials in the E, and Pee ee i. 


In the special case n =1 we find 


(9.7) D=44=HE+E), W=3(E~2%). 


Nn? 


are not necessarily skew. Of course, the D,, and A, are all symmetric. For 1 =2 
we find 


It follows that the spin W is skew (W’=—W). The higher spins W,, n> 14 


A, =E, +E +2E7E, 
(9.8) D, =4(A, —$4j) =4 [2,4 Ey —(E+£E7)*)+ETE, 
L[E, — Ey — E7(E —E*)]. 


II. The general constitutive equation 
10. Basic concepts 
A dynamical process for a body @ is defined as a motion # of # coupled with 
a system of forces for at each time T, subject to the principle of linear momentum 
and the principle of angular momentum (c/. [6], Section 5). A system of forces 
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can be split into body forces and contact forces. If sufficient continuity assump- 
tions are made, the contact forces are determined, for each time 1, by a field 
of symmetric stress tensors S(t). If a motion # and smooth fields of symmetric 
tensors S(t) are arbitrarily prescribed, one can always find a dynamical process 
such that S(z) is the corresponding stress tensor field at time t. Appropriate 
body forces can be chosen, for example, by putting the mutual body forces equal 
to zero and by letting the density b of the external body forces be given by 


(10.1) b=b—— divS, 


where # is the acceleration of the motion #. Then Caucuy’s law of motion ([6], 
(5.25)) holds, and the principles of dynamics are satisfied. 

If we subsequently talk about a process for a body # we mean a pair {#, S}, 
where # is an arbitrary motion of Z and S an arbitrary time family of symmetric 
stress fields. We can always adjoin suitable body forces to a given process so 
as to make it a dynamical process. 

A local process for a particle X CZ is defined as a pair {O, S}, where O is a 
local motion of X and S a symmetric tensor function of a real variable. The 
localization at X of a process {#, S} is defined by replacing @ by its localization 
O at X and S by its values for X. We shall use the same symbol S for these 
values and write S(t) =S(t, X). No confusion can arise because, in all local 
considerations, the particle X will be fixed. 

A constitutive assumption is a restrictive condition on the possible dynamical 
processes a body can undergo, and it characterizes the material properties of 
the body. Restrictive conditions on the possible motions alone are particular 
constitutive assumptions; they are called constraints. Examples of constraints 
are rigidity (every possible motion is a rigid motion) and incompressibility 
(every possible motion is isochoric). In this paper, we shall assume that there 
are no constraints, 7.e. that all motions # are possible. We shall investigate 
constitutive assumptions in the form of functional relations between the stress 
S and the motion # of a process {#, S}. Such relations will be called constitutive 
equations. 


11. The principle of objectivity of material properties 
A change of frame is a transformation of space and time specified by a point- 
valued function ¢ of a real variable, a function Q of a real variable whose values 
are orthogonal transformations, and a real constant a. It transforms a pair 
{#, T} consisting in a point # and a time 7 into the pair {x’, t’} given by* 
(11.1) ; x’ = c(t) + Q(t) (x — O), 
(12) it hp 


where O is an arbitrary point, the same for all transformations. Vectors u CV 
transform according to 


(11.3) Ww = Or) aU. 


* TouPin has analyzed these transformations recently [11]; he calls them ‘“‘Eucli- 
dean transformations of space-time’’. 
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If G is a linear transformation and uc VY, then GucY. Hence, by (11.3), 
(Gu)’=Q(t)Gu. Defining the transform G’ of G by (Gu)'=G'w’, we get 
Q(t) Gu=G'Q(t)u, and hence the following law of transformation for second 
order tensors: 


(11.4) Gl BO (ayG Ot) 


Two dynamical processes are called equivalent if they are related by a change 
of frame as made precise in Definition 7 of [6]. Under such a change of frame 
the contact forces transform in an objective manner, 7.e., according to the law 
(11.3). Hence the stress tensor must transform according to the law (11.4). 
The motion # determines the position of the particles and hence transforms 
according to (11.1). Disregarding the body forces, we can say that two processes 
{0, S} and {#, S"} are equivalent if they are related by a change of frame in the form 


w'(Z, v') =e(t) + Q(t) [8(Z, 7) — O], 
(11.5) S'(Z, 0) = Q(t) S(Z, 2) Q7(z), 
T=T—a@. 
Constitutive equations are subject to the following invariance requirement: 


Principle of objectivity of material properties. If a process {%, S} is com- 
patible with a constitutive equation, then also all processes {0', S'} equivalent to it 
must be compatible with the same constitutive equation. 

This is a special case of the general principle of objectivity stated in [6]. 
Its physical meaning is simply that the material properties of a body should 
not depend on the observer, no matter how he moves. 


12. The principle of determinism for the stress 

We ask for guiding principles which will enable us to find the most general 
form of a constitutive equation. The following two are implied by physical 
experience : 

(7) The stress at a particle X should depend only on the physical state of an 
arbitrarily small neighborhood of X. The state of parts distant from X should 
have no direct influence on the stress at X. This condition is implicit in the 
concept of a contact force. 

(17) The physical state of a body at a time ¢ should depend only on its past 
history, 7.e., on what happened to it at times tS</, and not on its future, 7.e., 
on what will happen to it at times t>?. This condition expresses the causality 
of natural processes. 

The physical history of a body consists of several components: its kinematical 
history, its thermodynamic history, its electromagnetic history, its chemical 
history, etc. In reality, each of these components may influence the stress. 
But, in this investigation, we shall disregard all non-mechanical influences. 


Principle of determinism. The stress S(t) at a particle X and at time t is 
determined by the past history of the motion of an arbitrarily small neighborhood 
of X. 
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In precise terms, this principle states that 
(a) the stress S(t) at X is a functional §, of the motion @, 


(12.4) S(t) =8,(Y) 


the domain of %, being the class of all possible motions and its values being 
symmetric tensors, 

(b) for any two motions # and # which coincide in some neighborhood of X 
for times t<t the value of the functional %, is the same, 7.e., 


(12.2) 3, (8) = 3, (8) 
whenever 
(12.3) 8(Z, t) =8(Z, 7) 


for t<t and Z in some neighborhood W(X) of X, however small. 


13. The general constitutive equation 
Not any arbitrary functional %, with the Property (d) is permissible in (12.1), 
because (12.1), as a constitutive equation, is subject to the principle of objectivity 
stated in Section 11. This principle states that (12.1) must hold equally if the 
process {#, S} is replaced by any equivalent process {#’, S’}, 1.e., that 


(4324) S(O = Be (ee) 


must be valid for any process {#’, S’} related to {#, S} by a transformation of 
the form (41.5). Choosing Q0=J, a=0, and e(t) =O—|0,(X) —O] in (41.5); 
we get S’=S, 7r’=1, and &, =O+4/6,—0,(X)]. It follows that S@)—=—S’7) = 
&:(0’) holds and hence that S(t) depends only on the vector-valued function yp 
defined by p(Z, t) =#(Z, t) —0(X, t). Moreover, the Property (b) of %, implies 
that S(¢) can depend only on the local behavior of y in an arbitrarily small neigh- 
borhood of X. This means that S(t) is determined by the localization O at X 
of the motion # as defined by (4.3). Thus, (12.1) reduces to the form 


(13.2) S(t) = &:(0) 

where @, is a functional with the property that 
(13.3) 5:(0) = 3:0) 
whenever 

(13.4) O(t) =O(t) for tt. 


Now we consider another equivalent process by choosing e(t) =O, Q=I, 
and a=t in (11.5). We then have 


Vata tt otto, SH) = Sior—si, 


and, for the corresponding local motion 0, 
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It follows from (13.1) that 
(13.5) S(t) = $"(0) = 3o(0'). 


By (13.4) we have $(0’) =@ (@’) whenever O’(t') =O'(t’) for t’<0, which 
shows that the value %)(0’) depends only on the restriction of the function 
O'(r’) to t’S0. But this restriction is nothing but the kinematical history 
O'C@¥ up to time ¢ as defined by (7.4). Omitting the index 0 in (13.5), we see 
that the constitutive equation (13.2) reduces to the form 


(13.6) S() =3 (0), 


where © is the kinematical history of X for the local motion @. The form of 
the functional % is independent of ¢ and depends only on the particle X. Its 
domain is @*. 

We finally consider an equivalent process by choosing e€(t) =O and a=0 
in (11.5), leaving the orthogonal tensor function Q arbitrary. The principle of 
objectivity then implies that the functional % of (13.6) must satisfy the relation 


(13.7) Qo3(O*) Oo = 3(Q*0O*), Qo = O*(0), 
for all kinematical histories 9* € @* and all orthogonal tensor functions Q* € O*. 

It is not hard to see that, conversely, the principle of determinism and the 
principle of objectivity are automatically satisfied for any equation of the form 
(13.6) provided the functional %} has the property (13.7). Hence (13.6) with (13.7) 
is the most general constitutive equation. An equation of this form restricts 
the class of all possible local processes {9, S} for the particle X and characterizes 
the local material properties of X. 

If we subsequently speak about a particle X, we always assume that a func- 
tional of the type discussed above is associated with it. We call it the functional 
of the particle X. 

14. Material isomorphisms 

The nature of the domain @¥ of the functional % of the particle X varies 
with the particle X. Hence there is no direct way to compare the functionals 
of different particles. It is desirable to render such a comparison possible, because 
only then can a precise meaning be given to the statement that two different 
particles consist of the same material. 

A body was defined in [6] as a certain mathematical structure. As in the 
case with any such structure, it is meaningful to talk about isomorphisms _ be- 
tween bodies. An isomorphism of a body & onto a body F is a one-to-one mapping 
y of B onto # such that 


(a) the configurations @ of B axe of the form 


ad —1 
(14.1) P=P°y, 
(b) the mass distributions m and m of Y and & are related so that 
(14.2) m(€) =m (y(@)) 


for all Borel subsets @ in &. 
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Assume that X and X are two particles and consider the isomorphisms y, 


if any, of neighborhoods 4(X) onto neighborhoods (X) which map X into 


X =y(X). We define an equivalence relation among all such isomorphisms by 
the condition that y~? if and only if y(Z) =y(Z) for Z in some neighborhood 
of X, however small. The resulting equivalence classes J” will be called the local 
isomorphisms of X onto X. If I’ is such a local isomorphism then the local 


configurations ® of X are of the form 

ny —i 
(14.3) DE Pale 
where @ is a local configuration of X. 


Definition 1. A local isomorphism I" of a particle X onto a particle a will 
be called a material isomorphism of X onto X provided the functionals & and § 
of X and X are related by 


(14.4) 3(O*) =F (OFT) 


for all kinematical histories O* € €¢. 
We shall say that two particles consist of the same material if they are 
materially isomorphic to each other. 


15. Constitutive functionals 


Assume that a local motion @ and a local reference configuration ® of a 
particle X are given. Let @’ and @ be the corresponding kinematical history and 
deformation history, respectively, as defined in Section 7. By (4.4) they are 
related by 


(15.4) G =O oD, 
Hence the general constitutive equation (13.6) may be rewritten in the form 
(4522) S(t) = § (Qo D). 


Definition 2. A functional & whose domain is the set Z* of all deformation 
histories and whose values are symmetric tensors is called a constitutive functional 
provided it has the following property: For all deformation histories Q*€Q* and 
all orthogonal tensor functions Q* €O* the relation 


(15.3) Qo 8 (Q*) Qo = G(Q*02*), Qo = O*(0), 
holds. 
It is clear that the functional ®, defined by 


(15.4) Gg (Q*) = G(Q*; D) = F(O* o D) 
is a constitutive functional, because the property (4 5.3) for © =Gg is equivalent | 


to the property (13.7) for %. Hence the most general constitutive equation may | 
be formulated in this way: 
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Theorem 1. Given a particle X and a local reference configuration ® of X, 
there is a constitutive functional ®g such that, for any local process {0, S} of X, 
the stress S 1s related to the local motion © by 


(15.5) S(t) = Go (2) = G(Q'; 6), 


where 92° is the deformation history of X defined in Section?. The constitutive 


functionals Sg and Sg corresponding to the local reference configurations ® and D 
are related, for all deformation histories Q*CQ*, by 


~ 


(15.6) @(Q* oA; B) = G(Q*; D), 


es | es 
where A=Do® 1s the deformation from @ to ®. 


The relation (15.6) is an immediate consequence of (15.4). We call G, the 
constitutive functional of the particle X relative to the local reference configura- 
tion @, 

Since the nature of the domain Y* of constitutive functionals is not related 
to particles, as was the case with the functionals % in (13.6), they characterize 
materials in a manner that is independent of the particular particle. The following 
theorem is a consequence of Definition 1 and (15.4): 


Theorem 2. Two particles X and X are materially isomorphic, i.e., consist 
of the same material, if and only if there ts a local configuration ® of X and a local 
tsomorphism I’ of X onto X such that 


7 cal 
(1527) (5(O*) = OGs(Q0%))) "OS Dias” 
for all Q*CD*. 
16. The local isotropy group 


Let I’ be a material automorphism of X, 1.e., a material isomorphism of X 
onto itself. By Theorem 2, there is a local configuration ® such that 
ES Ex —1 
(16.1) G (Qt; @) = G(Q*;G), O=DPol!. 
n=l Be 
The material automorphism J" and the deformation J =@®o®@ from ® to @ are 
related by 


ale al Sai sail 
(16.2) A=@Golo®@, LT=OoAo®@. 


Since J’, an isomorphism, preserves the mass distribution, it follows that A must 
be an isochoric deformation. 


In addition to (16.1) we have the relation (15.6) between the two functionals 
®, and 3. Combining these two relations we obtain 


(16.3) Gy (Q*o A) = Go(O*). 


=i) =! 
Conversely, if (16.3) holds for an isochoric deformation 4 €%, then ['=®oAo® 
can easily be seen to be a material automorphism of X. 
Arch. Rational Mech. Anal., Vol. 2 15 
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Definition 3. Let ® be a constitutive functional. The group G' of all tso- 
choric deformations AC ¥ with the property that 


(16.4) (6 (Q*o A) = G (Q*) 


holds for all Q* € Y* is called the local isotropy group of ©. 

G' is a subgroup of the group ¥ of all isochoric deformations. The material 
automorphisms of a particle also form a group. The relations (16.2) establish 
an isomorphism of this group with the local isotropy group Y of Gg. The local 
isotropy groups GY, and G;, corresponding to different local configurations ® 


and ® will in general be different. But they are all isomorphic to each other 
because they are isomorphic to oe group of material automorphisms of X. 
In fact, it is not hard to see that G', is conjugate to ee 


(16.5) G, = No Gio At, A=. 


III. Simple materials 


17. Simple constitutive functionals 


It may happen that the values © (Q*) of a constitutive functional © are not 
affected if Q*(s), for each s<O, is changed by a small term of order one. By 
(3.6), Q*(s) differs from V2Q*(s) only by such a small term. 


Definition 4. A constitutive functional & is said to be simple if 
7A) & (.2*) = G6 (VQ*) 
for all deformation histories Q* €D*. 


The condition (17.1) is equivalent to the following: If V.Q*(s) = VQ*(s) 
for all s<O, then 
(17.2) ( (Q*) = B(Q*). 
The values of a simple constitutive functional © are determined for all Q* € D* 
if they are known for the histories of linear transformations F*€ ¥Y*. In other 
words, & is determined by its restriction to #*. 


From now on we assume that & is simple, and we use the same symbol & 
for its restriction to ¥*. By (15.3) ® satisfies the relation 


(17.3) Qo (F*) Qo =G(Q*F*), Qo = Q*(0), 
for all F*€.#* and all Q*E€ O*., 


Let GE be an arbitrary regular linear transformation with the polar 
decomposition 


(17.4) Go Pare Ee lca Gre 
For any I’*¢€ #* we define 


(17.5) F* = F*GA—R*U*, R*cO*, OC ¥?, 
where R* and U* are determined by the polar decomposition of F*, We then have 


F* = F*G = R* U* PT = (R* P) (PT U* P)T. 
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Substituting R* P for Orand she U* P)T for F* in (17.3), we see that 
(17.6) ® (F*) =R, P@((P7O* P)T) PTR, Ry, = R*(0) 


We now define a functional 8 with domain %* x FY, by 
(17.7) OT Vie OUST) 

for all U*C Y* and TE Y,. Equation (17.6) then shows that 
(17.8) @(F*) = Ry P ®(P? U* P; T) PT RI. 

The functional ® has the property that 

(17.9) @(U*;T) =R(U*;T) if U*T=UFT. 


Conversely, if ® is any functional with domain Y%* x Y, and with the property 
(17.9), then (17.8) defines a simple constitutive functional &. The variables on 
the right side of (17.8) are defined in terms of F* and an arbitrary GEL. For 
the special choice G =I we get 


(17.10) Oey = Roe Us Ll) kas 1 hy = (0)s 
where U* and k* are defined by the polar decomposition 
(17.14)) Pea Ue ee OF ORGS. 
In the special case when G=F,y=F*(0) (17.8) becomes 
(17.12) O(A*) = RoR(REG TRS UA) RE 


where R,=R*(0), Uj =U*(0), and where UC Y* is determined by the polar 
decomposition 


(17.13) BE cme eer Fee Uris Te a Carer eg Gt 4 
We define a functional , with domain Y* by 

(17.14) RVC Ai QUAM) —-for-= OF EL 

Then (17.10) takes the form 

(17.15) 6 (F*) = R, &,(U*) RE. 


Conversely, if §, is an arbitrary functional with domain .Y*, then (17.15) defines 
a simple constitutive functional. 


We also define a functional §, as the restriction of & obtained by allowing 
for its first variable only those functions U,.¢ Y* whose value for s=0 is the 
identity, 

(17.16) OF (O\rea i: 

We denote the set of all U,* with the property (17.16) by 4*,. The domain 

of &, is then S*, x F,, and it is identical with ® in this domain. The function 
15* 
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U* defined by (17.13) has the property (17.16), and the same is true for 
Ro UX R,. Hence (17.12) has the form 


(17.47) @(F*) = R, R, (RE U# Ry; Uy) RF. 


Conversely, if &, is an arbitrary functional with domain GE eT TORR er 
defines a simple constitutive functional ©. 


18. Simple materials 


We say that the material at a particle X is simple or, briefly, that X is simple 
if the constitutive functional Gg of X, for some local configuration Dol 2ois 


simple. We assume from now on that this is the case. Let ® be another local 


et a 
configuration of X and let 4 =@®o® be the deformation from ® to ®. For 
every deformation history Q*€Q*, by (3.7) and the chain rule (3.4), we have 


(18.1) V( 2 oA VV 2 eAl: 
It follows from (17.2) that 
&(Q* oA; bd) = G(VQ*oA; D), 
and hence from (15.6) that : ‘ 
& (Q*; O) = G(VQ*; ®), 


which shows that ®g is also simple. Therefore, if the material at a particle is 
simple then the constitutive functional g is simple for all local reference con- 
figurations @. 

It follows from (15.6) that the functionals Gg and &g corresponding to two 
local configurations ® and @ are related by 


x~ 


(18.2) W(F*G; OD) = G(F*; ®), 


where G Go ®) is the deformation gradient from ® to @ defined byat$. 1) 
It follows from (18.2) that 


OG(F*;@) =G(F*;®) if G(b,%) =I. 


This means, according to (5.8), that G(*; ®) depends only on the gradient 
M=V® of the local configuration ©. We can therefore define 


(18.3) Gy (F*) = G(F*;M) =G(F*; 6) if M=V@. 
The following theorem is a corollary of Theorem 1. 


Theorem 3. Let X be a simple particle and let M be a configuration gradient 
of X. Then there is a simple constitutive functional ®y such that, for any local 
process {O, S}, the stress is related to the local motion by the constitutive equation 


(18.4) S() = 6y(F') = G(F*; M), 


where F'€ £* is the history of the deformation gradtent of the local motion © relative 
to M as a reference. The functionals ®y, and ® xz corresponding to two configuration 
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gradients are related by 


aA “aw 


(18.5) O(FLG TM) =O (F*> MM).  G = MM, 
for all F*C L*. 
The results of the previous section enable us to put the constitutive equation 


into various other forms: 


Theorem 4. The constitutive equation for a simple particle X may be written, 
with reference to a fixed configuration gradient M of X, in one of the following four 
forms: 


(a) There is a functional &, with domain S* such that the stress is given by 
(18.6) S(t) =R()R(U) RY’, 


where R(t) is the rotation tensor and U' the history of the right strain tensor relative 
to M as a reference. 

(b) There is a functional & with domain S* x Sand with the property (17.9) 
such that, for any configuration gradient M as a reference, the stress is given by 


(18.7) S(t) =R(t) PR(PTU'P; T) PTR(O!, 

where R (¢) 1s the rotation tensor and U! the history of the right strain tensor relative 
to M as a reference, and where P 1s the rotation tensor and T, the right strain tensor 
of the deformation from M toM, MM1= PT. 


(c) In the special case when M is the gradient of the configuration at some 
tame ty the equation (18.7) takes the form 


(18.8) S(t) = Ry, (0) B (to) 8 (RB (to) * Uy, R (to) 5 U(te)) B (to)? Bi, (1), 
where the rotation tensors R and R, and the right strain tensors U and U,, are taken 
relative to M and the configuration at time ty as a reference, respectively. 

(d) There is a functional &, with domain S*,x SF, such that the stress is 
given by 
(18.9) S(t) = R(t) &,(R? (i) UJ RO); UH) RO, 


where R, U, and U, are defined as before. 


The forms (18.6) and (18.9) have the advantage that the functionals ®, and 
R, are not subject to any restrictive condition. 


19. The isotropy group 
If ®& is a simple constitutive functional, then it follows from (17.2) that 
®& (Q*) =G(Q*o A) holds for all null-deformations 4€.Y, 7.e., whenever VA =I. 
Hence, by definition 3, the local isotropy group Y of & contains the group / 
of all null-deformations as a normal subgroup. The quotient group Y/W is 
isomorphic to a group of unimodular transformations. 


Definition 5. The group GY of all unimodular transformations H with the 
property that 
(19.1) @(F* H) = © (F*) 
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holds for all F*€L* is called the isotropy group of the simple constitutive func- 
tional ®&. 

G is a subgroup of the group Y of all unimodular transformations and it is 
isomorphic to G/N. 

Theorem 5. An orthogonal transformation Q is an element of the isotropy 
eroup GY of a constitutive functional & if and only if one of the following conditions 
as satisfied : 


(a) For all F*€L* 


(19.2) Q&(F*) Q7 = 6 (QF* Q*). 
(b) For all U*CS* and all TE SF, 
(19.3) QR(U*; T) 07 = R(0U* 07; OT QO’), 


where & 1s defined by (17.7). 

(c) For all U*ES* 
(19.4) Q &, (U*) QF = 8,(QU*Q"), 
where &, 1s defined by (17.14). 

(d) For all US CL5, and all TE F, 


(19.5) Q&, (Ug; T) Q7 = 82(Q UF OF; QT Q") 
where &, 1s defined in Section 17. 

Proof. If we substitute the constant tensor Q for the tensor function Q* 
in (17.3), and QF* for F* and Q for H in (19.1), we see that (19.2) holds if and 
only if Q€Y. The equivalence of the condition (a) with (b), (c), and (d) follows 
directly from the definitions of ®, ®,, and &. 

If Y contains the full orthogonal group @, then (19.2), (19.3), (19.4), and 
(19.5) are valid for all orthogonal Q. 


The isotropy group Y, of the constitutive functional ,, of a particle will, 
in general, depend on the choice of the configuration gradient M. But, as in 
the case of local isotropy, the groups @, and Y corresponding to two con- 
figuration gradients M and M are conjugate: 


(19.6) G.=G6G,Gi, G=MM1 
It follows from (18.5) that HE G, if and only if 

(19.7) @(F*; M) = @(F*; AM) 
for alli Cee, 


20. Isotropic and anisotropic solids 


We say that a constitutive functional ®& defines a solid if its isotropy group is 
a subgroup of the orthogonal group, 7.e., if YCO. A particle X is said to be a 
solid particle if there is a configuration gradient M of X such that & m defines 
a solid. A solid is called an isotropic solid if the isotropy group of its defining 
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functional is the full orthogonal group, 7.e., if G=@0. Let X be an isotropic 
solid particle. Any configuration gradient M such that Y%,—0 is then called an 
undistorted state of X. 


The following theorem follows immediately from Theorems 4 and 5. 
Theorem 6. Let X be an isotropic solid particle. Its constitutive equation 


may be written, with reference to an undistorted state M of X, in one of the following 
forms: 


(20.1) S(t) =&, (RO RT), 

(20.2) S(t) =RY) &(0'; T) RW’, 
(20.3) S(t) =R,,() (Uf; Vii) BR, (0), 
(20.4) S(O) = 8,(U5 VG). 


The notation of Theorem 4 applies here. In addition, V=RUR?’ is the left strain 
tensor relative to M as a reference and T, 1s the left strain tensor of the deformation 
from M to M. The functionals &, 8, and &, satisfy the conditions (b), (c), and 
(d) of Theorem 5 for all orthogonal transformations Q. 

A solid is called anisotropic if the isotropy group of its defining functional 
is a proper subgroup of the orthogonal group, 7.e., if CO and Y+0. Material 
symmetries in anisotropic solids, such as orthotropy, transverse isotropy, and the 
various types of crystal symmetry are defined according to the special nature 
of the isotropy group Y. 

21. Fluids 

We say that a constitutive functional defines a fluid if its isotropy group 9 
is the full unimodular group %, 1.e., if G=%. A particle X is said to be a fluid 
particle if, for some configuration gradient M of X, the corresponding constitutive 
functional defines a fluid, 7.e., if G,=W. Let M be any other configuration 
gradient. By (19.6) G; is conjugate to %,=—W. But W is a normal subgroup 
of ¥ and hence coincides with all its conjugates. It follows that, if X is a fluid 
particle, then @,=Y for all configuration gradients M of X. 

For a fluid it follows from (19.7) that 


(21.1) @ (F*; M) = @(F*; M) 

wx all 
whenever H =MM is unimodular. But, by (5.17), this is the case if and only 
if the densities gy, and @, coincide. If follows that the value & (F*; M) can depend 
only on the density 9,,. Therefore, we can define a functional § with domain 
L*x AB, (#, =set of positive real numbers) such that 


(21.2) G(F*; M) = 9(F*; em) 


for all F*€ F* and all configuration gradients M of X. Let F*—kR*U* be the 
polar decomposition of F*. Substituting A* for Q* and U* for F* in (17.3) and 
using (21.2) we see that 


(21.3) @(F*;M) =RyO(U*i 0m) Ro, Ro=R*(0), 
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which shows that § is determined by its restriction to Y* x #,. We use the same 
symbol § for this restriction. Since Y,=®% contains the orthogonal group it 
follows from (19.2) that § satisfies the relation 


(21.4) Q.(U*;d) Q7 = §(Q U* QF; d) 
for all U*C F*, all QEO, and all d>0. 

Theorem 7. The constitutive equation for a fluid particle X may be written 
in one of the following forms: 


(a) There is a functional § with domain S* x B, and with the property (21.4) 
such that the stress 1s given by 


(21.5) S(t) =R() 9(U; om) RO", 
where R(t) is the rotation tensor and U' the history of the right strain tensor relative 
to an arbitrary configuration gradient M with mass density oy. 


(b) In the special case when M 1s the configuration gradient at some time ty, 
the equation (21.5) takes the form 


(21.6) S(t) = Ri, (0) S (UL; 0 (to) Bi, ()- 

(c) There is a functional $, with domain S*, x B, and with the property that 
(21.7) Q 91 (Us; 4) OF = 91(0 UF QF; a) 
for all UZ CSF, yee all QEO, and all d>0, such that the stress 1s given by 


(21.8) S(t) =9: (U5 e(d)- 


Proof. The part (a) and its special case (b) follow from (21.3) and Theorem 3. 
The part (c) follows from (6) by choosing fj=¢ and by defining , to be the 
TESlLICliON LOLI) 10.4 a oe 

We note that an arbitrary functional §, with domain Y*, x &, and with 
the property (21.7) may define a fluid. The constitutive equation (21.8) is 
intrinsic in the sense that it does not depend on the choice of a reference con- 
figuration. The constitutive equation of a solid cannot be put into such an 
intrinsic form. 


22. Constitutive equations involving the Cauchy-Green tensors 
As we pointed out in Section 6, it is often better to use the Cauchy-Green 
tensors C and B instead of the strain tensors U and V. In order to do so, we 
define new functionals &, %,, &,, § and §,, of the same nature as the correspond- 


ing functionals without the superposed bars defined in Sections 17 and 20, by 
the following formulas: 


(22.1) ee T?) = U,R(U*; T) Up, 
cea pe AUNT: 
(22.3) Ry ws Tae We Pye 
(22.4) §(U*2; d) =U, 8(U*; d) Us, 
(22.5) 6(U SD) SoU sa, 
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valid for all UE L* (U=U*(0)), all USE LS, (Uf (0) =J), all TEY,, and 
all d>0. 

It is not hard to see that the general constitutive equations (18.6), (48.8), 
and (18.9) of Theorem 4 for simple materials then take the form 


(22.6) F(t)? S(t) F(t) = &,(C),* 
(22.7) F(t)" Sf) F(t) = R (to) 8 (R (to) Ci, R (to) 5 C(to)) B (to), 
(22.8) Si) =RW)R, (RO? CLR); C®) RW. 


In the case of isotropic solids, the functionals @, 8,, and , satisfy, for all 
orthogonal Q, the functional relations obtained from (19.3), (19.4), and (19.5) 
by superposing bars. Moreover, the simplified constitutive equations (20.3) and 
(20.4) for isotropic solids take the form 


(22.9) F, (i)? S(@) F,() =R (Ch; Blt), 
(22.10) 5 (=n Ce By) 
The functionals § and §, have the same properties (21.4) and (21.7) as the 


corresponding functionals without the superposed bars. The constitutive equa- 
tions (21.6) and (21.8) for fluids take the form 


(22.14) F(t)? S(t) B,) = 9 (C3 e(t)) 
(22.12) SHW=H1(Cs 0). 


IV. Special classes of materials 


23. Materials of the differential type 
The value &(F*) of a simple constitutive functional is determined by the 
values F*(s) of the tensor function F* for s<0. It may happen that &(f*) 
depends only on the values F*(s) for s very near to zero. If F* has sufficiently 
many continuous derivatives then /'*(s) may be approximated, for small values 
of s, by its Taylor expansion up to some order n. This Taylor expansion is deter- 
mined by the value of F* and its derivatives up to the order m at s =O, 1.e., by 


(n) (n) 


(23.1) R=F*(0), K=F*(0),....& =F*(0). 


Definition 6. A simple constitutive functional & is sard to be of the differential 
type tf 


(23.2)  (F*) = 6 (F*) 
whenever 

() (k) 
(29.3) Et (Ojo 2 (0), [Rm=Oe4l, coo nite 


We have seen that every simple constitutive functional & has a representation 
of the form (17.17) in terms of a functional %,. Since S, is just a restriction of & 


* This form, in other notation, has been proposed independently by GREEN & Riv- 
LIN [4]. 
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which is defined by (17.7) in terms of &, we have 
(23.4) Ry (UF; T) = G (US T) 


forall Ue F*, and all Te.7, IP is of the differential type then its value 
depends only on the values at s =O of its argument and its first m derivatives. 
Hence, by (23.4), the value &,(U,f; 7) depends only on the first derivatives 
of U* at s=0, since U#(0) =/. It follows that there is a function f of m+1 
symmetric tensor variables such that 
. * (nN) 

(23.5) Wy (UZ; T) = t(UF(0), UF (0), --., Ue (0); T) 

for all UX e %, and all TE Y,. In the case when Uf =U; is the history up to 
time ¢ of the right strain tensor relative to the configuration at time ¢ of a local 


(k) ; 
motion, the derivative U/(0) coincides by (8.8) with the k™ rate of strain D, (t). 
For simplicity we use the notation 


(23.6) HD. Dk Dee 7) =U cee 


We say that a particle X is of differential type if its constitutive functional 
Gy, for some configuration gradient M, is of the differential type. It is not 
hard to see that Gy, for any other configuration gradient M, is then also of the 
differential type. The following theorem is a consequence of (18.9) and the 
remarks made above: 


Theorem 8. The constitutive equation of a particle X of differential type may 
be written, with reference to a configuration gradient M of X, tn the following form: 


There 1s a symmetric-tensor-valued function £ of m +1 symmetric tensor variables 
such that the stress 1s given by 


(23.7) S() =ROEUR'() D,Q RO; UO) RY’, 
where D,(t) is the k™ vate of strain and where the rotation tensor R (t) and the right 


stvain tensor U(t) are taken relative to M as a reference. 


In a material of the differential type the stress depends only on the immediate 
past of the motion and not on its course at times long ago. 

Since the rates of strain D, can be expressed as polynomials in the accelera- 
tion gradients E, and Ej} as shown in Section 9, it follows that the constitutive 
equation (23.7) is a relation involving the displacement gradient G=RU, the 
acceleration gradients, and the stress. 


In the case of isotropic solids, it follows from Theorem 6 and (23.5) thater 
must satisfy the relation 


(23.8) QED; T) QT =t(QD, Q7; QT Q?) 


for all Di€ F, kh=1, 2,05, all QC, and all TESF,. A function with this 
property is called an isotropic tensor function. By (20.4) the constitutive equation 
(23.7) reduces for isotropic solids to 


(23.9) S() =1(D, (i); Vii)). 
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The constitutive equation for a fluid of the differential type has the form 
(23.10) S(t) =) (PD. ); e() 


where } is an isotropic function of m symmetric tensor variables and one positive 
scalar variable, 7.e., it satisfies the relation 


(23.11) Q5 (Dz; 4) O° =b(Q D, Q7; 4) 


for all D,C FY, k=1, 2,...,”, all d>0, and all QE€O. This is an immediate 
consequence of Theorem 7, (21.8). 

If we use the alternate forms (22.8), (22.10), and (22.12) of the general con- 
stitutive equations we arrive at the following forms for materials of the dif- 
ferential type: In the general case, 


(23.12) Si) =R() ER? A,O RO; CY) RW. 
For isotropic solids *, 

(23.13) S(é) =t(A,(); BQ). 

For fluids, 

(23.14) S(t) =§(A, (5 e@). 


In these equations f and § are of the same type as f and 9; in (23.12) f may be 
arbitrary but in (23.13) it must be isotropic. Of course, is also isotropic. A, (2) 
is the k™ Rivlin-Ericksen tensor, defined by (8.9) and related to the acceleration 
gradients by (9.2). 
24. Materials of the rate type 
The general constitutive equation of a simple material in the form (18.8) 
may be rewritten as 


(24.1) Si, () = & (Ui; U(to)) 

where 

(24.2) S,,(t) = R (to)? Ri, ()7 S(t) Ri, B eo) 
and 

(24.3) Uj. (s) = U,,(¢ +s) =R (ty)? U,, (6+) R(t). 


Keeping ¢, and U(t)) fixed, we may interpret (24.1) in the following manner: 
Assuming that the function U, defined by 


(24.4) U,, (t) =R (t)-G;, © Rib)? 


is given, the function a defined by (24.2) is completely determined. In other 
words, (24.1) defines an operation on functions U; with values in -%, which 
gives functions S, with values in % It may happen that this operation is defined 
by the process of solution of a differential equation for S,, in the form 


(m) (n) 


(24.5) £(Sy (t), a (2) ear o(2))5 U,, (Zt), U,, (é), no Ce (2) U(to)) ==0); 


* This form was first derived by RivLin & ERICKSEN [8]. 
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where f is a symmetric-tensor-valued function of m+u-+3 symmetric tensor 
variables. For simplicity we use a notation similar to (23.6), so that (24.5) 


becomes 
(24.6) #(5,,(); U;,(t); Ult)) =0. 


If the function U, is given, then (24.6) is a differential equation of order m for 
the function S,. We assume that the form of { is such that there is a unique 
: ® 


solution S, which assumes eves initial vanes S:(4), 4=0,1,-..,m—1, no 


matter how we choose U(ty), ee t,, and so (4,). 


Since f, is arbitrary, in (24.6) we may make the special choice y=? obtaining 


(24.7) F(S,(t); U,(); U() = 0. 
By (24.4) and (8.8) we have : 
(24.8) U, (t) = R? (t) D, (t) R(d). 


The tensor function S defined by 


(24.9) S(t) =R,, (7 S(O) R,, (6) 


will be called the /™ invariant stress rate. If we carry out the differentiation in 
(24.9) according to the product rule, and if we observe (8.7), we see that 


(24.10) Si 
PG FAO Dna 
p+qtr=l 
thus S; can be expressed explicitly in terms of the stress S, its time derivatives 
@ e 
S up to the order J, and the spins W, up to the order 7. We have Sg=S. For 


1=1, we get the invariant stress rate* 
(24.11) S=S=]S— WS ESW. 
Observing (24.2), (24.9), and (24.8), we see that (24.7) has the form 
(24.12) F(RO™SO RO; ROD) RW; UW) =0. 
A material with a constitutive equation of this form will be called a material of 


the rate type. 


It must be noted that (24.12) is not really a complete constitutive equation. 
The stress is not determined by the local motion alone but only when, in addition, 
initial values S)(t,), 7=0,1,..., m—4, for some initial time t,, are given. These 
initial values, on the other hand, should be determined by the history of the local 
motion up to the time ¢,. A constitutive equation of the type (24.12) characterizes 


not a single material but a family of materials depending on m symmetric tensor 
parameters. 


* It was introduced by ZaremBa [15]. 
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In the case of isotropic solids it follows from Theorem 6 that the tensor func- 
tion j in (24.12) must be isotropic and that the constitutive equation reduces to 


(24.13) (Si); Dil; VQ) =0. 
Fluids of the rate type are described by an equation of the form 
(24.14) 9 (Si(t); D.(O; e() =0, 


where g is an isotropic tensor function of m+n-+1 tensor variables and one positive 
scalar variable. However, not every equation of the form (24.14) defines a class 
of fluids of the rate type. It may also define a class of isotropic solids, because 
the stress may depend on the right strain tensor relative to some reference con- 
figuration through the initial values. 

Starting from the general constitutive equations (22.7), (22.9), and (22.11), 
one can easily derive alternate forms for the constitutive equations of the rate 
type: 

In the general case we obtain 


(24.15) {(R(t)? S,() R(); RW? A, () Rd); C() =0 


where A, (t) is the &t® Rivlin-Ericksen tensor (9.2) and where S;(t) is the J stress 
flux defined by 


0) 


(24.16) $0) =F, 07 SE, 


fo=t+ 


We find that Ss may be expressed explicitly in terms of the stress S, its deri- 
vatives up to the order m, and the acceleration gradients £,; up to the order m 
by the formula 


~ I! (q) 
@, Poh HO) canst 
ptqtr=l 
For /=0 we have S,=S, and for /=1 we get the stress flux* 
(24.48) SS =SH=SHETSISE. 
In the case of isotropic solids f is isotropic and (24.15) reduces to 
(24.19) (Sil); 4s); BO) =0*™*. 
For fluids we get 
(24.20) §(S:; Axl); @) =O, 


where § is an isotropic tensor function. 
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* It differs from the flux introduced by Caucuy (cf. the discussion given by 
TRUESDELL in [12], Section 55 bis), in which there are minus signs on the right. 
*k This form was first derived by Cotter & RIVLIN [J]. 
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228 WOLFGANG GERISCH: 
1, Einleitung 

Wir wollen iiber zwei Orthonormalsysteme berichten, die sich besonders zur 
Behandlung von Aufgaben aus der Theorie der an den Randern fest einge- 
spannten, rechteckigen (schiefwinkligen), elastischen, isotropen (anisotropen) 
Platte eignen, weil sie die fiir diesen Problemkreis maBgeblichen Randbedin- 
gungen erfiillen. 

Unsere Ausfiihrungen werden so beschaffen sein, daB sie zunachst bei der 
Anwendung der bekannten numerischen Methoden auf die angedeuteten Probleme 
von Nutzen sind. 

Es wird jedoch noch ein weiterer Zweck verfolgt. Bei der Untersuchung 
gewisser, hierher gehériger Operatoren im Hilbertschen Raum L? (Skalarprodukt 
von der Form ff f(x, y) g(x, vy) dx dy) trat ndmlich die Frage auf, ob durch 
geeignete Orthonormalsysteme ein Ubergang zum Hilbertschen Folgenraum /? 
(Skalarprodukt von der Form 2'2’A,,,, B,,,,) méglich und insbesondere realisierbar 
ist. Die erwahnten Operatoren haben hier die Gestalt von Matrizen (Kern- 
matrizen). Aufgabenstellungen liegen in der Form von unendlichen Gleichungs- 
systemen vor, die abschnittsweise explizit zuganglich sind (Reduktionsverfahren). 
Wir werden zu gegebener Zeit zeigen, daB unser zweites Orthonormalsystem allen 
hier angedeuteten Anforderungen geniigt, wahrend sich das erste System als fiir 
diese Zwecke zunachst ungeeignet erweisen wird. Weiterhin werden wir im An- 
schluB an die eingangs erwahnten Probleme geschlossen darstellbare Kernmatrizen 
angeben. Eine dieser Matrizen wurde bereits zur Untersuchung rechteckiger, 
isotroper, an den Randern durch gleichférmige Schubspannungen beanspruchter 
Platten herangezogen. Wir wollen systematischen Betrachtungen vorgreifen und 
schon in der vorliegenden Arbeit an Hand dieses Beispiels einen Ausblick auf die 
numerische Anwendbarkeit der hier zusammengestellten mathematischen Hilfs- 
mittel geben. 

AbschlieBend weisen wir noch auf einen wichtigen Sachverhalt hin. Ist q,, (x) 


b 

ein volistandiges Orthonormalsystem in L? mit f{ (x) 7(x) dx, dann ist natiirlich 
a b 

¢,(y) ein vollstandiges Orthonormalsystem in L? mit f k(y)/(y) dy. Ferner ist 

Din (X; ¥) =Pm(*) + @,(y) ein vollsténdiges Orthonormalsystem in L? mit 


bb 
J J f(x, y) g(x, y) dx dy*. Man kann sich also bei der Untersuchung der eingangs 
besprochenen Probleme auf die Konstruktion passender, vollstandiger Ortho- 
b 


normalsysteme in L? mit f h(x) i(x) dx beschranken. 


a 


2. Grundlagen 
Es sei ) ein abstrakter Hilbertscher Raum. Zwei Elemente f, ge heiBen 
orthogonal, wenn (/, g) =O ist. Die Elemente heiBen orthonormal, wenn zusitzlich 
l/l =[|g|| =4 ist. Sind die Elemente y,, ys...€ § paarweise orthogonal, so liegt 
ein Orthogonalsystem vor. Dann gilt (Yn> Yn) =4m* Onn Die Elemente ,, 
Yz--.< bilden ein Orthonormalsystem, wenn die Beziehung (Pp, Pn) =O» be- 
steht. Aus nicht orthonormalen, linear unabhangigen Yr», Yo.--€H kann man 
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durch Anwendung des Orthonormalisierungsverfahrens ein Orthonormalsystem 
Y1, P2--.€ konstruieren. Man hat zu diesem Zwecke die Determinanten 
A,,= det ||(y,, Y,)||gn) und B,, zu bilden, wobei B,, aus A,, entsteht, wenn man 
die Elemente der m-ten Zeile durch y,, y,... ersetzt. Die q,, y...€ § stehen 
dann mit den Wis Ws re 5) wie folgt in Verbindung: 9,=y,- Aj}; 9,,=B5,,: 
(A,,-1:4,,) 2, ™ =2,3.... Wir werden dieses Verfahren im nachsten Abschnitt 
auf einen konkreten Fall anwenden. 

Ein abzahlbares nicht notwendig orthonormales System von Elementen 
XY %2--.©€H heiBt in H volistandig, wenn fiir jedes /€ eine Beziehung der Form 


M(e) 
/-> Cm Xm 
1 


x», »1m Mittel™ beliebig gut approximiert werden kann. Man schreibt dann 


<e besteht; d.h. wenn jedes f€ § durch Linearkombinationen der 


f=>iCm %m- Bilden die y,, ein Orthonormalsystem, so heiBen die c,, Fourier- 
1 


Koeffizienten von / beziiglich des Orthonormalsystems y,,. Infolge der Ortho- 
normalitat ergibt sich dann c,,=(f, v,,). Das System der y,, heiBt abgeschlossen, 
wenn es kein vom Nullelement verschiedenes g € § gibt, das zu allen y,, orthogonal 
ist. Aus Vollstandigkeit folgt Abgeschlossenheit und umgekehrt. Die beiden 
Eigenschaften tibertragen sich im Falle nicht orthonormaler y,, durch den Ortho- 
normalisierungsprozeB auf das zu konstruierende Orthonormalsystem. Wegen 
des hierher gehérigen Begriffes ,,dicht“ verweisen wir auf [2] *. 

Gelten alle Axiome von-{ bereits fiir eine Teilmenge §,¢ , so heiBt , 
Unterraum von {. Die Differenz $, der Mengen und §, ist auch Unterraum 
von §. wird als orthogonale Summe wie folgt geschrieben: ) =), 9% Dp. 
Jedes Element /€ § ist eindeutig zerlegbar in f=g +h; gE; hE Do. 


3. Das Funktionensystem ,,,(~) =sinax-sinmsx 


3.1. Vorbemerkung 


Wir stellen einleitend fest, daB das genannte Funktionensystem in nicht 
orthonormiertem Zustand bereits von M. M. FILONENKO-BoropiTscu [3] und von 
E. METZMEIER [4] zur Untersuchung von Aufgaben aus der Festigkeitslehre be- 
nutzt worden ist. Auf diese Arbeiten werden wir an anderer Stelle noch aus- 
fiihrlicher zu sprechen kommen. 


3.2. Anwendung des Orthonormalisterungsverfahrens auf das Funktionensystem 
7 
Y», (x) =sinx-sinmx im Raume L? mit f h(x) (x) dx 
0 
Die unter 2. angefiihrten Formeln liefern im vorliegenden Falle das folgende 


Orthonormalsystem 


™m 


/ 8 : 1 : 
(62) = | SS 56 . x m ungerade, 4 
Pin (X) \ CLEONG sin x Or. sin | «| m ung e (1) 
o m 
; es 7 ae see Es % ms é A 
Grlk) = Va eye aA 2 iH sin | x} m gerade. (2) 
* [2] 5.9 
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Durch Anwendung der Additionstheoreme erhalt man: 


a m—1 
2 1 | oe 
, (x) =|/—-= -(4+2-> cosux—m-cos(m-+1) x| m ungerade, (3) 
Pm (%) /2 |/m (m-+-2) px i ( 
m—1 
2: > cossux — m-cos(m +1) x 
Le Oiaen 


m gerade. (4) 


, / - 
Pin (x) al \ IU |m (m i 2) 


Fiir (1)—(4) geben wir die folgenden geschlossenen Darstellungen: 


ps 1 sin m % e HH) 
P(x) = \2 ; Ym (m+ 2) 4 Site sea! a0 (5) 
Oya \2 . TES -{sin(m +1) x -ctg x — (m + 1)-cos(m + 1) ee (6) 


3.3. Die Darbouxsche Differentialgleichung 
Durch formale Rechnung beweist man, daB (1)—(6) Lésungen der Differen- 
tialgleichung 


i eae | eon Be 3 x 

Pm (x) | sin? x (mm si 1) if Pn (X) 0 (7) 
sind. (7) ist eine spezielle Darbouxsche Differentialgleichung. Differential- 
gleichungen dieser Art sind von W.v. KopPpENFELS {[5| untersucht worden. 
Daselbst werden fiir (7) mit Hilfe Jacobischer Polynome G,,(f, g, x) unter anderem 
folgende Losungen angegeben: 


Pig (4) SI 0? Gy (2, psi) m ungerade, (8) 
Pm(x) ~ sin? x - cos x-Gy—2 (3,8,sin?x)  m gerade. (9) 


Unter Benutzung der Definitionsgleichung fiir die G,,(p, g, x)* wollen wir noch 
andeuten, wie die von uns gefundenen, handlichen Darstellungen (3) und (4) aus 
(8) und (9) ermittelt werden kénnen. Aus (8) ergibt sich: 


Q(x) ~ sin? x - Go(2, she sin? x) = . (10s 2%), (10) 
@3(x) ~ sin? x - G,(2, : ‘ sin? x == SNe 4 < [1 - eae sin? x) 
\ < , 2 (11) 
eke ib (1 -+ 2cos2% — 3cos4x) usw. 
Entsprechend erhalt man aus (9) : 
P2(%) ~ sin® x - cos x - Go (3, 3, sin? x) = sin? x- cos x-4 
(12) 


= 4: (2cos2cos x — 2cos3x) usw. 


Fin allgemeiner Beweis zu (10)—(12) z.B. durch den InduktionsschluB ist mit 
einem auBergewohnlichen Aufwand an trigonometrischen Umformungen ver- 
bunden. Wir haben daher von diesbeziiglichen Uberlegungen Abstand genommen. 


SAIGON O07 
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a 


3.4. Abgeschlossenheit (V ollstandigkeit) 
Behauptung. Das durch (1)—(6) dargestellte at aa P(x) ist 


abgeschlossen (vollstandig) beziiglich des Raumes L?2 mit bi h(x) t(x) dx, 


Beweis 1. Wir zerlegen Te mit: f h(x)7(x)dx wie folgt in zwei Unterrdume: 


L217 @ 13 *- Dabei ist f(x) € rome fast tiberall f(x) =f (a — x), und f(x) €L3, 
wenn fast iiberall f(«) = Ceo Dann weisen wir nach, daB (1) beziiglich 
Li und (2) beziiglich L3 abgeschlossen sind. Zunachst wird jedoch angenommen, 
daB (1) beziiglich Li und (2) beziiglich L3 nicht abgeschlossen sind. Demnach 
existieren in L{ und L} Funktionen g,(x) und g,(x), die zu allen Gliedern von 
(1) bzw. (2) orthogonal sind. Fiir die genannten g (x) setzen wir mit den in ti 


bzw. L} abgeschlossenen Orthonormalsystemen [ites \2 COS M%, W=254, 6... 


COS PH K, IW ANS 25 25, wie folet, an waga(% ) ~|/2. (“ +m cos me 
£9 (x) ~ me Diem cosmx. Nach dem oben Gesagten ibe See dispute 


J ml rales ieee m=1,3,5. 1S ml) %) d4—=0, m—=2, 4, 6 .. gelten. 


Wir zeigen in einem Hilfssatz, daB aus den zuletzt aufgefiihrten Beziehungen 
auf das Verschwinden der Fourier-Koeffizienten c,, geschlossen werden kann. 
Die gemachten Ansatze fiir g,(x) und g.(x) liefern dann g,(x) ~O und gy(x) ~0, 
und wir erhalten einen Widerspruch zu der eingangs gemachten Annahme iiber 


die Nichtabgeschlossenheit von (1) bzw. (2). Damit ist die Behauptung bewiesen. 


Hilfssatz zu Beweis 1. Sind a,, und b,, Fourier-Koeffizienten zweier Funk- 
tionen h(x), 7(x) € L? beziiglich eines in L? abgeschlossenen ee 


so besteht die verallgemeinerte Parsevalsche Gleichung Lae 4) ai pa (rbedlite 
Wir wenden diese ae zur Umformung der petamineen ae £1 (% ) Ci), 
MVNA, 5”. ath On (X) gaa e 0) 12,4, 6 veswans Debes beriicksichtigen 


wir, daB (1) und (2 (2) auch als durch die erwahnten trigonometrischen Orthonormal- 
systeme dargestellt aufgefaBt werden kénnen. Man erhalt die Gleichungssysteme 


m—1 


(m (m+ 2))~# (2 Ex. 2 ON mM - bee 0, mM ==41,3,5%..; (m(m+2))-4- 
m—1L 

2,» Cy M - Cn a0 oe 4,6.... Diese Gleichungssysteme liefern der 
hohe ue : : 

Reihe nach c= |/2 Cy, Cg=|/2 Co---; Cg=C1, Cs =%--.. Man gewinnt nach abzahl- 

bar vielen Schritten die geradzahligen Fourier-Koeffizienten als Funktion von 

cy und die ungeradzahligen als Funktion von ¢: ¢g=cy=ce=---=|/2 9; 

C4=C,=C,=+:» =c,. Wir hatten bei diesen Uberlegungen vorausgesetzt : g, (x) EL}; 


&o(x) CL3. Entwickelt man solche Funktionen, so muB auf alle Falle die Bessel- 


sche ee » (igs </ei gi(x) dx 00s ce fe g3(x) dx<oo erfiillt werden. 
il, ‘ae 


Das ist eae nur uGalicHS went cg=c,=0 poset wird. Dann ist aber c,,=0, 
m=O0,1,2... wie zu zeigen war. 


* [7] Sy 2 
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Beweis 2. Wir zeigen zunachst, daf fiir die Vollstandigkeit von zusammen- 
i ee Funktionensystemen der Form y,,(«) =/(%) + %m(*) beztiglich L? mit 
f h(x)i(x)dx die folgenden ca: Bedingungen hinreichend sind: a) y,,(%) ist 


i beziiglich L? mit Sai h(x)i(x\dx; b) f(x) und y,,(x) sind stetig in 
Se 4 ANATS) Jt ; besitzt Agesloct hdchstens hebbare Unstetigkeiten. Mit Riick- 


sicht auf einen Satz von F. G. Tricom1* und wegen Bedingung a) kénnen wir 
uns darauf beschranken, die Vollstandigkeit von y,,(x) beziiglich y,,(«) nachzu- 
weisen. Da das Produkt stetiger Funktionen stetig ist, kann aus den Bedingungen 


b) und c) gefolgert werden, daB /?(x) und dom) ° tn(*) in axx<b stetig sind. 


P(x) ; 
Der zuletzt genannte Ausdruck hifert fiir m =n: "th ax= Hin (#) AxX< 0; 
R f?(*) f(x) 
nt € L* mit Fh) (x) dx. Unter Zahiitedabiic der pueee a) kO6nnen wir 
y M 
xm () dmieren: [ (#=()_y¥' 
ie wie folgt ne (x) 2 ten) Kile )) dx< 5 7 “_ (M Schranke 


von /?(x)). Der Integrand ist eine stetige Funktion, da auch die endliche Summe 
stetiger Funktionen stetig ist. y,,() ist genau ese pies Seg beziiglich y,, (x), 


wenn 7,,(%) wie folgt approximiert werden kann: pean -¥ Dum Pu l% (9) ne 


) a 


: Mn 2 
| (im) —f (x) +2) Bum: tu(2)) dx<e. Wir schreiben verte Ausdruck um und 
a w=1 3 


schatzen unter Bezugnahme auf unsere obigen Stetigkeitsiiberlegungen mit dem 
erweiterten Mittelwertsatz der Integralrechnung wie folgt ab: 
b 


fre (2 = Sua w)) ds 


a 


fe Mi, 4 Mn, 
< 72 ° Xm(*) -_ b < ; : 2 VM. ieee pod 2 s 
sre. { (208 Dom tale)] da <M. f dom y,(x)) ax 
Wie oben erwahnt kann der zuletzt genannte Ausdruck durch passende Wahl 
der b,, kleiner ¢/M gemacht werden. Damit ist der erste Teil dieses Beweises 
beendet. 


Offenbar erfiillt unser Funktionensystem y,,(«) =sin x - sin es die Bedin- 
gungen a)—c) in 0S *< a und ist demnach vollstandig in L? mit f h(x) 1(x) dx. 


Die Vollstandigkeit iibertrégt sich durch den Orthonomiolisierunesprozes auf 
(1), (2). 

Zusatz. Das genannte Funktionensystem kann auch als Differenz zweier 
cosinus geschrieben werden. Im Anschlu8 an diesen Sachverhalt hat V. STEPA- 
NoFF [9] notwendige und hinreichende Bedingungen fiir die Vollstandigkeit von 
Funktionensystemen aufgestellt, deren einzelne Glieder Linearkombinationen 
eines volistandigen Orthonormalsystems sind. 


* [8] S. 27 
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3.5. Die Integrale Fonts) x) p, (x) dx und Soh x) Dp, (x) dx 
Man erhalt: 


J ants): (id 0 (m—s)=1 (mod 2), (13) 


7% 


ee (m —s) =1 (mod 2), (14) 


Pony m(m +2) 
Pine) 9.2) de = (owe) pe 


Mas (m—S) =0 ~*~ (mod 2), 


“MH 4) Oy, 5 


Bil 
3 (15) 


f Pm (%) yp; (x) dx=(m-+1)- ee . Joi (3 m?+ 6m-+1) +(m-+1)8-6 


F s(s+2) 15 ™*) (16) 
MSs (m—s)=0 (mod2), 

Sant w) dx = fo (x ) Pm (2) 4%, (17) 

Soke x) dx = fg %) Pm (x) dx. (18) 


4. Die zugeordneten Kugelfunktionen II‘, (x) 
4.1. Definition 


Die in —1X*X1 auf 1 normierten Losungen des Eigenwertproblems 


(1 — x2) - f(x) — 2%- f(x) + {m(m are as f(x) =0, 
m=4,5...; f(—1),#(+1) endlich 


(19) 


heiBen zugeordnete Kugelfunktionen /7},(x). Verzichtet man auf eine Nor- 
mierung, so entstehen Lésungen von (19), die unter dem Namen zugeordnete 
Legendresche Polynome bekannt sind. Wir bezeichnen diese mit P(x). Die 
IT4 (x) stehen mit den P(x) wie folgt in Beziehung: 


4 RE \n. ee Lm=3 5 
LT) Ae (x) | @ | S) Fond sy m4. (20) 
Die zugeordneten Legendreschen Polynome P(x) lassen sich durch Legendresche 
Polynome P,,(x) ausdriicken: 
4 qdmr4 
PA(x) = (1 — 298 Sim = yt 


a“? — 4)", (21) 


4.2. Orthonormalitat — Abgeschlossenheat 


Die eno uceen IT} (x) bilden ein abgeschlossenes Orthonormalsystem beziig- 
lich L? mit zs AY x) i(x)dx. Da sie definitionsgemaB den Pi(x) proportional sind, 


braucht man nur Orthogonalitat und Abgeschlossenheit fiir das System der Bj (x) 
nachzuweisen. GemaB den unter 4.1. gemachten Bemerkungen iiber Normierung 
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erwarten wir fiir die P4(x) die folgende Orthogonalitatsrelation : 
sen 
f PA(x) BS (x) dx = 
an 
Eine ausfiihrliche Herleitung findet sich bei J. LENSE*. Bei R. Courant und 
D. Hirpert** werden mit den Methoden der Variationsrechnung allgemein 
gehaltene Aussagen tiber Vollstandigkeit von Systemen von Eigenfunktionen 
gemacht. Das hier zur Debatte stehende Eigenwertproblem *** ist so beschaffen, 
daB aus den erwahnten Sdtzen auf die Vollstandigkeit der P(x) beziiglich des 


Raumes der in —1<%<1 stetigen Funktionen geschlossen werden kann. Dieser 
tar 
Funktionenraum liegt in L? mit f h(x) 7(x)dx dicht. Das hat zur Folge, daB 


a 


2 LP(m+5) 
2m+1 T'(m— =3y' Orga (22) 


die P4(x) auch beziiglich L? mit ‘° Pe: (x) dx ein vollstaéndiges (abgeschlossenes) 
Orthogonalsystem bilden. ot 
ails ” 
4.3. Die Integrale f{ II4(x) 14 (x) dx, fu (x) dx und [ko ea bel Cattle 
= 


Man hat sich bei der Bestimmung er aufgefiihrten Terenas hauptsachlich 
auf die Formeln 


de? ] +1 rT’ —- , ”r 
Pax) = Tet Bal) —8- TIES Bh ale) +24 B(x), (23) 
qi T'(2m+n+2) 


m=>Q (Trinity 1907), (24) 


dgmet tmen()! = Seto £11 Pn) 
[PO(x)-P(x)dx=0 fir 1>0, m<s (25) 


zu stiitzen und kommt nach umfangreichen formalen Rechnungen zu folgenden 
Resultaten: 


+1 , 
[Th (*) I(x) dx =0 — (m—s)=0 (mod 2), 120) 
=a 
rw, tlio, 
J HE (x) LIS 0) dx = — [ITS (x) Dp (x) ax, (27) 
= —1 
+1 , ros 
Th Oe ar 1 OP CaS) s=3) 
ai (x) 1S (x) dx -\(m+4 (s Teac P(m—3) | P(s+5) (28) 
MSs (m—s)=1- (mod 2), 
ye! Ay 
SAR) LTE (at) diac) (m — s)=1 (mod 2) (29) 
ao ahys 
LAR) LIM &) deme 0. i = s\ee ts (mods (30) 
+1 ” Se hi 
J Hn (3) ALS (x) dx = J TTS (x) 1p (2) dx (34) 
eae = 
* [10] S.151 


“$8 Ai]| Sy SKole: 
*kK 61S. 443 
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+1 Iv 
J 174 (x) - I74( () dx = fT) TT4 (x) d (32) 
oH ” 1 1 1\ P(m+s) I(s—3) 
4 ame 
NPC TONS ces de Toroeee 
x {5s(s +4) —9m(m +1) + 60} (33) 
mss (m—s)=0 (mod2), 
+1 Iv 


[ITS (0) TTS) dx = py (m+ 4) (04+ 4) FOES ERD 


x {7(s — 2) (s +3) - ((m — 2) (m— ei 
— 5(m — 4) (m +5) - ((m — 5) m(m +4) (m + 6) + 252)} 
M=s (m—s)=0 ree 


(34) 


5. Anhang: Numerisches tiber freischwingende, rechteckige, isotrope, 
an den Randern durch gleichférmige Schubspannungen beanspruchte Platten 
5.1. Abkiirzungen 


Durchbiegung als Funktion von x, y, 


= aS Plattensteifigkeit, 


Elastizitatsmodul, 
Poissonsche Zahl, 
Plattendicke, 
Dichte, 

2af Kreisfrequenz, 


SS os Se) les les 


o 
—— 


Seitenlangen, 


Schubkraft am Rande pro Langeneinheit. 


aq 


0.2. Differentialgleichung — Randbedingungen 
Biegeschwingungen rechteckiger, isotroper Platten, die an den Randern durch 
gleichférmige Schubspannungen beansprucht werden, lassen sich durch die Dif- 
ferentialgleichung 


Ot w a\2 dw a\to4w-, 2% a [a\* Pw QW aE Se pet 
+2-(4)-soeat($) gatas (5) 5: (F)-o% w=0 635) 


aa Ox? Oy oy QO 
beschreiben. Zu (35) treten in unserem Falle die Randbedingungen fiir Ein- 
spannung 
i) = 08 ia ir gee cele  @) == 0). eS) firs. y= (36) 
Ox Oy 
Es empfiehlt sich, folgende Abkiirzungen einzufiihren: 
2 el a\s 2. _ 2t (a\?, es: 
A= &§~.(5} CS Gr Dae Oe (37) 

(35) geht dann itiber in 

otw oD 4. ow aa ee aa 

ae” aa ay? “pe By +O % re Aw =0. (38) 
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5.3. Operator — Matrix 
In der Differentialgleichung (38) ist der Differentialausdruck 


ot My fe —- a. tg tate (39) 
ax c aE 2 oy4 Ox Oy 


maBgeblich. Mit Hilfe von (39) und (36) laBt sich im Hilbertschen Raum p Be 
+ +41 

mit dem Skalarprodukt f f f(x, y) g(x, y)d«dy ein Differentialoperator defi- 
al il 


nieren, den wir mit A bezeichnen wollen. Neben Eigenschaften wie Hermitizitat, 
Halbbeschranktheit und Hypermaximalitat bleibt zu beweisen, daB A fiir #=+0 
durch zwei Unterréume reduziert wird. 

Der Ubergang von L? nach /? (Skalarprodukt von der Form 22A,,,,B,,n) 
ist mit der Schrédingerschen Formel formal sehr einfach zu vollziehen, wenn 
man in L? ein fiir A zuldssiges vollstandiges Orthonormalsystem W,,,, (x, y) 
kennt. Dem Operator A in L? entspricht in /? die Matrix 


+1 41 
ea ee Widxdy. (40) 
—1 


Setzen wir voraus, daB unser Orthonormalsystem ®©,,,,,(x, y) =IT4(x) - ITA (y) 
fiir A zulassig ist, so fiihrt (40) auf die Integrale (26) —(34). Mit den Abkiirzungen 


of Adis Bear! 
turn Ir 2) (8+ 3 a) a Tine) Fes) (44) 
asm m>s 
0) (m—s)=1 (mod 2) 
7 (s—2) (+3) x 
x ((m—2) (m—A) (m+-3) (m+5) +60) 
bo, .= 
ms sGnsea Ones mss (42) 
x ((m—5) m(m-+1) (m+6) “+ 252) (m—s)=0 (mod 2), 
Os m>s 
0 (m—s)=1 (mod 2) 
ons = 5s(s +1) —9m(m +1) +60 mss (43) 
, (m—s)=0 (mod 2), 
sm m>s 
lt mM<s = 
oe mee (m—s)=1 (mod 2) (44) 


me) (m—s)=0 (mod 2), 


erhalt man folgendes Ergebnis: 


1 
Ann st = 517 Bn 5° De On: PP 20- 


1 
315 Ams oms 315 me 
HO Ons: ro Ant Oye FOO: Ams Ant t4ins Uys M,S,N,tZ4. 


Pee “S7z9C SL‘6S 4 
aie $79'% J || sesoe- 
SL‘6S 4 
oot ae @- S64 + Sz 1E0‘0— 2 Sz 1€0'0— 
1 € 25 ie 
(Zr) @-SEOASTIEO'O G-LE-LPEASTZ 900'0— ak SZ ¥£0'0 
@-s9s 4 — @-SLEIZ'T a 54 £60'0 
oe @-S8SA + Sz FE0'0— @-SL£60'0— 
‘ ‘ $ ‘ 
@-S8SA°Sz1E0'O— @-SLE68'0— 6+ Ez /L' $4 EvE'O 
@-SL£60'0— Jl secveo @: SL EpE‘O— 
aie Sth'9 $L£6'8— 
Sz4'9 epee SL'Sph 
SAS SL‘Sbh SZ 89v'EzE 
S€6A + Sz 1EO‘O— = -$2 €60‘0— : ‘$t F£0'0— 
@-SE6A SZ FEO‘O ey 0) “Ter Z4£0'0 
, 4 S 
(9r) ¢ “SLE6D'!0—_-@ - SEA“ $Z 1EO'O— <I. -$z 1£0°0— 
@-SLEITE— SL BIE @-SL£68'0 
é s , s 
@-S8SA-S7IEO'O G+ ET /} SL £60'0— “ey /L $4 £ef'0— 
, , S 
¢: s *$2£60'0— @- S854 -SzLEO'O "ey /L 54 £rf0— 
@ SZ £600 @- SL £600 @-SLEPE'O 


=| ‘ BL | eee 
Bi/L $40 ol Seek 


eos “th bez 


ese ZSE 
e.g a &-S£6A*SZ FEO'O CASSS (Sa Ue arp @- SBS) + SZ FEO‘O— @-SL£60°0— 
S 
@-S6A-STIEEO'O—  G-LE-EbbA-SZ900'0— G-SLEITE &-SESA+ Sz HEO‘O— @-SL€68'0— rer “SLEPE‘O 
a *S$Z bE0'O as “StHeo'O)— G = -$2£60° -$L€60'0— ASSETS -SLEpEo— 
OE) CVE = ray z1E0 Oey 0) ¢ Grey v£‘O @ 
$$ -£4 Ese Haniee end Ny 9) 
Se ae Aull 8 ail SLE‘0 61/L 2 SPt 6h ¢ £1 scc60% $Z £600 
: th $s 
ao 4g $2 898'9L08 Sa SKA! Pall |e aie *S6‘081 SoS 


SLE9‘O£ — 


SLEEE ooh SMe 


€h : $f. | [eenaee ae ‘ ‘ ‘ 
ll -sisn mit l all SA? SL 8tL pee all $z9't $z90'S 
eae eae: 5 L£9'0f — IB ‘ 2a: 
ai ize 08h 9'0 rm tS ee (0) 
? $$ : i te S e 
$££60'0 Ap Ver 1/L' S292 $z90'S 7 “$2°0 SL g9%‘€9 
é $$ ‘ Peer ( “ $ ; ‘ 
@-SLEA“STHEO'O— #77 /[ $2 t60'0— @ SLE I — @-SBSA-STHEO'O E+ ey /| SL £60'0— @-S2£60'0 
ae S£€60'0O—  # - SEGA SZ HEO‘O— @-SL BIZ eal: $£ €60'0— @+S8S/ -Sz FOO @- SLE 600 
ru StIeo‘o— all Sz E0‘0— r-SLE6Z'O ¢- ae SLEPe‘o— o-g/|-seereo- @-SLEPEO 
; $$ $$ F ab sso ae 
Se tel ge 7 S78°9 i “SOL raze Bh SL 896° 
$s ss ; Ue ieee an SS ieee 
ae a Sa $Z3'9 seal POe Vez “SL°ESL HE A-sos 
$$ ; $$ ( : 8. : a3 : Q 
Gi “$789 Zh Sz3'9 SL 867 S| $s‘8 al $L86'91 — S76‘7 
Hs: aaa il es: SL gtz‘o£ $z6! -$£g9'h— 
127 a 127, io eh 8 weet Coe 7 89" 
ath. gt Hl esc'egs $1869) — Sa 6°65 ff sy: 
vez) L' 9° [rad 8 i BO COE BID uO. et 9 
Sy lee eae 
; *SL896'b— Lh S‘OL $z6'% 35299. “$Z°9 SZ 896°LE 
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Mit A wird auch a,,,., fiir ®+:0 durch zwei Unterraéume reduziert. Den im 
Unterraum 1 liegenden Teil von 4,,,,,,ethalt man durch die Bedingung (m — n)=0; 
(s —t) =0 (mod 2). Der im Unterraum 2 liegende Teil ergibt sich analog durch 
(m—n) =1; (s—t) =1 (mod 2). 


5.4, Eigenwerte 


Wir wollen hier fiir eine Platte vom Seitenverhaltnis 2:1 jeweils den ersten 
Eigenwert der beiden Teile von 4,,,,; als Funktion des Schubparameters # 
bestimmen. Es zeigte sich, da8 man brauchbare Naherungen erhalt, wenn man 

sich auf die Betrachtung des neunten Ab- 

schnitts .des jeweiligen Teiles von 4@,,,.1 
beschrankt. Wie die Gleichungen (46) und 
(47) zeigen, lassen sich die Koeffizienten 
dieser Abschnitte aus (41)—(45) verhaltnis- 
maBig einfach und ohne Benutzung von 
Rechenmaschinen ermitteln. Aus (46) und 
(47) zu bildende Determinanten sind weiter- 
hin zu vereinfachen. Man kann namlich durch 
passende Erweiterung samtliche daselbst auf- 
tretenden Wurzeln beseitigen. Zur Auswertung 
wurde das Eliminationsverfahren angewandt, 
wobei zur Vereinfachung der Rechnungen die 
eingangs formulierte Problemstellung umge- 


0 30 100 i? 200.8 20 kehrt wurde. Wir geben eine Ubersicht iiber 
Abb. 1. Abschnittseigenwerte ,2((6) und ,a“(9) die erhaltenen Resultate in Tabelle1 und 
Abb. 14. 


Abb. 1 zeigt, daB fiir d< 196 ,A(9) (9) <4 (8) ist. Fiir > 196 ist ,A!9 (9) > 2! (8). 
Bei } = 196 iiberschneiden sich die Kurven ,4!) (9) und ,A{® (9). Es steht zu ver- 
muten, daB der erste Eigenwert von a,,,,,; fiir ®~196 die Vielfachheit zwei hat. 


Tabelle 1 
o 0 A18:22 170,90 206,50 
(6) | 592 480 240 | 0 
o 0) 119,79 | 166,20 | 204,69 
ie () as 1008 640 320 0 


Von physikalischem Interesse ist derjenige Wert # =®,,i,, fiir den der kleinste 
Eigenwert von a,,,,, Zu Null wird. Wir beschranken uns wiederum auf die Be- 
trachtung von ,qj7,,, und a! , und stellen fest (Tabelle 1; Abb. 1), daB 5A (9) 
der maBgebliche Eigenwert ist. Zu A) (8) =0 liefert uns Tabelle 1 diesen sog. 
Beulwert fiir die von uns untersuchte Platte vom Seitenverhaltnis 2:1 naherungs- 
weise ZU oAii¢= 204,69. Al (8) wird zu Null fiir = 206,50. Dieser Wert ist 


weniger wichtig, da d,.i¢<,0 ist und die untersuchte Platte bereits bei RU 
ausbeult. 
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Beulwerte fiir Platten vom Seitenverhaltnis 2:1 wurden von S. Icucut [11], 
W. Mouerr [12] und B. Bupransxy und R. W. Connor [13] angegeben. S.IGucHI 
benutzt Reihenansaitze, W.Moueir geht mit dem Differenzenverfahren vor, 
B. BupIAnsky und R.W. Connor wenden die sog. Lagrangesche Multiplikator- 
methode an. In Tabelle 2 vergleichen wir die Resultate der einzelnen Autoren 
mit den von uns gefundenen nach (37) passend umgerechneten Werten. 


Tabelle 2. Ndherungsweise bestimmte Beulwerte fiir eine Platte vom Seitenverhdltnis 2:1 


S. IeucHi W. Mone1rT ee Ls a ae Eigene Werte 
; 10,96 : | 
vy) Bears. _ (0) 10,¢ 
: (extrapoliert) | ny O48 
11,80 | 
, | 
2 rit = 10,42 | 10,34 | 10,37 
(extrapoliert) | 
echen- : ; a ¢ 
Gleichungen 18 Gleichungen 9/11 Gleichungen 9 Gleichungen 
aufwand ? so} : beets ? i ee a 8 


5.9. Ergenelemente 


Wir beschraénken uns auch hier auf die Betrachtung von ,a!) , und .a\), |. 


Fir 0 = 206,50 und 4,4, = 204,69 erhalt man aus (46) und (47) die in Tabelle 3 
zusammengestellten, auf 1 normierten Abschnittseigenelemente in /?. Diesen 
Abschnittseigenelementen entsprechen in L? die Elemente 


if (%, y) A Din (%, y) und af (%, y) =22 oA. Din (%, y). 


mn mn 


Bei der Auswertung stand die Tabellierung [14] der zugeordneten Legendreschen 
Polynome zur Verfiigung. Die Abb. 2 und 3 enthalten das Ergebnis. 


Tabelle 3 


9) (9) (9) (9) 4 (9) 
Aus 1446 A Ga 14 66 144g 


0,003 3623 


| 

ie 0,655 13 —0,048585 | —0,59071 | 0,095503 | 

fiir 

Fe | 9 

18 = 206,50 AY) Ae | Aw | A® | 
= — as | ona —_ 

| | | | 

0,076551 | —0,41844 | 0,018718 |  0,17048 


(9) (9) (9) (9) 
ose 2456 2476 25s 


0,004 979 


ae 0,885 38 —0,081199 | —0,18517 0,037922 | 

fiir | | 

2 rit = 204,69 Age 2452 2A 
ss a= (= ee = | = 

0,010894 | 0,24118 | = 0,33966 | —0,010534. | 


Beulfiguren fiir Platten vom Seitenverhaltnis 2:1 finden sich bei S. Icucui{ 13 | 
(und [76]) und bei W. Mouerr [12]. Wir geben die beiden maBgeblichen Figuren 
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in Abb. 4 wieder und stellen fest, daB die von S. Icucut publizierte Figur bis auf 
zwei Knotenlinien mit unserer Figur aus Abb. 2 iibereinstimmt. Die von 
W. Moueir angegebene Figur wird bis auf Feinheiten durch unsere Figur aus 


=—— 
= S = t a 
| | L] 
a = 


Abb. 2. Naherungsweise bestimmte Eigenfunktion. Seitenverhaltnis 2:1; ,# = 206,50; ~AM=o0 


Abb. 3 bestatigt. Wie bereits bei der Bestimmung der zugehérigen Werte von @ 
angedeutet worden war, ist unsere zu ,# gehdrige Figur (Abb. 2) nicht stabil, 
da > oti, ist. Dieser Sachverhalt ist auch fiir die von S. IcucHI berechnete 


Abb. 3. Naherungsweise bestimmte Eigenfunktion, Seitenverhaltnis 2:13 .Pkrit= 204,69; AMO (Beulfigur) 


Figur zutreffend, Stabil ist die zu .0,,i, gehdrige Figur (Abb. 3). DaB die Platte 
tatsachlich das Bestreben hat nach dieser Figur auszubeulen, wird von W. MoHEIt 
auBerdem noch durch Versuche erhartet. In den Arbeiten [/1]*, EW etaral le (ee 


* [11] S. 8,9 
** [15] S. 80 
*** [16] S. 820, 823 
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und [/2|* finden sich gewisse diesbeziigliche Unklarheiten. Wir hoffen, daB 
wir insbesondere durch Abb. 1 und die dort gemachten Bemerkungen zur Auf- 
klarung der Verhaltnisse beigetragen haben. 


CR 


Abb. 4. Naherungsweise bestimmte Eigenfunktionen. Seitenverhaltnis 2:1. (Zum Vergleich spiegelbildliche Wiedergabe; 
reziprokes Seitenverhaltnis bzw, andere Richtung der Schubkrafte) S. Icucnr ,0’~ 10,96; W. Monerr yd .3¢= 11,80 


5.6. SchluBbemerkung 
Das gegebene Beispiel zeigt, wie nutzbringend sich die in der Einleitung 


angedeuteten Gedankengange anwenden lassen, wenn man sich auf die gegebenen 
Hilfsmittel aus der Theorie der Orthonormalsysteme stiitzt. Wegen des sehr 
einfachen Aufbaues der Kernmatrix (45) lieBe sich die Anzahl der fiir die Nahe- 
rungsrechnung herangezogenen Gleichungen ohne Schwierigkeiten erhéhen. Zur 
Auflésung kénnte man Rechenautomaten verwenden. Fehlerabschatzungen sind 
moglich. 
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Uber den Einfluf der Stromlinienkriimmung 


auf die Stabilitdt laminarer Strémungen 


HERMANN WITTING 


Vorgelegt von H. GORTLER 


Einleitende Ubersicht 


Betrachtet man die Strémung einer zahen Fliissigkeit, die zwischen zwei 
koaxialen Zylindern bei einer konstanten Rotation des inneren Zylinders ent- 
steht, so beobachtet man oberhalb einer gewissen Drehzahl am auBeren Zylinder 


die Bildung charakteristischer Streifen (vgl. die 
in Abb. 1 wiedergegebene Versuchsapparatur). 
G.I. TayLor [/] hat das Auftreten dieser Streifen 
durch das Vorhandensein von Wirbeln erklart. 
Unter der Annahme von Zylinderradien, die groB 
sind im Verhaltnis zu ihrer Differenz, konnte er 
namlich zeigen, daB die laminare Drehstromung 
zwischen den Zylindern oberhalb einer kritischen 
Drehzahl instabil ist gegeniiber Wirbeln, deren 
Achsen in die Umfangsrichtung zeigen. Die Breite 
der Streifen entspricht dann gerade dem Abstand 
dieser abwechselnd rechts und links drehenden 
Wirbel. 

Spater hat dann H. GOrTLER [2] fiir einen 
ahnlich gelagerten Fall, ndmlich die Grenzschicht- 
str6mung an einer konkaven Wand, die Instabili- 
tat gegeniiber Langswirbeln, wie sie in Abb. 2 
dargestellt sind, theoretisch nachgewiesen. Diese 
Grenzschichtstromung laBt sich als eine Taylor- 
Stroémung auffassen, deren innerer rotierender 
Zylinder durch eine ,,auBere Str6mung™ ersetzt ist. 
Aus methodischen Griinden mu8B H. GORTLER [ 2) 
dabei die Voraussetzungen machen, daB die Wand 


Abb. 1. Versuchsapparatur mit Taylor- 
Gortler-Wirbeln (nach F. Scuuttz-Grunow 
und H. Hern [9]) 


und alle Stromlinien in der Umgebung der betrachteten Stelle konstante Kriim- 
mungen haben und daB der Kriimmungsradius der Wand groB ist im Verhaltnis 
zur Grenzschichtdicke. Als Instabilitatsbedingungen ergab sich dabei, da der 
Parameter Re, |/9/R (Rey die mit der Impulsverlustdicke # gebildete Rey- 
noldssche Zahl, R der Kriimmungsradius der Wand) gréBer sein mu als ein 
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gewisser, von der Form des Grenzschichtprofils praktisch unabhangiger Wert. 
Derartige Wirbel wurden dann auch durch H.W. LIEPMANN [3], [4] experimen- 
tell nachgewiesen. 

Eine physikalisch ahnliche Situation liegt bei einer Grenzschichtstromung 
immer dann vor, wenn die Stromlinien konkav sind zur Seite wachsender Werte 
der Geschwindigkeit. H. GOrTLER [5], [6] hat deshalb die Vermutung ausge- 
sprochen, daB auch iiberall dort Instabilitat gegentiber solchen wirbelartigen 
Stérungen vorliegt, wenn nur eine weitere, der obigen entsprechende Instabili- 
tatsbedingung erfiillt ist. In der vorliegenden Arbeit soll diese Frage naher 
untersucht werden fiir den Fall, daB es sich bei den konkav gekriimmten Strom- 
linien um die Taler einer wellenférmig gestérten Grenzschichtstro6mung handelt. 

Diese Frage ist z.B. im Hinblick 
auf den laminar-turbulenten Um- 
schlag von groBem Interesse. Be- 
kanntlich hatte W. TOLLMIEN [7 | 
zeigen kénnen, daB die laminare 
Plattenstro6mung oberhalb einer 
kritischen Reynoldsschen Zahl in- 
stabil ist gegeniiber wellenférmigen 
Stérungen. Spater wurden diese 


Abb. 2. Taylor-Gértlersche Langswirbel in der Str6mung an einer 


konkaven Wand (nach H. G6rr.eEr [2]) Wellen durch G. B. SCHUBAUER 


und H. K. SkRAMSTAD [8] experi- 
mentell nachgewiesen und als Anfangsstufe des laminar-turbulenten Umschlags 
erkannt. Auf die Frage, wie der weitere Ubergang der wellenférmig gestérten, 
aber immer noch laminaren ebenen Strdmung in die turbulente Str6émung er- 
folgt, konnte aber bisher keine befriedigende Antwort gegeben werden. 

Die Untersuchungen der §§9 und 10 zeigen, daB eine gegen Tollmien-Schlich- 
ting-Wellen instabile Blasiussche Plattenstromung in den Talern dieser Wellen 
wiederum instabil sein wird gegeniiber St6rungen vom Typ der Taylor-Gortler- 
Wirbel, wenn nur die Wellenamplitude in der GréBenordnung von 104 Ver- 
drangungsdicken der Plattengrenzschicht liegt. Sobald also die Amplituden der 
aus irgendwelchen zufalligen Stérungen entstandenen Tollmien-Schlichting- 
Wellen die GréBe von 10-4 Verdrangungsdicken erreicht haben — und dieses 
wird in Anbetracht der auBerordentlich kleinen kritischen Amplitude sehr bald 
der Fall sein —, wird die Strémung gegentiber Langswirbeln sekundar instabil 
sein und somit in eine dreidimensionale Strémung tibergehen. Da erfahrungsgemaB 
auch Taylor-Wirbel sehr bald nach ihrer Entstehung instabil werden — vel. etwa 
die Experimente von F. Scuuttz-Grunow und H. HEIN [9] — liegt die Ver- 
mutung nahe, da die auf Grund der Sekundarinstabilitat entstandenen Wirbel 
sehr bald gegen weitere Stérungen instabil werden, so daB sich méglicherweise 
der laminar-turbulente Umschlag auf dem Wege iiber Instabilitaten héherer Ord- 
nung erklaren laBt. 

Aber auch fiir andere Strémungsprobleme ist diese Frage von Interesse. Wir 
weisen hier nur auf das folgende Beispiel hin: Auf Grund theoretischer Uber- 
legungen ist durch Absaugen cine Stabilisierung der Grenzschicht gegeniiber 
Tollmien-Schlichting-Wellen zu erwarten. Saugt man nun aber periodisch — etwa 
durch Schlitze quer zur Strémungsrichtung — ab, so bekommen die Stromlinien 
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schwache Welligkeiten. Damit liegt aber eine fiir die Wirbelinstabilitat giinstige 
Situation vor. Unsere Untersuchungen zeigen nun, daB die durch eine periodische 
Wandwelle gestérte Plattenstromung — deren Stromlinien sind denjenigen bei 
periodischer Absaugung sehr ahnlich — in den Wellentalern gegeniiber den wirbel- 
artigen Stérungen instabil ist, wenn nur die Amplitude in der GréSenordnung 
von 10-* Verdrangungsdicken der Grenzschicht liegt. Die Stabilisierung gegeniiber 
Tollmien-Schlichting-Wellen wird also erkauft durch eine Tendenz zur Instabilitat 
gegeniiber sekundaren Langswirbeln. 

Die mathematischen Untersuchungen der vorliegenden Arbeit lassen sich 
weitgehend im Rahmen einer linearisierten Stérungstheorie durchfiihren. Dabei 
werden nur spezielle St6rungen einfachster Bauart betrachtet. Erst ihre linearen 
Uberlagerungen im Sinne einer Fourier-Synthese geben die allgemeinen Stérungen 
wieder, die den Untersuchungen sowohl als primare wie auch als sekundare 
Stérungen eigentlich zugrunde liegen. Dabei entspricht die Linearisierung der 
primadren Stérung gerade der Taylor-Gortlerschen Linearisierung hinsichtlich 
groBer Kriimmungsradien. Zur Rechtfertigung derartiger Linearisierungen wird ein 
den Methoden der St6rungsrechnung nachgebildetes iteratives Verfahren benutzt. 

Der Ansatz wirbelf6rmiger Sekundarstérungen fihrt auf ein Eigenwert- 
problem mit der Wellenamplitude — genauer mit deren Faktor «¢ beziiglich der 
in einer bestimmten Normierung gedachten Amplitudenverteilung — als Eigen- 
wert. Diese sekundaren Langswirbel werden in den Talern der gegeneinander 
phasenverschobenen primaren Wellenstérungen angesetzt. Es zeigt sich, daB der 
EinfluB der Phasenverschiebung auBerordentlich klein ist, und daB sich dieser 
als St6érungskorrektur des entsprechenden nicht phasenverschobenen Problems 
sehr leicht berechnen 1aBt. 

Diese Arbeit ist durch eine Reihe anregender Diskussionen beeinfluBt worden. 
Besonderen Dank schulde ich meinem verehrten Lehrer, Herrn Prof. Dr. H. GOrT- 
LER, durch den ich auch zur Beschaftigung mit dieser durch ihn so geférderten 
Fragestellung angeregt wurde. Ebenso wurde mir bei der Durchfithrung der 
numerischen Rechnungen von mehreren Seiten Hilfe zuteil. Besonderen Dank 
schulde ich hier Frl. H. Kompe, die die Hauptlast der Rechnungen getragen hat. 


I. Neubehandlung der Taylorschen und Gortlerschen Instabilitaten 
§ 1. Das Eigenwertproblem der Taylor-Wirbel 

Den Ausgangspunkt fiir alle Uberlegungen bilden die in der Einleitung be- 
schriebenen Taylor-Wirbel. Ihre Theorie soll hier deshalb kurz skizziert und 
zugleich in einer fiir das folgende zweckmaBigen Weise erweitert werden. 
G. I. Taytor [1] selbst betrachtet namlich nur den Grenzfall, in dem die Zylinder- 
radien groB sind im Verhaltnis zu ihrer Differenz. Wir werden hier die Lésung 
in Form einer Reihenentwicklung angeben, die fiir alle geniigend groBen Radien 
brauchbar ist und bei der die Taylorsche Losung als Ausgangsnaherung auftritt. 
Die Koeffizienten lassen sich dabei nach den Methoden der sog. St6érungsrechnung 
bestimmen. 

In diesem Zusammenhang ist es zweckmaBig, auch den Taylorschen Grenz- 
fall mit anderen, fiir die numerische Durchfiithrung einfacheren Methoden zu 
gewinnen als dies in [1] geschehen ist. Wahrend namlich Taytor [J] die Differen- 
tialgleichungen durch einen Reihenansatz mit Bessel-Funktionen lést, werden 
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wir im Anschlu8 an die Vorgehensweise von H. GORTLER [2] und G. HAMMER- 
Ln [10] bei dem verwandten Eigenwertproblem der Gértler-Wirbel (vgl. § 3) 
das Eigenwertproblem in ein System von Integralgleichungen verwandeln und 
den kleinsten Eigenwert dann mit Hilfe dieser Integraloperatoren iterativ be- 
stimmen. Wir fiithren diesen Lésungsweg auch deshalb hier an diesem einfachsten 
Beispiel einer Wirbelinstabilitat durch, weil er sich in geeignet abgewandelter 
Form ebenso zur Lésung des komplizierteren Eigenwertproblems der wirbel- 
artigen Sekundarstérungen eignet. 

In den §§1 und 2 seien alle Lésungen mit der Differenz D der Zylinder- 
radien, alle Geschwindigkeiten mit der Drehgeschwindigkeit U des inneren Zylin- 
ders, der Druck mit @U? und die Zeit mit D/U dimensionslos gemacht; ferner 
sei eine Reynoldssche Zahl eingefiihrt durch Re = UD/y. Bezeichnen dann 7, 
und z wie tiblich Zylinderkoordinaten, u, —v (diese Wahl des Vorzeichens ist 
zweckmaBig im Hinblick auf die spateren Grenzschichtuntersuchungen) und w 
die Geschwindigkeiten in y,7 und z-Richtung und R den Radius des auBeren 
Zylinders, so folgt aus den Navier-Stokesschen Gleichungen und der Kontinuitats- 
gleichung fiir die Geschwindigkeits- und Druckverteilung dieser stationdren 
Grundstrémung! R_1 Re 


M SSbals Tips Sal iN ); 
(4.4) ares 
w=, 


r 


p =polr) = [ “ar. 
Zum Nachweis der Instabilitat dieser Stromung bei geniigend groBen Rey- 


noldsschen Zahlen setzt G.I. TayLor [1] wirbelf6rmige Stérungen der Grund- 
str6mung (1.1) an in der Form 


I 


u,(r) cosa ze”, 


> 
| 


v, (7) cos az e", 
w, (7) singz e”’, 


pb = py (r) cosoze”. 


@y 
| 


Dieser Ansatz ist mit den linearisierten Stérungsdifferentialgleichungen vertrag- 
lich, falls die Stéramplituden dem folgenden System von Differentialgleichungen 
geniigen® (Striche bedeuten Ableitungen nach dem Argument 7) 


1 2 
Uy — Oo Uy, 


Uy 
i 
y2 


y Rew, — Re(up +S) = uf +4 
(4.3) y Rev, + Re ae u, — Rep, = 1 +— 


Wh 1 2 
yA a Ay ape as ST 


” 1 
y Rew, —o Rep, =v, +— wv, — 0 w,, 
v 


(1.4) vy + — v,—ow,=0. 


1 Vgl. H. SCHLICHTING [11], S. 63. 
; 2 Die Gleichungen (1.3) folgen aus den Navier-Stokesschen Differentialgleichungen 
fiir Zylinderkoordinaten, die man bei H. SCHLICHTING [ZOOS Zeyh aa Up=U, U,= —v 


und v,= w findet; ebenso die Gleichung (1.4), die aus der dort inul 
pea eie Grates etal), er dort angegebenen Kontinui- 
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Hinzu kommen die Randbedingungen fiir das Verschwinden der Stérungen an 
den Zylinderwanden 
(1.5) =), =v, =O tir y—R und +7 —R —1. 


Aus (1.3) und (1.4) lassen sich #, und w, leicht eliminieren. Die beiden rest- 
lichen Differentialgleichungen lauten dann im Falle neutraler Stérungen (y =0) 
nebst den zugehérigen Randbedingungen 


a fhe al 1 2 — ’ Ug 
ede y dv # — 08) 4, = — Re(ug +“) u, 


ye) a? tl o@ 4 2 U, 
Lal peas ope Ye. 
(Gat+ sa a) a, 20° Re tty, 
(1.7) “0 — vO ir y= ind 7 =k — 1. 


ZweckmaBigerweise fiihrt man nun 7 =R —, als neue unabhangige Verander- 
liche ein. Schreibt man noch 2, statt Re v,, so folgt aus (1.6), (1.7) (Striche be- 
deuten jetzt Ableitungen nach dem Argument 7) 


2 fy? Uo 
es [x6 R(1—n/R) Be 
(1.8) ne . 
ES ES — eee 3 o =A Uy, 
(1.9) ty Iu eae 0 Tits ae Ol md... 7 4 


mit dem Differentialoperator 


(ie 1 d 1 1 


—— —_— — — — — -_—- — co 
dn? R(i—nfR) dn R® (1—n/R)? 


(4.10) 


und dem Geschwindigkeitsprofil der Grundstroémung 


i LAR 4 | U] | 1 
we Molt) = So aiR RS 1-fR) 
Der Ansatz wirbelartiger St6érungen (1.2) fiihrt also auf Differentialgleichungen 
und Randbedingungen, die linear und homogen sind und somit ein Eigenwert- 
problem darstellen. Als Eigenwert kénnen wir die Reynoldssche Zahl Re auf- 
fassen, und zwar in Abhangigkeit vom Parameter o. 


Praktisch interessiert nur der kleinste Eigenwert, und zwar besonders dessen 
Minimum hinsichtlich o; denn bei der zugehérigen Reynoldsschen Zahl sind erst- 
malig zeitlich neutrale Wirbelst6rungen (1.2) mdglich. Diese Wirbel haben eine 
ganz bestimmte Periode in z-Richtung, namlich 2 = 2z/o*, wenn o =o* die Stelle 
des Minimums fiir den kleinsten Eigenwert ist. A ist also der Abstand der Achsen 
zweier gleichsinnig drehender Wirbel. Einen entsprechenden Abstand haben die 
Staulinien, die man am duBeren Zylinder sichtbar machen kann und die die 
Breite der Streifen festlegen. Man kann (fiir jedes o>0) zeigen, daB gréBere 
Reynoldssche Zahlen die Eigenwerte angefachter Stérungen sind und daB kleinere 
Reynoldssche Zahlen den Eigenwerten abklingender Stérungen entsprechen, so 
daB also die Eigenwertkurve neutraler St6rungen den Instabilitatsbereich von 
dem Stabilitaétsbereich trennt. 


ype 
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In dem von TAYLoR behandelten Grenzfall R>>1, bei dem 1/R gegen 1 ver- 
nachlassigt werden kann und somit nach (1.11) (7) durch ™9(7) =1 und L 
durch d?/dy?—o? zu ersetzen ist, vereinfacht sich das System (1.8), (1.9) zu 


Mt 2 tA y, 
| Uy9 — 0 U9 = U0 %40, 
1.12 bey 
tt € 
019 — 207049 + 04%) = — a. O Ugo M10; 


(1.13) hie elo= ye 0 tr =O) mid 74. 


Die Daten der Grundstrémung gehen also nur in der Kombination 2 Re?/R ein. 
ZweckmaBigerweise faBt man deshalb diese GrdBe als Eigenwert (9=2 Re2/R 
fiir den Grenzfall R>>1 auf?. 

Der betragsmaBig kleinste Eigenwert laBt sich numerisch am einfachsten 
nach dem Iterationsverfahren* ermitteln. Hierzu ist es zweckmabBig, das Eigen- 
wertproblem (1.12), (4.13) zuvor in ein System von Integralgleichungen zu tiber- 
fiihren®. Bezeichnen namlich G(y,7) bzw. H(y,7) die Greenschen Funktionen 
zu den Differentialausdriicken der linken Seiten von (1.12) unter den Rand- 
bedingungen (1.13) und ist K(y,7) =o? H(y,7), so sind die Gleichungen (1.12) 
und (1.13) aquivalent dem System von Integralgleichungen 


Myo(n) = - fe (4,7) Ugo (7) V9(7) an, 
(1.14) 
V9(%) = Hof K(n,7) Uo0(M) %10(n) 47. 


Die in 4 und 7 symmetrischen Greenschen Funktionen lauten dabei 


Sinan Sina(y —1), "oN 


a4 eel Tem 
eer sing Cit on Sit o(y — 1), ho, 
(A, (0)) fa(o (Q—1)) —6 A (2) 4,6 (—1))] felon) + 
(1.16) H(y,4)= [ (0) f(o (4 —1)) +o Hy eens. (7—1))|A(on), n<%, 
(i (0)) fa(o 7) ie beara ai ee 
sea (0) hon) —o A, (0) f(07)] Alo(4—1)), n>1 
mit 
H,(o) = 1 Sin o 


20 o? Gin? c—o*’ 


H,(o) = sited See 


2 o2 Sin? o—o* 


* Fir das formal gleichlautende Eigenwertproblem der Gértler-Wirbel (§ 3) wurde 
dieser Parameter bereits von H. GOrTLER [2] als Eigenwert eingefiihrt. In seiner 
dimensionsbehafteten Schreibweise lautet der Parameter Uo=2 Re? DR. 

* Vel. hierzu L. Covratz [12], H. WieLanprt [13] oder H. BUCKNER [14]. Fiir das 
Eigenwertproblem der Gértler-Wirbel wurde es bereits von G. HAMMERLIN [10] 
verwendet. 

ne Fur das Eigenwertproblem der Gortler-Wirbel hat dies bereits H. GORTLER [2] 
getan. 
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und TAG) == 7 Oe 


fo(t) = Gin t — t oft. 


Die iterative Bestimmung des kleinsten Eigenwertes laBt sich nun wie folgt 
durchfiihren. Wir definieren eine Iterationsfolge durch 


1 
uro(m) = — J G (1,7) Moo() Yo (%) 47, 
rn 0 
(1.17) : 
obo (1) = J K(q,7)) Moo (h) 10 (7) 47. 


Bezeichnet dann L{v, 9] irgendein lineares Funktional von V9, ZB. den Funk- 
tionswert von v49(7) an einer Stelle 7 =7*, so konvergiert L[v%)]/L[v%t?] fur 
y—> co gegen diesen kleinsten Eigenwert, falls nur die Ausgangsnaherung v{} (7) 
an der zum kleinsten Eigenwert gehérenden Eigenfunktion v,9(y) beteiligt® ist. 
Diese Bedingung kann bei einer numerischen Durchfiihrung des Iterationspro- 
zesses als erfiillt angesehen werden. 

Obwohl sich wegen % (7) = die Iterationsfolge analytisch auswerten lieBe, 
falls nur v{}(y) hinreichend einfach gewahlt wiirde, haben wir des geringeren 
Aufwandes wegen, nicht zuletzt auch im Hinblick auf die spateren Verallgemeine- 
rungen, die Integrale numerisch berechnet, und zwar unter Benutzung der Trapez- 
regel’?. Dies wurde fiir o=0,; 2,0; 3,0; 3,4; 4,0 und 5,0 durchgefiihrt, und zwar 
mit der Schrittweite  =0,05. Die Ergebnisse sind in Abb. 3 in § 2 aufgetragen. 
Es gibt also einen kleinsten Eigenwert ,,- 3400 fiir op 3, was einem Abstand 
gleichsinnig drehender Wirbel von etwa zwei Spaltbreiten entspricht. Dieser Wert 
laBt sich experimentell bestatigen. 

Es sei noch bemerkt, daB sich die Existenz eines kleinsten positiven Eigen- 
wertes in voller Strenge beweisen la8t. Durch Einsetzen der ersten in die zweite 
Gleichung laBt sich das System (1.14) noch auf eine einzige lineare homogene 
Integralgleichung zweiter Art zuriickfiihren, deren Kern in 0<y, 7<1 positiv 
ist. Nach R. Jentzscu [15] gibt es daher® fiir jedes o>0 einen reellen positiven 
Eigenwert j)(c), der einfach ist und kleiner als die Betrage aller anderen Eigen- 
werte; fiir diesen wechseln die zugehérigen Eigenfunktionen %19(7) bzw. v9(n) 
das Vorzeichen nicht. 

Zum Schlu8 erwahnen wir noch, daB sich der kleinste Eigenwert fiir o—0, 
wie o” verhalt. Man folgert namlich aus (1.15), (1.16) 


leita? n<%, 
Gia n>; 
AV lal che es aN), <7) 
ign ten—snl. Yon: 


6 Zum Begriff der Beteiligung vgl. H. WIELANDT [13]. 

7 Die genauere Simpson-Regel ist wegen der Unstetigkeit in den Ableitungen der 
Greenschen Funktionen nicht durchgangig zu benutzen; aus Griinden rechnerischer 
Einfachheit haben wir deshalb hier und spater mit der Trapezregel gearbeitet. 

8 In gleicher Weise hat bereits G. HAMMERLIN [10] die Existenz des kleinsten 
Eigenwertes fiir das Eigenwertproblem der Gortler-Wirbel gezeigt. 


(1.48) lim G(n, 7) = Go(y, 4) = 


o—>0, 


(19) Jim Hf) =Holnf) =| 
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Damit reduziert sich das Eigenwertproblem (1.14) fiir o—0, auf ein (nach 
R. Jentzscu [15] lésbares) Eigenwertproblem® mit dem Eigenwert dim (9 (a) 0”). 


§ 2. Lésung des Eigenwertproblems durch Storungsrechnung 


Das Eigenwertproblem der Taylor-Wirbel wurde in §1 nur fiir den bereits 
anderweitig von TAYLOR [1] untersuchten Grenzfall R>1 gelést. Wir setzen 
nun zur Lésung von (1.8), (4.9) bei beliebigem R den Eigenwert uw =2 Re?/R 
und die Eigenfunktionen ™, und v, als Reihenentwicklungen nach Potenzen von 
4/R an: | ne i 

eee OM bas T Me a ee 


1 


gin 

(2.1) Hy hu ae r M9 fase 7 
4 1 

% = %07 "11 Cie ea a 


Ebenso denken wir uns den Differentialoperator (1.10) und das Geschwindigkeits- 
profil (1.11) nach Potenzen von 1/R entwickelt 


Sy Ee d Te, 
La a “ R sc h “ie 
(2.2) 
Uo (9) = Uoo(H ) + 02 (1) + M92") a5 
mit 
Uo() =, 
2 
Moi (%) = oe 
3 2 
sie EE ae 


Der Koeffizientenvergleich der nullten Potenz fiihrt dann gerade auf das System 
(1.12), (1.13), so daB wy im folgenden den durch (1.17) ermittelten kleinsten 
Eigenwert bedeutet. Der Koeffizientenvergleich bei der ersten bzw. allgemein 
k-ten Potenz liefert die Systeme 
Uy — 07 yy = M90 %11 + Mor %10— Moo%10 + “o> 

(2.3) 1} — 207 v4, +047, = — [ig F? Ugg U1 — fly O? Up 9% 9 — 

— flo Up My 9 — fg F? Myo M194 + 29 — 207 V49; 
BA eo ioe Wh eee ee 

Vip — 207 014 + OF Vy, = — [ig O? Ugg Uy _ — My O? Mpg Mg + °°*5 


(2.5) My ,=Uy,=Vj,=0 fir »=0 und y=1. 


® Die Normierung der Eigenfunktionen denken wir uns etwa durch Hee V109(y) =1 
<1 
fiir o>0 cuecbsernit: Dann existiert offenbar Jim P10 (yn) und nee (1.14) auch 
dim vso(n) sowie Jim fH (7,4) too (9) trol) a8. 
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In (2.4) sind dabei alle Entwicklungskoeffizienten von (2.1) angegeben, die den 
Index & tragen, also alle GréBen, die aus (2.4) ermittelt werden miissen. (2.4) 
stellt ein inhomogenes System von Differentialgleichungen fiir die Funktionen 
u,,,(y) und v,,(7) dar, das unter den homogenen Randbedingungen (2.5) zu lésen 
ist; der inhomogene Teil dieser Differentialgleichungen enthalt dabei die noch 
nicht naher festgelegte GréBe 1. Nun ist das zu (2.4) gehérige homogene System 
von Differentialgleichungen [es ist von k unabhangig und mit dem System nullter 
Ordnung (1.12) identisch] unter den homogenen Randbedingungen lésbar; dieses 
ist ja gerade durch die Bestimmung des kleinsten Eigenwertes erreicht worden. 
Also! ist das inhomogene Randwertproblem (2.4), (2.5) im allgemeinen nicht 
lésbar, sondern nur dann, wenn der inhomogene Teil? einer weiteren Bedingung 
gentigt. Diese Lésbarkeitsbedingung laéBt sich durch geeignete Wahl der noch 
freien GrdBe mw, stets erfiillen und umgekehrt ist , durch diese eindeutig fest- 
gelegt. Die Funktionen «,,(y) bzw. v,,(y) kann man dann anschlieBend (bis auf 
additive Vielfache von (7) bzw. v49(y)) bestimmen. Die GréBe yu, ergibt sich, 
ohne zugleich auch u4,,(4) bzw. v,;(7) ermitteln zu miissen; diese Entwicklungs- 
koeffizienten der Eigenfunktionen benétigt man erst zur Berechnung von /q,+4. 
Man bezeichnet diese Methode der Lésung eines Eigenwertproblems durch Zuriick- 
fiihrung auf ein benachbartes, leichter lésbares Eigenwertproblem als Stérungs- 
rechnung?. 

Zur Durchfiihrung der St6rungsrechnung miissen wir noch das zu (1.14) 
transponierte System lésen 


So (7) = %o0(7) JK, 7) %o(n) 4n, 
(2.6) : 


ro(n) = — Mo woo(n) J Gm 7) Sq (%) 47. 


Dieses hat nach bekannten Satzen aus der Theorie der Integralgleichungen* 
dieselben Eigenwerte wie das System (1.14). Die hier bendtigte Eigenfunktion 
79(y) des kleinsten Eigenwertes errechnet sich also wieder iterativ entsprechend 
(1.17). Insbesondere wechselt 7(7) das Vorzeichen nicht und es ist somit 


1 
(2:7) I sh 79 (1) %10(y) an = O. 


Zur einfacheren Durchfiihrung der St6rungsrechnung fiir k21 ersetzen wir 
auch die Gleichungen (2.3), (2.4) und (2.5) durch aquivalente Integralgleichungen. 
Dies ist wegen der gleichen Differentialausdriicke auf den linken Seiten von (2.3), 


1 Wir machen hier und im folgenden Gebrauch von der bekannten Alternativaus- 
sage tiber die Lésbarkeit inhomogener linearer Randwertprobleme; vgl. etwa 
E, KamxKeE [16], S. 184. 

2 Zu diesem gehért neben den in (2.4) nicht mehr angeschriebenen Gliedern auch 
der u,(7) und v,,(y) nicht enthaltende Term ug 67 Ugo Uy9- 

3 Vel. zu dieser etwa H. BUCKNER [1/4]. Die Storungsrechnung zur Berechnung 
des Parameters uw und damit der Reynoldsschen Zahl Re nutzt also die Vorteile der 
»geometrischen™ Linearisierung [gemeint ist hier die Vernachlassigung 1/R gegen 1 
bzw. spater (§§ 5ff.) die Entnahme der Primarwellen aus einer linearisierten Theorie] 
und sucht den dadurch bedingten Fehler nachtraglich durch Anbringung von Korrek- 
turen wieder auszugleichen. Ware andererseits Re bekannt, so kann man dieses 
Verfahren als eine iterative Methode zur Bestimmung der GréBe 1/R auffassen. 

4 Vgl. etwa W. SCHMEIDLER Lb So ZOse 
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(2.4) und den gleichen Randbedingungen (2.5) wie bei dem System nullter Ord- 
nung (4.12), (1.13) mit denselben Greenschen Funktionen (1.15), (4.16) méglich. 
So ergibt sich fiir k= 1 


1 
ty ,(n) = — J G(n,%) Ugo (7) V1 %(%) aici 
(2.8) or pic uae eat yaa 
Vi) lo LEON) Uoo (7) Uy p (1) aH+ te J K(9,%) 4o0(%) U9 (7) dn+--:. 
0 


Durch Einsetzen der ersten Gleichung in die zweite erhalt man eine inhomogene 
Integralgleichung fiir v, ; (7) 


A 

V4(9) = — Hod J K(n, i) Moo (7) G (H, ) M00 () Yah) 47 AH + 
(2.9) : i ; un th 
tia J K(n, 7) U0 (7) My 9 (7) dn = K(9,%) &%-1(4) 4H 


mit einem durch die vorhergehenden Naherungen bestimmten, hier aber nicht 
naiher angegebenem g,_,(7). Die bereits friiher erwahnte Losbarkeitsbedingung?® 
besteht nun darin, daB der inhomogene Teil dieser Integralgleichung fiir v,; (7) 
orthogonal ist zu der (bis auf einen willkiirlichen Faktor) eindeutig bestimmten 
Lésung 7(7) des transponierten Systems (2.6); daher muB gelten 


(2.10) Pad J 79 (9) K(n, 7) 4oo(%) %10(n) 477 Fine 79 (n) K(y, 7) &x-1(4) 44 an =0. 


Der Faktor von , ist dabei nach (1.14) und (2.7) Z/u9==0; also ist uw, durch diese 
Lésbarkeitsbedingung eindeutig bestimmt zu 


ale ial 
(2.41) n=" [ fro K(n, %) &e—1(H) 4H an. 
00 


AnschlieBend kann v,,(7) und damit ,,(7) nach einem Verfahren§® fiir derartige 
inhomogene Integralgleichungen oder naherungsweise durch Lésung des (2.9) bei 
Benutzen der Trapezregel approximierenden linearen Gleichungssystems berech- 
net werden. 

Die Gleichung (2.11) laBt sich auch so herleiten: FaBt man (2.9) fiir alle 
Werte 7 mit O<7<1 auf als Bestimmungsgleichungen fiir die Werte der Funk- 
tion v,,(7) und «,, so hat man aus diesen Gleichungen v,,(7) zu eliminieren, 
um jy, berechnen zu kénnen. Diese Elimination la8t sich durch Multiplikation 
von (2.9) mit 7)(7) und anschlieBender Integration erreichen, wobei man von der 


Symmetrie der Greenschen Funktionen (1.15), (1.16) und von (2.6) Gebrauch 
zu machen hat. 


Speziell 1aBt sich die erste Korrektur 4, des Eigenwertes durch reine Qua- 
draturen aus den Eigenfunktionen des Systems nullter Ordnung (1.14), ihren 
Ableitungen und den Eigenfunktionen des transponierten Systems (2.6) bestim- 
men, ohne zuvor die Korrekturen 1, (7) bzw. v,, (7) der Eigenfunktionen berechnen 


a Vgi. etwa W. SCHMEIDLER [17], S. 269. Wir werden diese Lésbarkeitsbedingung 
durch direktes Berechnen von y,; bestatigen. 


SV elretwa ja tORN (18).S. 72. 
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zu miissen. Die nach (2.11) zur Berechnung von m, benétigte Funktion g, (7) 
lautet nach (2.3) 


8o(%) =o Moo (%) cf G (7,7) bs 29 (7)) %o0(7) + 19 (7) + V9 (7) M91 (7) | dn — 


~ fo 41071) (tox (1) +11 %00(0)) +e U10 (7) — 2040 (n)- 


(2.12) 


Wir haben uns bei der numerischen Durchfiihrung auf die Berechnung der 
ersten Korrektur s4,(a) des Eigenwertes an den in §1 genannten o-Stellen be- 
schrankt. Die Ergebnisse sind in Abb. 3 mit eingezeichnet. Das Minimum dieser 
Kurve liegt dabei in unmittelbarer sy 
Nachbarschaft des Minimums von 
Mo(c); der Abstand der Wirbel, die 
sich einstellen, ist also fiir groBe R 
praktisch unabhangig von dem spe- 
ziellen Wert von R. 


4000 


Fiir o—0, verhalt sich auch gy 
(44 (co) wie o-?, wie man aus (2.11), 
(2.12) unter Beriicksichtigung von 
(1.18), (1.19) sowie des asymptoti- 2g 
schen Verhaltens von w,(o) sofort 
folgert. Die entsprechende Eigen- 
schaft laBt sich auch rekursiv fiir ;gg 
4, (co) beweisen, wie nach der Ab- 
hangigkeit des Eigenwertproblems 
(4.8), (1.9) von o nicht anders zu ; * Taam, 7 7, a 
erwarten war. oe 

Der Vorteil der St6rungsrechnung Abb. 3. Kurven p9(a) und 44, (c) fiir das Eigenwertproblem der 

ae Taylor-Wirbel 

zur Losung des von 1/R abhangenden 
Problems (1.8), (4.9) liegt neben der Einfachheit darin, daB wir so in dem Bereich 
ausreichender Approximationsgiite von (2.1) eine analytische Abhangigkeit von 
4/R erhalten, wahrend eine direkte numerische Lésung immer nur fiir den ein- 
zelnen Wert 1/R eine Lésung liefert. 


§ 3. Das Eigenwertproblem der Gortler-Wirbel 


Wie in der Einleitung bereits erwahnt, hat H. GORTLER [2] unter geeigneten 
Vernachlassigungen gezeigt, da eine Grenzschichtstr6mung an einer konkav 
gekriimmten Wand gegeniiber Langswirbeln instabil ist, falls nur Re, |/9/R 
gréBer ist als eine kritische Zahl. Da die Fragestellung der vorliegenden Arbeit 
als Verallgemeinerung bzw. Verbesserung dieser Untersuchung angesehen werden 
kann, wollen wir in den §§ 3 und 4 eine auf das Spatere zugeschnittene Darstellung 
der Theorie der Gértler-Wirbel geben. Wir denken uns dabei alle Langen mit 
der Grenzschichtdicke 6 (definiert durch den Wandabstand, wo sich die Ge- 
schwindigkeit um 1% von der Geschwindigkeit der AuBenstr6mung unter- 
scheidet), alle Geschwindigkeiten mit der AuBenstromung U, den Druck mit 
oU? und die Zeit mit 6/U dimensionslos gemacht; ferner sei wie in §1 eine 
Reynoldssche Zahl eingefiihrt durch Re = U 6/r. 
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Dem Instabilitatsnachweis von H. GOrRTLER [2] liegt das gegentiber den tat- 
sichlichen Verhiltnissen vereinfachte Modell zugrunde, nach dem es sich um die 
Grenzschichtstrémung an einer konkaven Wand konstanter Kriimmung 1/R 
handelt, deren ungestértes Geschwindigkeitsprofil in der Umgebung der betrach- 
teten Stelle von der Wandbogenlange unabhangig ist!. Bei Verwendung der 
gleichen Koordinaten wie in § 4 — die konkave Wand entspricht ja gerade dem 
AuBeren Zylinder, der Rand der Grenzschicht dem inneren — lautet damit die 


stationére Grundstromung 
(3.1) u=U(y), YV=0, w=0, p=Pfo(n) 
und der (1.2) entsprechende Ansatz wirbelartiger Storungen 


i =u,(n) cosa ze”, 


=> 


=, (nos 02 2", 
® =w,(n) sino ze”, 
pb =f, (y) cosa ze”. 


Diese Ansatze fiihren wie in §1 auf die Differentialgleichungen? 


! Uo 
53) Ns [x0 R(1—n/R) Ag 
Die, wees oa o — uy 


mit dem Differentialoperator 


d? 1 d 1 1 : 


B-4) 1 = Gq — Ra—al®) an FB a/R 


Da die Wirbel iitber den Rand der Grenzschicht hinausgreifen werden, hat man 
sinngema48 das Verschwinden der Stéramplituden erst fiir 7— co zu fordern; 
demgema8 lauten die Randbedingungen 


(3.5) Uy = 0, = 0, = 0 Tir a= oO ind? yee 
Der Ansatz wirbelartiger Stérungen fithrt also wieder auf ein Eigenwertproblem 
mit 4 =2 Re?/R als Eigenwert. 

H. GORTLER [2] behandelt wie G.I. TayLtor [J] nur den Grenzfall eines 
groBen Wandradius’; er vernachlassigt also, und zwar unter Berufung auf grenz- 


1 Dieses sind die einschneidendsten Voraussetzungen, da sie dem elliptischen Cha- 
rakter der Differentialgleichungen widersprechen, aber das Instabilitatsproblem auf 
ein rein lokal zu betrachtendes reduzieren. Die analoge Voraussetzung wird — wie 
bereits in der Einleitung erwahnt — auch in der Theorie der wellenférmigen Instabili- 
taten gemacht; vgl. W. ToLLMIEN [7] oder H. ScHLicHTING [19]. — Man erhalt auf 
Grund dieser Annahmen die gleichen Differentialgleichungen wie bei der Taylor- 
Stromung. Insbesondere kann man die Stéramplituden w,, v,, w, und py, gemaB (3.2) 
als Funktionen von 7 allein ansetzen und erhalt nur ein Eigenwertproblem eines Sy- 
stems gewohnlicher Differentialgleichungen; vgl. auch §§ 7 und 10. 

* Die Singularitat der Differentialgleichungen fiir 7 =F ist durch die beiden in 
der Anm. 1 erlauterten Annahmen bedingt. In dem von H. GOrTLER [2] allein be- 
trachteten Grenzfall R>1 tritt die Singularitat jedoch auBerlich nicht zutage; vegl. 
(3.6). Andererseits ist aber auch der unten folgende Ubergang von (3.3) nach (3.6) 


nicht durchweg gerechtfertigt, sondern nur fiir solche » mit 7<R, also zum Beispiel 
in der Grenzschicht; vgl. hierzu auch Anm. 3 und 4, 
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schichttheoretische Abschatzungen®, alle Glieder der Ordnung 1/R gegen solche 
der Ordnung 1. Diese Annahmen vereinfachen nach (1.12) die Differential- 
gleichungen (3.3) zu 
Vj9 — 207 Vy9 + G4 9 = — [gO Ugo M10, 
wobei %9(7) das Geschwindigkeitsprofil der Grenzschicht in nullter Naherung 
ist. Hierfiir soll im folgenden dasjenige der Blasiusschen Plattengrenzschicht 
genommen werden. 

Mit Hilfe der zu den linken Seiten von (3.6) bei (3.5) gehdrenden Greenschen 
Funktionen 


fiipuae f 
J [5° 'Sinon, eas 
(3.7) G(n,7) = 


“e-" Ginof, 1, 


fem chk ghee . 
a ee Neen) A 


= e 8) fa (7 —n)+ 4¥], 7” <7, 
gos [09 20? np +0 +H) +4} 
SBS Gi Ostet hash) Bi ot 

K(9,7) = 0° A(n, 7) 


iberfiithrt dann H. GOrTLER [2] das Eigenwertproblem in das System von 
Integralgleichungen 


(3.8) = A(y,n) = 


und der Abkiirzung 


U1 9(4) = eS G(,%) U0 (7) V1 0(%) dy, 
(3.9) : 
V49(4) = Ho J K (9,7) Mo0(%) 10 (4) 477. 


Gleichzeitig gibt er eine Naherungslosung an*. Eine erste strenge Lésung stammt 
von G. HAMMERLIN [10]. Sie beruht auf der Anwendung des Iterationsverfahrens 


3 Diese Vernachlassigungen, die gerade dem Ubergang zu dem nullten Glied einer 
Reihenentwicklung nach 1/R entsprechen, werden aber auch au®erhalb der Grenz- 
schicht vorgenommen. Sie fiihren darauf hinaus, da alle Stromlinien, also auch 
diejenigen auBerhalb der Grenzschicht, die gleiche Kriimmung erhalten. Da es nun 
die Stromlinienkriimmung ist, die die Wirbelentstehung bewirkt, begiinstigt das 
Gortler-Modell, genauer das Eigenwertproblem (3.5), (3.6) die Wirbelinstabilitat, und 
zwar um so mehr, je weiter die Wirbel tiber die Grenzschicht hinausgreifen. Dieses 
ist besonders fiir kleine Werte von o (groBe Abstande der Wirbelachsen) der Fall. 
So strebt dann auch der kleinste Eigenwert “)(c) von (3.6), (3.5) fiir o>0, monoton 
fallend einem endlichen Grenzwert zu, was der physikalischen Wirklichkeit nicht 
entspricht. 

4 Diese Naherung besteht im wesentlichen darin, daB das Integrationsintervall der 
zweiten Gleichung auf die Grenzschicht reduziert wird; bei der ersten Gleichung ist 
dieses nach (3.10) ohnehin der Fall. Damit wird die Stromlinienkriimmung auBerhalb 
der Grenzschicht gar nicht beriicksichtigt, d.h. die Situation fiir die Wirbelintensitat 
gegeniiber (3.6), (3.5) — vgl. Anm.3 — wieder verschlechtert. Insbesondere gehen 
die im Gértler-Modell, genauer im Eigenwertproblem (3.6), (3.5) enthaltenen, der 
physikalischen Wirklichkeit nicht entsprechenden Annahmen somit in die Gértler- 
sche Naherungslésung nicht ein. So verhalt sich auch der kleinste Eigenwert dieser 
Naherung wie o? fiir o>0,. 
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zur Bestimmung des kleinsten Eigenwertes und wird im folgenden (nur gering- 
fiigig abgewandelt®) wiedergegeben. Gleichzeitig weist G. HAMMERLIN [10] die 
Existenz eines kleinsten positiven Eigenwertes /“)(0) auf Grund eines Satzes von 
R. Jentzscu [15] nach. Wir haben diesen Beweis bereits in § 1 — auf den Fall 
der Taylor-Wirbel iibertragen — angegeben. 

Bei der numerischen Losung der Gleichungen (3.9) machen wir wesentlichen 
Gebrauch von der Eigenschaft eines Grenzschichtprofils, daB es stets eine Zahl M 
gibt, so daB mit ausreichender Genauigkeit gesetzt werden kann 


(3.10) Uyelh) =15 Yoo) SOnnte iy ews, 


Dann reduziert sich namlich das Integrationsintervall der ersten Integralgleichung 
auf O<7<M, und das Integral in der zweiten Gleichung kann fiir M7 0 
(bis auf einen nur numerisch zu ermittelnden Faktor) exakt berechnet werden. 
Definiert man analog (1.17) die Iterationsfolge 


M 
uyo(n) = re Gq, %) Yoo (7H) Yr0(%) a7, 
(3-11) oo 
vie (7) =JS K(q, 7) Moo(H) “x0 (iq) 47), 


0 
so ist nach (3.7) wegen 7 >7 


(3.12) Utd (see 9" irene 
mit 

M 
(3.13) E® = — [| © Gino F woo (i) vf (7) aH. 


0 
Die zweite Gleichung (3.11) vereinfacht sich dann fiir 7 <M (nur fiir diese Werte 


wird ja vYo (7) zur Berechnung von wf" (n) bendtigt) zu 


M 
vo” (7) ih K (1,7) Moo(%) e0(%) di) + 


a (ACRE BEV es EEG nel ween 
mi 
A” = Eo (2 1 Neate 
4 4o : 
(3.15) BY — —” (5% siiias ) e- 20M 
8o 160? ; 


Cc” — @” “ely ee 2oM 
8o 


Berechnet man die noch auftretenden Integrale numerisch (etwa nach der Trapez- 
regel mit der Schrittweite 70,1), so lassen sich die Iterationen durchfiihren. 
Der Kleinste Eigenwert ergibt sich wieder als Grenzwert des Quotienten 
v1 (ny) [vt (ym) fiir vy co. Wir haben fig (G) pfin 6 = 055 0)3 200,550,758 Oho 
1,5 berechnet und dabei eine Neuberechnung [20] des Blasius-Profils paces 


> Wir ftthren die _Iterationen nicht analytisch sondern numerisch aus. Dann 
braucht das Grenzschichtprofil %99(7) nicht durch eine analytische Funktion approxi- 
miert zu werden und das Verfahren kann auch fiir kleine Werte von o benutzt werden. 
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gelegt; fiir diese ist M=1,6. Die Ergebnisse sind in Abb. 4 in § 4 aufgetragen. 
Der kritische Wert des Gortler-Parameters ist also sehr viel kleiner als im Falle 
der Taylor-Wirbel. Dies liegt daran, daB der Strémungsbereich und damit der 
Bereich, dessen Stromlinienkriimmung die Wirbelentstehung bewirkt, sehr viel 
groBer ist als im Taylor-Experiment. 


§ 4. Gortler-Wirbel bei groBem Achsenabstand 
Die in Abb. 4 wiedergegebene Kurve des kleinsten Eigenwertes s,(c) des 
Systems (3.9), das wir entsprechend der ersten Gleichung (3.11) sowie (3.14) 
auch so schreiben kénnen 


M 
Uy (4) = st G(q,%) “oo(7) %10(H) 47), 
M 
Vy 0(7) =o J K(n, n) Mo (1) Uy (7) dint (Ana Byer (CyB) en, q=SM 


mit analog (3.15) definierten Konstanten A, B und C, fallt mit kleiner werden- 
dem o monoton und strebt insbesondere fiir o->0, einem endlichen Grenzwert 
zu. Diese zunachst tiberraschende Existenz von lim Mo(c) wurde erstmalig von 


(4.1) 


G. HAMMERLIN [10] auch in Strenge bewiesen, und zwar unter Benutzung eines 
Potenzreihenansatzes fiir die Eigenfunktionen in 0<7<1. Wir wollen diese 
Aussage hier ohne Zusatzforderungen ahnlich wie bei dem in §1, Anm. 10, 
gegebenen Beweis aus (4.1) folgern. 

Hierbei wird einerseits Gebrauch gemacht von den aus (3.7), (3.8) leicht her- 
zuleitenden Beziehungent 


2 : , Nese 

(4.2) lim G(y, 7) = Go(y7) =). i 
o>0, 1, NAN 

Se ea et 5 

4. lim H(n,4#) = Hy(n,7 =| cues i 
(4.3) lien Ea eee eto 50) anf tig, n>. 


Zum andern benutzen wir die fiir die Eigenfunktion u,9(7) des kleinsten Eigen- 
wertes giiltige Abschatzung 


(4.4) | “10(7)| S|E| e~7", 720 
mit 
M 
(4.5) L= — [ £Sin oF wolf) V4 0(7) a7). 
0 


Zu deren Beweis verifiziert man zunachst mit Hilfe von (4.1), (3.7) 
M 
(4.6) My 9() = Bee" — f © [Sin o (y A) wool) %10(%) 44. 
n 


In (4.6) hat der zweite Term, da v,9() nach dem Satz von R. JenTzscuH [16] 
das Vorzeichen nicht wechselt und y<7# ist, das entgegengesetzte Vorzeichen des 
ersten. Da aber auch ~,,(7) immer einerlei Vorzeichen hat, und zwar dasjenige 


1 Diese Grenzfunktionen erfiillen also nicht mehr die Randbedingungen fiir 7 > cv. 
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des ersten Summanden, kann der zweite Summand in (4.6) nur den Betrag des 
ersten verkleinern, d.h. es muB (4.4) gelten. Hierbei kann noch E eal gesetzt 
werden, da die Eigenfunktionen nur bis auf einen Faktor bestimmt sind. 

Die Existenz von jim fo(o) ergibt sich nun einfach wie folgt. Nach (4.4) ist 
u(y) fiir alle c>0 (fiir diese gibt es ja eine Lésung von (4.1)) stetig und be- 
schrankt. Nach (4.3) existiert dann auch 


M 
(4.7) J H(n, 7) Uoo(%) %10(%) 47 


fiir jedes o> 0 und ist ebenfalls beschrankt; der erste Summand in der Klammer 
der zweiten Gleichung von (4.1) strebt also, fiir c>0, gegen Null. Die beiden 
restlichen Summanden streben, wie man unter Benutzung von (3.15) folgert, 
gegen —(4)7?. Die zweite Integralgleichung reduziert sich also auf 

. 2 
(4.8) lim 2 9 (77) =— MoE ae 


o> 0, 


und damit die erste wegen (4.2) auf 


a] M 


Hol [aPainlid ag +2 [Pwo aah, 1 <M, 
(4.9) jim URAC == ey " 
Mog | iP M00(H) ai n>M. 
0 


Andererseits ist nach (4.6) wegen E =1 
| lim Uy 9(yn) = 1, 72M. 


Somit muB auch lim ju) (a) existieren, und es ist 


o—+0, 


(4.10) lim g(a) = 


o> 0, 


8 


i : 
J 7? Uo 9 (7) 27] 

Die in Abb. 4 wiedergegebene Kurve des kleinsten Eigenwertes juy(o) von 
(3.9) bzw. (4.1) besagt, daB Gértler-Wirbel nach dieser Theorie um so gefahr- 
licher sind, je gréBer der Abstand ihrer Wirbelachsen ist. Dieses widerspricht 
aber physikalischen Uberlegungen wie auch den Beobachtungen. Tatsdchlich 
gibt das in §3 benutzte Modell die Wirklichkeit fiir kleine Werte von o nicht 
richtig wieder. 

Wir haben schon in § 3, Anm. 3, darauf hingewiesen, daB das Modell (3.6), 
(3.5) die Verhaltnisse fiir die Wirbelentstehung dadurch begiinstigt, daB alle 
Stromlinien die gleiche Kriimmung 1/R erhalten, wahrend in Wirklichkeit die 
Kriimmungen nach auBen abklingen. Da sich die Maxima der Greenschen Funk- 
tionen G(7,7) und K(y,7) fiir o+0, mehr und mehr in den Bereich groBer 
Werte von 7 und 7 verschieben, d.h. in den Bereich, in dem eine zu groBe Strom- 
linienkrimmung angenommen wird, gibt das Modell die Neigung zur Wirbel- 
instabilitat fitr diese o-Werte falsch, und zwar zu giinstig wieder. Insbesondere 
hat der endliche Grenzwert Jim 4o(@) keine physikalische Realitat. 
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Rein auBerlich zeigt sich also der Fehler des Modells (3.6), (3.5) schon darin, 
da8B die Kurve pip (o) fiir o>0, monoton fallt und somit kein Minimum fiir einen 
positiven Wert o* hat. In einer spateren Arbeit hat G. HAMMERLIN [2/] gezeigt, 
daB die Eigenwertkurve neutraler Stérungen fiir o->0* itber alle Grenzen wachst 
und damit insbesondere ein Maximum fiir einen Wert o =o*>0 vorliegt, wenn 
man nur die bei H. GORTLER [2] vorgenommene Vernachlassigung des Terms 
Uyo/R gegen uo in der ersten Gleichung (3.3) riickgangig macht. Diese Vernach- 
lassigung ist namlich nur in der eigentlichen Grenzschicht gerechtfertigt (vgl. 
auch §3, Anm. 2), da fiir 7>M nach (3.10) %9(y) =0, aber u9(7)/R =1/R>0 
ist. Obwohl diese Anderung der Differen- 249 
tialgleichung nur klein ist, macht sie sich 
fiir o—>0, stark bemerkbar. 

Da sich die anderen Vernachlassigungen 57 
fiir nicht zu groBe Werte von 7 auch auBer- 
halb der Grenzschicht noch rechtfertigen 
lassen, lost G. HAMMERLIN [21] das um den jg 
Summanden mu /R gegeniiber (3.6), (3.5) 
modifizierte Eigenwertproblem 


ta — o% 4, = (woo — “22)0,, 
(4.11) R 


Mr 5 
vl — 20%, +044, = — [607 Ugo %, 


(4.12) Uy =, =; =0 fiir n=0 und noo, 0 


und zwar wieder durch einen Potenzreihen- Abb. 4. Eigenwertkurve nullter Ordnung j49(9) fiir 
ansatz der Eigenfunktionen. Einfacher und Pee, Giagaikiien wet Gn ae Ae 
zugleich mit einer analytischen Abhangig- ee LS a ea KS mee 
keit von 1/R erhalt man den _ kleinsten ery 
Eigenwert j«(c), indem man das neu hinzugekommene Glied als Stérung des 
Eigenwertproblems (3.6), (3.5) auffaBt und eine Stérungsrechnung mit dem Stor- 
parameter 1/R durchfiihrt. Wir wollen die Formeln hierfiir nicht im einzelnen 
angeben; die entsprechend (2.1) definierte Korrektur j,(c) ist in Abb. 4 mit auf- 
getragen. Die Durchrechnung zeigt, daB sich y,(o) fir o—+0, verhalt wie o%, 
also gdnzlich anders als in einem die Kriimmungsverhaltnisse richtig wieder- 
gebenden Modell, wie es etwa das in den §§7 und 8 entwickelte ist. Abb. 4 
enthalt weiter die Kurve (o) =) (o) + (4/R) 1, (o) fiir die Werte 1/R = 6/R =10% 
und 1074, die die GréBenordnung der von H. W. LIEPMANN [3], [4] bei seinen 
Experimenten benutzten Daten haben und fiir die die Glieder zweiter und héherer 
Ordnung vernachlassigt werden kénnen. Abb. 4 zeigt, daB die Lage des Minimums 
von dem Wert 1/R abhangt, und daB speziell fiir 1/k—>0 die Kurve (co) in die 
Kurve {49 (a) tibergeht. 


II. Wirbelinstabilitat einer durch eine periodische Wandwelle gestérten 
Schichtenstromung 
§ 5. Problemstellung und Lésungsmethoden 
Wir wenden uns nun der eigentlichen Fragestellung dieser Arbeit zu: der 
Untersuchung der Instabilitat wellenférmig gestérter Strémungen gegentiber 
sekundaren Langswirbeln. Bei der Herleitung der Differentialgleichungen fiir 
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derartige Wirbelstérungen wollen wir uns auf ebene Strémungen als ungestérte 
Grundstrémungen beschranken, deren Geschwindigkeitsverteilungen tiberdies als 
lokal von der Koordinate in Hauptstrémungsrichtung unabhangig angesehen 
werden kénnen. Wir erfassen damit einerseits die ebenen Schichtenstr6mungen 
wie die Poiseuille-Stromung, bei denen eine solche Abhangigkeit iiberhaupt nicht 
vorliegt, zum anderen aber auch naherungsweise eine jede Grenzschichtstromung, 
deren Geschwindigkeitsverteilung sich auf der Strecke einer Wellenlange der tiber- 
lagert zu denkenden Wellenstérung nur vernachlassigbar wenig andert. Diese 
Einschrankung der betrachteten Grundstrémungen ist einerseits daftir erforder- 
lich, daB die Stabilitatsuntersuchung auf ein Eigenwertproblem gewoéhnlicher 
Differentialgleichungen fiihrt, andererseits aber auch durch die Tatsache bedingt, 
daB nur fiir solche Str6mungen Amplitudenverteilungen von Wellenstérungen 
berechnet vorliegen. Trotz dieser das Problem wesentlich vereinfachenden An- 
nahmen diirften jedoch die erzielten Ergebnisse von allgemeinerer Giiltigkeit sein. 

Im folgenden seien alle Langen mit einer charakteristischen Lange, etwa der 
Grenzschichtdicke 6 (wie in §3 definiert durch denjenigen Wandabstand, bei 
dem die Geschwindigkeit den Wert der AuBenstr6mung bis auf 1% erreicht hat), 
alle Geschwindigkeiten mit einer charakteristischen Geschwindigkeit, etwa der 
Anstrémgeschwindigkeit U, die Zeit ¢ mit 6/U und der Druck # mit @ U? dimen- 
sionslos gemacht; Re bezeichne die mit den charakteristischen GréBen gebildete 
Reynoldssche Zahl U6/v. Dann lauten die die Str6mung beschreibenden Navier- 
Stokesschen Gleichungen und die Kontinuitatsgleichung fiir den Geschwindig- 
keitsvektor p und den Druck #: 


Ts 1 
(SrA) i grad p + Ab, 
(5.2) divt==.0: 


In den tiblichen kartesischen Koordinaten x, y, z mit Geschwindigkeitskompo- 
nenten uw, v, w haben dann die Geschwindigkeitsverteilung und die Stromfunktion 
der ungestérten Grundstrémung die Gestalt 

y 
(5.3) Uy =U(¥), Uy=0, W=0, Y=S u(y) dy. 


0 


In Anbetracht der zum Auftreten neutraler sekundarer Langswirbel erforder- 
lichen auBerordentlich kleinen Wellenamplitude kann die primire Wellenstorung 
einer linearisierten Theorie entnommen werden!. Wir brauchen daher nur 
Partialstérungen zu betrachten und kénnen unter Beschrankung auf zeitlich neu- 
trale St6rungen deren Stromfunktion ansetzen in der Form 


(5.4) epr(x, ¥, 1) = € R{p(y) e799} = eg(y) cos (ax — Bt +4(y)) 
mit 


i) q(v) = Vor (y) +93 (y), (vy) =aretg SE 


‘In Strenge hatte man die primadre Wellenstérung aus den nichtlinearisierten 
Storungsgleichungen durch einen Ansatz epr(x, y, t) + & Or(%, y,t) +--+ mit « als 
Entwicklungsparameter zu ermitteln und dann diese Stromfunktion den folgenden 
Uberlegungen zugrunde zu legen. i 
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a, B reell, wobei p(y) =¢,(y)+7¢;(y) Lésung der (linearen homogenen) Sommer- 
feld-Orrschen Differentialgleichung ist. Dabei ist p(y) in einer beliebigen aber 
festen Normierung zu denken. Der Tatsache, daB die Amplitudenverteilung ¢(y) 
einer neutralen Tollmien-Schlichting-Welle (8 > 0) als Lésung einer linearen homo- 
genen Gleichung bei homogenen Randbedingungen nur bis auf einen Faktor fest- 
gelegt ist, tragen wir durch die zunachst unbestimmte GréBe ¢ Rechnung. Bei 
der durch eine schwache Wandwelle erzwungenen Wellenstérung (6 =0) setzen 
wir die Wandamplitude mit der spater variabel zu denkenden GréBe ¢ als Faktor 
an; dann kann auch hier die Stromfunktion als Lésung der Sommerfeld-Orrschen 
Gleichung (mit 6 =0) in der Gestalt (5.4) angesetzt werden, da die einzige inhomo- 
gene Randbedingung gerade die Wandamplitude und damit die GréBe ¢ als Faktor 
enthalt. 

Die Vermutung einer Instabilitat gegeniiber sekundaren Langswirbeln griindet 
auf der Kriimmung der Stromlinien. Nach den Ausfiihrungen der §§1 bis 4 
wird man deshalb zu deren Nachweis den Geschwindigkeitsvektor und die Navier- 
Stokesschen Gleichungen in Richtung parallel und senkrecht zu den wellen- 
formigen Stromlinien zerlegen, um auch hier die Zentrifugalbeschleunigung 
explizit in den Gleichungen auftreten zu lassen. Zugleich wird man neue Koordi- 
naten &, 7 einfiihren, deren Linien 7 =const die wellenférmigen Stromlinien und 
deren Linien § =const die Linien gleicher Phase sind. Speziell sollen dabei 7 =0 
die Wand und € =0 ein Wellental beschreiben?. 


Fiihrt man tiber diesem System krummliniger Koordinaten entsprechend (3.2) 
einen Ansatz wirbelartiger Sekundarsté6rungen ein, so zeigen die Stérungsglei- 
chungen eine groBe Ahnlichkeit mit denen der §§1 bis 4. Insbesondere ergibt 
sich wieder ein Eigenwertproblem, dessen Eigenwert jetzt aber — da die Rey- 
noldssche Zahl durch die vorgegebene Wellenstérung als bekannt anzusehen ist — 
die Wellenamplitude der Stérung, genauer die der Kriimmung 1/R in den §§ 14 
bis 4 entsprechende GréBe ¢ ist. 

Im Falle 6 =0, also einer durch eine periodische Wandwelle erzwungenen 
Primarstérung, liegt sogar eine fast vollstandige Analogie vor. Insbesondere 
laBt sich, da ¢ fiir das Auftreten neutraler Sekundarwirbel nur sehr kleiner Werte 
bedarf, auch hier das Problem dadurch wesentlich vereinfachen, daB man zundchst 
den Grenzfall «<1 betrachtet und dann die Glieder hoherer als erster Ordnung 
in ¢ durch Anwendung der bereits in § 2, Anm. 3, angedeuteten iterativen Metho- 
den beriicksichtigt bzw. deren Einflu8 in einer diesem Vorgehen entsprechenden 
Weise abschatzt. Insbesondere ist der linearisierte Ansatz (5.4) nachtraglich 
gerechtfertigt. 

Ebenso 1a8t sich im Falle B>0, also bei mit der Phasengeschwindigkeit B/« 
stromabwarts wandernden Tollmien-Schlichting-Wellen als Primarstérungen die 
Analogie noch weiter verfolgen, jedoch liegt hier zugleich ein grundlegender 


2 Der Wahl des Koordinatensystems kommt eine groBe Bedeutung zu, weil durch 
sie zugleich die sekundaren Wirbelstoérungen festgelegt sind; wir werden namlich 
diese spater in der formal einfachsten Form (7.1) ansetzen. — Die Einfiihrung einer 
die Zustande gleicher Phase charakterisierende GroBe & erscheint zweckmabig, da wir 
uns, wie in der Stabilitatstheorie iiblich, auf eine lokale Betrachtung der Wellentaler 
zu beschranken haben und somit diese durch konstante Werte einer Koordinate 
kennzeichnen miussen. 

Arch. Rational Mech. Anal., Vol. 2 18 


262 HERMANN WITTING: 


Unterschied vor. Geht man namlich in diesem Falle zu einem mit der Phasen- 
geschwindigkeit B/x bewegten System iiber — dieses erscheint zweckmabig, da 
sich nach (5.4) alle Zustande gleicher Phase mit B/« stromabwarts bewegen und 
dann insbesondere die Stromlinien nicht mehr von der Zeit abhangen — so gibt 
es in der Nahe der kritischen Schicht y, — definiert durch u(y.) =B/« — ge- 
schlossene und insbesondere sich schneidende singulaére Stromlinien der wellen- 
formigen Grundstrémung; vgl. Abb. 6. Damit erhalt das Eigenwertproblem der 
wirbelartigen sekunddren Stérungen in den Gliedern héherer Ordnung in é fiir 
y =y, eine Singularitat, wahrend es in nullter Ordnung mit dem fiir 6 =0 erhal- 
tenen iibereinstimmt. In der Umgebung der kritischen Schicht diirfen also die 
Glieder héherer Ordnung nicht mehr vernachlassigt werden und auch das diese 
Glieder beriicksichtigende iterative Verfahren kann nicht mehr in der tiblichen 
Weise angewendet werden. Da iiberdies im Falle 6 =0 die Koordinate 7 zweck- 
maBigerweise anders eingefiihrt wird als fiir 8>0, und hier auch die Rechnungen 
nur bis auf Glieder hdherer als erster Ordnung in ¢ durchgefiihrt zu werden 
brauchen, betrachten wir in diesem Kapitel zunachst nur den Fall 8=0 und 
verschieben die Ausfiihrungen zum Fall 6>0 auf das folgende 3. Kapitel. 


§ 6. Einfiihrung neuer Koordinaten und Herleitung der Storungsgleichungen fiir den 
Fall B =0 
Wir behandeln zunachst den Fall einer durch eine periodische Wandwelle 
erzwungenen Wellenstérung (8 =0). Die Stromfunktion der wellenférmig ge- 
stérten Grundstrémung lautet nach (5.3), (5.4) bis auf Glieder hdherer Ordnung 
in € 


(6.1) p(x, ¥) = poly) + €¢(y) cos (ax +7 (y)). 


Ihre Stromlinien sind also durch konstante Werte von w(x, y) gekennzeichnet. 
Setzt man diese Konstanten an in der Form w (yn), so sind, wie man durch Taylor- 
Entwicklung verifiziert, die Stromlinien bis auf Glieder héherer Ordnung in ¢ 
durch konstante Werte von 


(6.2) n=y+ eQ(y)cos(ax+ty(y)), OG) = Bicake 


gekennzeichnet. Diese Darstellung ist auch an der Wand (vy =0) gestattet, da 
hier g(y) wie auch alle entsprechenden Anteile der Glieder héherer Ordnung in ¢ 
(vgl. §5, Anm.1) Nullstellen haben. Offenbar ist gerade Q(y) die Amplitude 


der Stromlinienwelle und 7 =0 die Gleichung der Wand, wahrend x = — + iy) 
ein Wellental beschreibt. Wir benutzen deshalb - 


g=x+—j(y), 
n=V +8 Q(y) cos (ax +7 (y), 


ee 


(Eis 


(6.3) 
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als neue Koordinaten. Dann ist 


Og oa2G ath re) 
Cre are 
fa) ap od) 0 
j = Sa 5 = (1 éS) ? 
(6.4) Oy ca 0€é on 
a eee 
iy meee 
CS 
Of = Ot. 
mit 
r =0Q(y) sin (ax +) (y)), 


s = Q’(y) cos (a x +7 (y)) — O(y)7"(y) sin (a * +7 (y)). 


Zur Zerlegung des Geschwindigkeitsvektors » in Richtung parallel und senk- 
recht zu den wellenférmigen Stromlinien definieren wir ein neues orthogonales 
Dreibein ¢,, es, ¢s gemaB 


g=(1+es)i+ter], 
(6.6) @=—erit+(1+és)j, 
Mees & 


wobei wie iiblich i,j, £ die Einheitsvektoren des x, y, z-Systems sind. Aus der 
Vektorzerlegung 


(6.7) p=aetbetce,=ui+tvj+wt 


folgt dann die Umrechnungsvorschrift fiir die Komponenten des Geschwindig- 
keitsvektors 
a= (es) er, 


(6.8) b=—eru+(1+6s)v, 


c=—wW. 


Fiir die Grundstr6mung der sekundaren Wirbelstérungen, also fiir die wellen- 
formig gestérte Strémung, verifiziert man 


(6.9) Ag =Uy(y), O=0, CoO. 


Zur Herleitung der Stérungsgleichungen sind die Navier-Stokesschen Glei- 
chungen (5.1), (5.2) beziiglich des Dreibeins (6.6) zu zerlegen und dabei die 
Koordinaten (6.3) als unabhangige Veradnderliche zu verwenden. Zunachst ist 


nach (6.7) 


de db de dc 
= | peed) | 2 a 
Carrs Way ee ner omer eae 


bei Benutzen der aus (6.6) durch Differentiation nach ¢ entstehenden Ableitungs- 
gleichungen 


dey ds dy 

eee a ee 
Cegek a aaY ds 
ae. Sat i ogee 


18* 
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folgt fiir dv/dt die Zerlegung 


dv da ds dy _ (ab dy , ds \ dc 
eel er Siar eb a (5 ee ayaa eb 7) eat Fhe 


Entsprechend ergibt sich fiir J») mit der Abkiirzung V=grad 


Av =(Aa+e2VaVs—e2VbVr+eadAs—ebArieyt+ 


(6.10) 
+ (4b6+¢e2VbVs+e2VaVr+ebAs+eadr)e,+ Ace. 


Zur Zerlegung von grad f gehen wir aus von der Darstellung 


Op . Op . 0 
grad p = oe oS £28 
ersetzen wir die Ableitungen gemaB (6.4) durch Ableitungen nach den neuen 
Variablen und die Einheitsvektoren t, j, £ gemaB (6.6) durch die Vektoren e,, €2, es, 
so erhalten wir 


gradp =(1—es er’) ee 


fos op , Op), , oP 
(—7'(t éS) er) og T ees I €3. 


In gleicher Weise ermitteln wir auch den Ausdruck div vp fiir die Kontinuitats- 
gleichung, indem wir von der Darstellung 
Ou ov ow 


div) = 4 = 
a Gas | ek Oz . 


ausgehen und die neuen Geschwindigkeitskomponenten gem&B (6.8), die Ablei- 
tungen nach den neuen Variablen gem&B (6.4) einfiihren. 


Damit lauten die Bewegungsgleichungen (5.1), (5.2) 


d ad da +\ @ 
eon na (1 ester j') e4 
+ (da +e2VaVs — e2VbVr+eaAs—edAn), 
db ds adr PAs, o ) 
(6.11) ait oa =—[(Era éS) ere + st|+ 
+45 (4b + 627bVs+e2VaVr+ebAs+eadn, 
ae Op 4 
ap sonar 
(1+estert;' ae | uel ES) EY Eo 
@ 
(6.12) at) og c les 


+ (1 +203) = be (5 a) 04 ery 
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mit den Differentialoperatoren 


Fim Get |(t+e’ cosas) +b{2j"(1-+0Q' cosaé) —e+ (a O+07'2) sinadl| A+ 
+0(1+26 Q’cosa—26Qj'sinaé) Fer, 
eVrV=e[(a? 0+ Qj") cosaé+ Q7' sina £] + 
+e[aQj' cosaé-+aQ'sinag] =, 
eVs V=e| © '(Q"—a? Q—Q7"?) cosaé—* (a? Q'+20'7'2+ 0777") sina] + 
+ 6[(Q” — OF?) cosa E—(20"7’ + 07”) )sina] =, 
(6.13) 


ec set 
+2[4 7 (4 +e Q’cosaé&—eQ7’ Spiga ke es 
lee 
; rr o€ 
+5 + £[(0"— 98 0 — 07°) cosaé — (20°F +07") sin af] 2 


a6 an + 
+ (1+ 2¢€Q’ cosaé—2¢Q7’ sina §) 35+ 


und den Abkiirzungen 


ae =eaa?Qcosaé+t eb(a Qj’ cosaé+aQ’sinaé), 
eS = — ca(aQj’ cosas +40! sina é) + 
+ eb[(Q” — Q7'?) cosa é — (20’7’ + Q7”) sin« é], 
(6.14) &Ar =e0(20'7' + 07”) cosaé + ea(O'— 0? 0 — 07?) sina é, 


eAs = 2(Q'” — a2 Q’—30'j'2—307'j") cosaé + 
+ 6(Q7/3+ 0207’ —3Q"79 “Tt — 307" — Q) "’’) sina €, 
(2° — 2) — 0(Q” — a2 Q — Q7'2) cosaé — 0(20"j’ + Qj”) sina. 


Dabei sind Q und 7 nach (5.5) und (6.2) zunachst Funktionen von y, also gemaB 
(6.3) als Funktionen von 7 auszudriicken. 

Um die Instabilitat der wellenfoérmig gestérten Grundstrémung (6.9), die nach 
Voraussetzung Lésung der Bewegungsgleichungen (6.11), (6.12) ist, gegentiber 
Sekundarwirbeln nachweisen zu k6énnen, leiten wir zunachst die Differential- 
gleichungen fiir beliebige Sekundarstérungen her. Wir setzen also zur Lésung 
von (6.11), (6.12) 

a = a(y) + 4(E,7,6, 7), 


b b(é,n,6, 7), 
C= DUSSET AEB i 
b =ho(n) +£(6.0.6,7). 


| 


(6.15) 
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Unter der Annahme, da8 eine Linearisierung hinsichtlich a, b, é, gestattet ist, 
eae sich dann die Seca zu 
+ a(1+eQ’ cosa £) 2 e 74 ay(1 +22eQ’ cosaé — 2€Q7 sin a €) b 
+ € a)[(Q” — sane a2 0) cosa & — (2Q'7'+ Q7”) sin « €] b+ 
— €2a)(a Qj’ cosx& +a Q’ sina é) a 


i F 
—(1— €Q’ cosaé +6207’ sinaé) 22 + 
(6.16) + (da + 2 2VsVa — 02Vr1Vb+ eAsa —e Arb), 
e€ 
2  ag(t +eQ' cosag) 22 + 2 2aya? Qcosaea 
1 + gi op ap 
= -— st (1—€Q’cosaé+ €Q7' sing &) —eaQsinag| 3 an =f 
+o (Ab + 22VsVb+ e2Vr Vai + eAsb+eAra), 
C ; a¢ ep ens 
 +a(ited cos « €) S = oS ae Ac, 


(I+6Q'cosa8) SE + +\—4 ‘“dteQ’ pee sina £)—ea Qsinag| 22+ 


(6.17) 4 (14260! cosaé—2eQj' sinaé) 2 
+ e[(Q” — Q7’?— a2 Q) cosa é — (2Q"j’ + Qj”) sinx €] 6 + 9¢ — 0- 

Dazu kommen die Randbedingungen 

(6.18) a@=b=é 


an den berandenden Wanden bzw. fiir 7—> oe 


‘S) 


§ 7. Wirbelansatz und Formulierung des Eigenwertproblems fiir den Fall B =0 


Wir betrachten nun speziell Stérungen vom Typ der Taylor- bzw. Gértler- 
Wirbel. Die Gleichungen (6.16), (6.17) zeigen, daB der (1.2) bzw. (3.2) ent- 
sprechende Ansatz 


I 


a, (yn) cosa e”*, 


ao & 


( 
(7.1) aie cos a¢ e”*, 


sing e”*, 

p = b(n) cosa e”* 
zwar die Abhangigkeit von den Variablen ¢ und 7 befriedigt, nicht dagegen die- 
jenige von &. Um auch diese erfiillen zu kénnen, miiBten die Stéramplituden 


a, 6, ¢, und #, als Funktionen von 7 und & angesetzt und dementsprechend aus 
partiellen Differentialgleichungen bestimmt werden}. 


? Wir haben eine solche Vorgehensweise an anderer Stelle (bisher unver6ffentlicht) 
bis zu einem gewissen Grade curchgenairy indem wir die partiellen Differential- 
gleichungen an drei benachbarten Stellen €= —h, £=0 und €=+h durch gewohn- 
liche Differentialgleichungen und dabei die Differentialquotienten in &- -Richtung durch 
entsprechende Differenzenquotienten ersetzt haben; die hierzu durchgefiihrten Rech- 
nungen bestatigten die GroBenordnung der im folgenden mitgeteilten Resultate. 
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Lai einer exsten Lisung, die wir hier vornchmen wollen, kinnen wir jedoch 
dex Tatsache Gebrauch machen, da8 sich Hie Kocffizienten der Stérungs- 
ichunge ee, en ee catwcies snd namlich de 


iplizie ORE lagi angle eal elon odex die £-Abbaneie- 
wird durch den Faktor BLE vermittct, der sich im Wellental — hier 


achungen 73 Vinen — nur in zweiter Ordnung mit £ andert. Somit kinnen 
fir diese exste Lisung dic £-Abbangigkeit vernachlassigen. 
Des Approximationsgrad diesa Naberungsioung Ait sich am besten erkennen, 
con man Ge SArisnung in dnem Wellental als Modell einer Grenzschichtsts6- 
mung Yangs dnex konkay gckripnmten Wand auffait und sie vergicicht mit den 
uspringsichen Gistles-Modell ($43 und 4) der konkaven Wand konstanter 
Krismmung (und konstant gekrisnmter Stromlinien). Dann wird namlich jetzt 
das Abkiingen da Stromlinienkriunmung fir y—> 20 richtig wiedergegeben®. 
Damit chalten wit beste Resultate als beim Gortles-Modell*, zB. cin der 
physikalichen Wirklichkeit enteprechendes Vethalten dex Kigenwertkurve far 
0-0, hea Ge Abhanggkcit des ksitixhen o-Wettes von der Wellenlinge der 
Vernachizwigt man aso die Abhangigkeit von £ im Wellental, dh. setzt 
gman in den Koclfizienten von (6.16), (617)  =0, 9 ist der speziclle Stérungs- 
ansatz (71) mit den Glechungen (6.16), (617) vertraglich, falls die Stérampli- 
tuden dem folgenden System von Difierentialgicichungen gentigen (Striche be- 
deuten Ableitungen nach dem Argument 7): 


ty + OylS + 24D) by + 24 — OF? — 2° Q) by — 22g OF a, 
ie (N+ 26D) — a, + (0 — OF*-— 2 OG + 
$US — O79) hy — 222 OF YEH" — 29-3972 

BOF 1m — 22 VF +2071, 

yh, 46 2gt Om, = — $y + 4 [A+2eVK)  K— Fb, + 

| +MY —OF*- BYU, + EUS — OF VU, + 222 OF a + 

| 4 MY —2Y 30 F237 (y+ 22 OF +207) a), 

$4=Oh + fob +2) 4-084 + (0 — O7?- 2 Qa], 


(7.2) 


(73) A+22Q)K +e" — O7?-# Ob +04=0. 


* Es «3 hier noch wont, 425 auch bei dex Behandlung dex Wisbelinstabilitat 
nex chenen Staupunkists@mung durch H. Goertiex 61 die Krismmung richtig 
berticksichtigt wird. Da in Giesem Problem jedoch keine Vagleichslange vorhanden 
ist, liclerte Giese Untereuchung keinen kritischen Wert ines Parameters. 

2 Auch won dea Vernachisesigung dea Anderung dex Stromlinienkriimmung in 

£-Richtung wird jetzt asst an sche viel spatter Sicile, namiich erst unmittelbar vor 
dex numerischen Auswertung, Gebrauch gemacht. 
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Diese sind unter den aus (6.18) mit (7.1) folgenden Randbedingungen zu lésen 
(7.4) a4 =b,=¢,=0 fir 7 =0 und y=. 


Der Stérungsansatz (7.1) fiihrt also auf ein Eigenwertproblem mit der die Wellen- 
amplitude festlegenden GréBe ¢ als Eigenwert; die Reynoldssche Zahl ist ja 
durch Vorgabe einer speziellen Wellenstérung bereits festgelegt. Damit ist eine 
zunachst rein auBerliche Analogie zu den in den §§ 1 bis 4 behandelten Problemen 
der Taylor- und der Gértler-Wirbel hergestellt. 

Diese Analogie laBt sich aber, wie nach der gleichartigen Behandlung der 
physikalisch verwandten Fragestellungen zu erwarten war, noch wesentlich weiter 
verfolgen und somit zur Lésung von (7.2), (7.3), (7.4) ausnutzen. Beachten wir 
zunachst, daB nicht nur die Wellenamplitude, sondern auch die Kriimmung der 
Stromlinien durch ¢ mitfestgelegt ist; sie betragt z.B. im Wellental ea%?Q. Spe- 
ziell bezeichne ¢ «?Q,, den Maximalwert der Kriimmung, der an der Wand (7 =0) 
angenommen wird und der somit der Wandkriimmung 1/R der §§1 bis 4 ent- 
spricht. Da nach den Theorien der Taylor- und der Gortler-Wirbel schon relativ 
kleine Kriimmungen zur Wirbelinstabilitat ausreichen, wird man auch hier kleine 
Kriimmungen und damit kleine Werte von ¢ als Eigenwerte vermuten, was durch 
die Resultate der numerischen Rechnungen auch tatsachlich bestatigt wird. 


Wir werden deshalb dem Vorgehen in den §§1 bis 4 entsprechend zunachst 
die GréBen f, und c, aus (7.2), (7.3) eliminieren, Re 6, durch 0, ersetzen und die 
dabei auftretende und dem Gortler-Parameter 2Re?/R entsprechende GréBe 
fs =2¢e Re®a?Qy, als Eigenwert einfiihren*. Dabei denken wir uns Q,, auf 1 nor- 
miert. Dann werden wir zuerst den Grenzfall «<1 betrachten, d.h. alle Glieder 
hdherer Ordnung in ¢ gegen 1 vernachlassigen; w selber werden wir dabei als 
eine GroBe nullter Ordnung in ¢ ansehen. AnschlieBend werden wir auf Grund 
einer Reihenentwicklung nach Potenzen von « eine die Glieder hoéherer Ordnung 
in ¢ berticksichtigende Korrektur berechnen und so eine verbesserte Lésung fiir 
das Gesamtsystem (7.2), (7.3) und (7.4) gewinnen. 

__ Es soll jedoch betont werden, da8 dennoch ein Unterschied gegeniiber den 
Problemen der §§ 1 bis 4 besteht, und zwar deshalb, weil jetzt die Reynoldssche 
Zahl Re bekannt und damit nur noch eine GréBe, namlich die Wellenamplitude e, 
in unseren Gleichungen frei ist. W&ahrend friither die Reynoldssche Zahl Re 
Eigenwert (bis auf den konstanten Faktor \/2/R) und die GréBe 1/R Stérparameter 
war, ist jetzt die Wellenamplitude ¢ zugleich Eigenwert (bis auf den konstanten 
Faktor 2 Re?a?Qy) und Stérparameter. Bei unserem Verfahren zur vereinfachten 
Berechnung des Eigenwerts handelt es sich somit eigentlich nicht um eine Sté- 
rungsrechnung, sondern um ein der Stérungsrechnung nachgebildetes Iterations- 
verfahren, bei dem die Glieder héherer Ordnung in ¢ nach und nach beriicksichtigt 
werden. In gleicher Weise ware auch das in § 2 behandelte Verfahren zu inter- 


pretieren, falls man dort die Reynoldssche Zahl als bekannt ansieht; vel. 
Anm. 3, § 2. 


‘ Geht man nicht von der Analogie zu den §§1 bis 4, sondern von den in der 
zweiten Gleichung (7.5) auftretenden Koeffizientenfunktionen aus, so hegt es wegen 


Qa=4 nahe, die GréBe w= 2e Re®a2qy als Eigenwert einzufiihren. Wir werden 
dies im Fall B>0 tun; vgl. §§ 9 bis 10. 
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AuBerdem scheint dem Verfahren dadurch eine gewisse Willkiir anzuhaften, 
da8 ein Ausdruck mit dem Faktor ¢ als Glied nullter Ordnung beriicksichtigt 
wird, wenn er den Faktor Re? hat. Dieser Tatsache entsprechend hatte man einen 
Ausdruck mit dem Faktor ¢?Re? als Term zweiter Ordnung zu interpretieren, 
wenn man ¢?Re? in der Form (eRe)? schreibt, aber als Term erster Ordnung, 
wenn man ¢*Re? in der Form ¢(¢Re?) schreibt. Diese Unbestimmtheit liegt 
offenbar daran, daB die Reynoldssche Zahl Re als Konstante eine geeignete 
Potenz der variablen GréBe ¢ bei dem Koeffizientenvergleich zu kompensieren 
vermag. Wir beseitigen sie durch die Festsetzung, die GréBen ¢ und Re? nur 
bei dem einen (durch seine numerische GréBe entsprechend seiner physikalischen 
Bedeutung *) ausgezeichneten Term nullter Ordnung zusammenfassen zu wollen ®. 
Damit ist die iterative Beriicksichtigung aller Glieder eindeutig vorgezeichnet. 

Die durch Elimination von f, und ¢, aus (7.2), (7.3) entstehenden Stérungs- 
gleichungen fiir neutrale Stérungen nach Ersetzen von Re 0, durch 0, lauten: 


05 Fie a ant Ele hh 
bi’ — 20%by + 04d, = — wo? Qaya, + €07Ly[a, dy). 


Dabei bezeichnen L,[a,, >,] und L,[a,, b,| die folgenden linearen Differential- 
ausdriicke in a, und 0,: 


Ly [as, 04] = 4520" dy + a9(Q” — Qi? — a2 Q) by — 2aya Q7' Re ay — 
20! a — (Q” — Qj? a2 Q) ai — 2(0"— Q 7") ak + 
+2007 1b; — [(Q” — a Q — Q7")'—207j'2— 
~ 07'iJa, +2007 $2.07") ody, 


(7.6) 


Lala, by] = — 3, O'BY + 40" by — [8 Q" + 2(0" — Q7'?— a? Q)] br” — 
— 4 [60 +4(Q" — O72 a2 Q)'] by — 
— 4, [20% +3(Q" — O72 2 Q)"] 84 + 


(7.7) + [40”— 2072+ 2(Q"— 072 a2 Q)] of — 
— 4 [(@”— 972-2 Q)"— 


pe een 1 eae Da) ree Te Osada) Ou 
+ 2a Q7’ Rea, + Re(2« Q’7’ + Q7”) a. 


5 Dieser fiir die Wirbelentstehung entscheidende Term geht aus dem Beschleuni- 
gungsanteil der 2. Gleichung (6.11) hervor, und zwar gerade aus dem Summanden 


& nae der durch die Drehung des Vektors'e, bedingt ist, also immer dann auftritt, 
wenn die Stromlinie gekriimmt ist. Die Tatsache, daB die Stromlinien konkav sind 


relativ zur Richtung wachsender Werte der Geschwindigkeit, wird erst bei der Lésung 
des Eigenwertproblems ausgenutzt. 

6 Insbesondere soll auch das aus dem Glied nullter Ordnung gema8 Q = Q(y) = 
Q(n) — €Q(n) Q’(y) hervorgehende Glied mit dem Faktor «Re? als Glied zweiter Ord- 


nung beriicksichtigt werden. 
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Die Randbedingungen gehen bei dieser Elimination tiber in 
(7.8) a,=b,=6,=0 fir 4=0 und 97>. 


Um nun das Eigenwertproblem fiir den Grenzfall e<1 und die Gleichungen 
fiir die Stérungskorrektur zu erhalten, setzen wir zur Lésung von (7.5), (7.8) 


[= Mo + e+ Mee + 


(7.9) hho tare +he,e+:*: 
by = by + Oyet+h.e+---. 


Dabei soll jz) den unter der Annahme e<1 hergeleiteten Eigenwert bezeichnen. 
Durch Koeffizientenvergleich erhalten wir dann als Systeme nullter und erster 
Ordnung (hdhere Potenzen von ¢ wurden ja bereits vernachlassigt, kénnen also 
nicht mehr beriicksichtigt werden) : 


aa 2 pli , 
Ay — FA) = Ayo, 
(7.10) 2 <p : . 
bio — 207 dy 9 + 04 by9 = — My OA Mo, 


ay 4 — Ca,= a bi t+ Ly[%o, O10], 
(7.11) : 
biy — 207 by, + 64d, = — Myo? Q ay yy — 07 Day Ay 9 + 07? Le [A9, Pr 9]- 


Die Randbedingungen (7.7) gehen durch (7.9) iiber in 
CA) Og Oe 0 0. AUT) a es or 


Damit tritt die Analogie zum friher betrachteten Gortler-Modell vollends 
zutage: Einerseits stimmen die Differentialausdriicke der linken Seiten von (7.10) 
sowie die Randbedingungen mit denen von (3.6) iiberein, so daB dieselben Green- 
schen Funktionen zur Uberfiihrung des Differentialgleichungssystems in Integral- 
gleichungen benutzt werden kénnen; andererseits stimmen auch die rechten 
Seiten von (7.10) bei entsprechender Abanderung der Stromlinienkriimmung mit 
denen von (3.6) tiberein. Bemerkenswert ist weiter die Tatsache, daB die Phasen- 
verschiebung der Wellen erst in (7.11) eingeht und somit nur in der Stérungs- 
korrektur beriicksichtigt zu werden braucht. 


Mit den Greenschen Funktionen (3.7), (3.8) gehen die Gleichungen (7.10), 
(7.12) tiber in 


a0) =~ f Gln.) aA) bro lf) a4, 
(7.13) ee 
bs 9() =H Mo J K(n,%) Q (7) Ao (77) Ay 9 (77) dn 


und lassen sich nun, wie in §8 noch naher ausgefiihrt werden soll, leicht iterativ 
lésen. Die Gleichungen (7.11), (7.12) werden mit Hilfe derselben Greenschen 
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Funktion iibergefiihrt in die Gleichungen? 


co 


aaa (a) = — f Gln f) a5 (A) (4) a9 — i G(n,A) Ls[a0(A), bro(A)] 49, 
(7.14) ba3(q) =o Kona) Q (i) do (fi) a4 (4) df + 
+ af K(f) Q (8) 40) a0 (4) \dij— IK( (7,4) La[ayo (4), b19(4)] 49. 


Aus diesen Gleichungen lassen sich wie friiher in § 2 die Korrekturen a,,(7) und 
5,,(7) der Eigenfunktionen eliminieren. Wir benétigen hierzu noch die Lésung 
7%)(yj) des zu (7.13) gehdrenden transponierten Systems 

So(7) = Q(%) 4o(n) J K(q,%) 70%) 47, 
(7.15) 


¥o(q) = — Ho 40(n) J G(, 7) S9(4) dH. 


ey Sse: 


Setzt man namlich a,,(7) nach der ersten Gleichung (7.14) in die zweite Glei- 
chung ein und integriert dann diese — noch mit 7) (7) multiplizierte — Gleichung, 
so ergibt sich als Korrektur des Eigenwertes® 


y= [J frome (1.9) QA) 40(A) CGA) La [aro, bso] dn 49 ai + 


00 0 


=i He f fro K(y,%) Lo[419, 10) 44 47. 
0 0 


(7.16) 


Dabei bezeichnet wieder 


(7.17) di =f rol by9(7) 4n. 


§ 8. Lésung des Eigenwertproblems fiir den Fall einer durch eine periodische Wand- 
welle gestérten Blasiusschen Plattenstrémung 


Zur Lésung des in den §§ 6 bis 7 hergeleiteten Eigenwertproblems der sekun- 
daren Wirbelstérungen benotigt man die durch die periodische Wandwelle 


(8.1) y =eR{A e**} = — ecos(~x +4) 
erzeugte Wellenstérung (5.4), d.h. die Amplitudenverteilung Q(y) und die die 


Phasenverschiebung beschreibende Funktion j(y). Diese Funktionen werden von 
uns im Falle der Blasiusschen Plattenstr6mung fiir verschiedene Werte von « 


7 Diese Gleichungen erhalt man auch, wenn man das System (7.5), (7.7) mit den 
Greenschen Funktionen (3.7), (3.8) in Integralgleichungen iberfiihrt. 

8 Eine entsprechende Formel erhalt man fiir die Korrektur auf Grund der Beriick- 
sichtigung der héheren Potenzen von e. 
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und Re berechnet!. Hier sollen nun die kritischen Amplituden der Wandwelle 
(8.1) berechnet werden, bei denen Instabilitat gegeniiber wirbelartigen Storungen 


auftritt. . 
Da Q(y)>0 ist fir 0<7< 0, gibt es nach R. JENTZSCH [15] einen betrags- 


maBig kleinsten Eigenwert fug(a) von (7.13), der positiv und einfach ist, und fiir 
den die zugehérigen Eigenfunktionen ihr Vorzeichen nicht wechseln. Wir kénnen 


fo(o) also nach dem Iterationsverfahren berechnen. 
Bei der numerischen Durchfiihrung machen wir wesentlichen Gebrauch von 


den Eigenschaften eines Grenzschichtprofils 

(8.2) ao (ny =4)) “QGy=0 tir neu 

und der aus der Sommerfeld-Orrschen Differentialgleichung herzuleitenden Be- 
ziehung” 


(8.3) Qn) =Ke-*" fir n=M. 
Dann vereinfachen sich namlich die Gleichungen (7.13) zu 
M 
(8.4) 44 9() = — f G(%,%) 40(%) B10(%) 27, 
0 
(8.5) b(n) =o S K(m,7) (7) 40(%) 410(7) 47, 
0 


wobei noch das Integral der Gleichung (8.5) fiir M<7< oo exakt berechnet 
werden kann. Nach (8.4) und (3.7) ist namlich 


(8.6) ao) =—Be fury 
und damit nach (8.3) und (8.5) 


M 
bso(n) = Ho { K(f) QC) oli) A404) 4) + 


(8.7) ; ; 
+(An+ Be "+ (Cn — Berk fir 17<M 


mit geeigneten aus (8.10) bis (8.13) unmittelbar abzulesenden Konstanten E, A, 
B,C. Fir 72M wird b,9(y) nach (8.4) nicht bendtigt. 


1 Im Gegensatz zu H. GORTLER [22], der die Wellenst6érung aus den Grenzschicht- 
gleichungen berechnet, haben wir dabei die Navier-Stokesschen Gleichungen zugrunde 
gelegt. Wir hatten also die Sommerfeld-Orrsche Differentialgleichung zu lésen, und 
zwar unter einer inhomogenen Randbedingung, die durch die Wandwelle (8.1) bedingt 
ist. Nach Entwicklung eines geeigneten Rechenverfahrens wurden diese Amplituden- 
verteilungen auf den elektronischen Rechenmaschinen ERMETH an der ETH Ziirich 
bzw. IBM 650 an der TH Darmstadt berechnet. Den Herren Professoren Dr. E. STIE- 
FEL, Dr. H. RutisHauser und Dr. A. WALTHER, die mir ihre Maschinen hierzu freund- 
licherweise zur Verfiigung gestellt haben, sowie besonders auch Herrn Dr. Tu. Gets, 
der die Programmierung und Durchfiihrung der Rechnungen entgegenkommender- 
weise tibernommen hat, médchte ich auch an dieser Stelle herzlich danken. 


* Fir 72M lautet namlich diese Differentialgleichung 


glV—2a2 Pp’ +a! p—ia Re(y’’—a? yp) =0, 


y=|o?+i«Re. Da Re>o und damit Rr> « ist, ksnnen wir die zweite Losung gegen 
die erste vernachlassigen. Bei H. GOrrLer [22] wird (8.3) aus den Differential- 
gleichungen der reibungsfreien Stromungen gefolgert. 
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Zur Berechnung des kleinsten Eigenwertes j)(c) definieren wir die Itera- 
tionsfolge 


(8.8) a(n) = — J Cli) ali) B84) 4, 
(8.9) bes (n) = J K(q,4) Q(A) ao(f) a9 (A) dA. 


Die Funktion b9¢”(y) wird also zur weiteren Rechnung nur fiir 7<M, die 
Funktion aj}(y) dagegen auch fiir 7=M benétigt. Nun ist nach (8.8) und (3.7) 


(8.10) Hon) Ere firineM 


mit der numerisch zu bestimmenden GréBe 


M 


(8.41) EB = — [© Sino aa(7) HY(i) 7, 
0 


so daB auch a{’}(y) nur noch fiir 7 <M ermittelt zu werden braucht. Wegen 


(8.10) 1aBt sich das Integral (8.9) fiir M<7< oo exakt berechnen. Dann ist 


M 
(8.12) Ys? (n) = J K(y, 7) Q (7) (7) at0(4) aj + 
aS (Ay aE B”) evi (Cx = Bey eo" fiir n = M 


mit 
AO Eo R{ 26% M L 20° o ikea 
40 2o+a (20+a)? 2o+a : 
EY | o o { * 
(Cie ME (2o+a) M 
(8.13) 2 46 R(5 ' “(20+a)2 wera)? ; 
CO TRON Op (20-+0) 
46 2o+a 


Nach Vertafelung der Greenschen Funktion G(y,7) und K(y,%) fir OS, 7M 
mit einer geeigneten Schrittweite / lassen sich die (endlichen) Integrale nach 
der Trapezregel berechnen. Ausgehend von einer Funktion 01}(n) 3:0, die die 
gesuchte Eigenfunktion 0, (7) schon méglichst gut approximieren sollte, konnen 
somit die Iterationen durchgefiihrt werden. Der kleinste Eigenwert ergibt sich 
dann z.B. als Grenzwert der Quotienten bY} (7)/bYo ” (y) fiir »—> co, weshalb wir 
diese nach jedem Iterationsschritt fiir 7 =0(/) M berechnen. Die Streuung dieser 
Werte ist dann zugleich ein MaB fiir die Giite der Naherung. Bedingt durch eine 
giinstige Eigenwertverteilung laBt sich der kleinste Eigenwert schon aus zwei 
Iterationsschritten mit ausreichender Genauigkeit ermitteln. 

Wir haben den kleinsten Eigenwert in sechs Fallen ermittelt, namlich fiir die 
Parameterwerte « =0,5 und « =0,75, jeweils kombiniert mit Re = 2500, Re = 5000 
und Re =6000. Dabei wurde M =1,6 und/=0,1 gewahlt. Die numerische Aus- 
wertung wurde vorgenommen fiir o =0,; 0,3; 0,5; 0,7; 1,0 und 1,5. Die Werte 


des zugehoérigen Gortler-Parameters [p= 2 = Je, Re* eo 0 (0) ‘sind im Abb. 5 


aufgetragen. Sie zeigen, daB der kritische Wert o mit wachsendem « wachst; 
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mit wachsender Wellenlange der Primarstérung wird also der kritische Achsen- 
abstand der sekundaren Langswirbel gréBer. AuBerdem liegt offenbar sowohl 
fiir «—0,5 als auch fiir « =0,75 eine wenn auch schwachere Abhangigkeit des 
Gortler-Parameters von der Reynoldsschen Zahl vor. Wie nach den Ausfithrungen 
des §3 und am Anfang des § 7 zu erwarten war, ist nach diesem Modell Wirbel- 
Instabilitat erst fiir gréBere Werte des Gértler-Parameters als beim Gortler-Modell 


zu erwarten. 7000 
6 


2500 


7000 


500 ie 


0 a5 10 15 20 


Oo —_P 


Abb. 5. Eigenwertkurve nullter Ordnung in Form des Géortler-Parameters u,=2¢€) Re? «? Q(0) fiir verschiedene 
Parameterwerte: 


a=0,5, Re=2500 (Fall1); a=0,5, Re=5000 (Fall2); a=0,5, Re=6000 (Fall 3); 
a=0,75, Re=2500 (Fall4); «=0,75, Re= 5000 (Fall5); «=0,75, Re=6000 (Fall 6) 


In methodisch ahnlcher Weise ]a8t sich das transponierte System (7.15) 
lésen. Wir definieren hierzu ausgehend von einer Funktion 7?) (7), die die Eigen- 
funktion 79 (7) wieder schon méglichst gut approximieren sollte, die Iterationsfolge 


(8.14) so’ (7) = Q(n) a(n) J K(y.m) 10 (9) aH, 
0 

(8.15) 70" (9) = — ao(y) f (9,7) 0 (fi) a9. 

0 

Dann ist nach (8.2) 

(8.16) Asien baton yesh. 

so daB sich (8.14) vereinfacht zu 
M 

(8.17) so’ (n) = Q (4) ao(n) f Ky, i) 78 (H) ay. 
0 


Nach (3.8), (8.2) und (8.3) ist dann 
(3.18) $0 (n) =(EPn + EP) e"t9" fir n> M 
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mit den numerisch zu ermittelnden GroBen 
M 


v 17 1 —On A On 
EP=K [|i {20° + 0} — 7 ef {0} 
(8.19) aie ; 
G I Teo $s Ales As 
EP =K [| fo eifoq +1} + fe fo —1} 
0 


19 (ij) a7, 


19 (7) 47. 


Auch hier sind also beide Funktionen der Iterationsfolge nur noch fiir HSM zu 
berechnen. Dabei ]aBt sich in (8.15) wegen (8.18) noch das Integral fiir M<7< oo 
exakt integrieren. Dann gilt 


M 

(8.20) 78*?(n) = — a5 (n) { [EO 9) 8H) aA +E” Sinoy| fur <M 
0 

mit 


(8.21) Eo —|Ep(_* rey Dera eee 


“s g-(2o+a) M 
2o+«) ' (20+4a)%o 2 “o(20-+4) ; 


Sowohl s§”(7) als auch 7j'*”(y) kénnen dann mit Hilfe der bereits vorliegenden 
Vertafelungen der Greenschen Funktionen ermittelt werden. Auch hier reichen 
zwei Iterationsschritte aus. Wegen (8.16) ist dann 


(8.22) I = 79 (1) 10 (m) 407; 


dieser Wert 1a8t sich wiederum leicht nach der Trapezregel bestimmen. 

In der gleichen Weise lassen sich die in (7.16) auftretenden Integrale berechnen. 
Diese Aufgabe vereinfacht sich noch dadurch wesentlich, daB bei nahezu allen 
Summanden von L, und L, eine Ableitung der Funktion 7(7) als Faktor auftritt, 
die fiir 7=M verschwindet und iiberhaupt nur in Wandnahe zu beriicksichtigen 
ist®. Auch reduziert sich wegen (8.16) die Integration tiber 7 auf das Intervall 
O0<y<M. Wir wollen die Formeln hier nicht weiter angeben, sondern nur mit- 
teilen, daB die relative Korrektur yz, des Eigenwertes in allen durchgerechneten 
Beispielen nie mehr als 4°/o) betrug. 

Es soll nun noch der Grenzfall kleiner o-Werte genauer betrachtet und dabei 
insbesondere gezeigt werden, daB sich der kleinste Eigenwert jv9(o) wie im Falle 
der Taylor-Wirbel wie o? verhalt. Der Beweis schlieBt sich an die entsprechenden 
Uberlegungen in den §§1 und 4 an. Wir machen wieder Gebrauch von der 
Méglichkeit der Normierung der Eigenfunktionen durch E=1 und der Ab- 
schatzung (4.4) 


(8.23) Moly) Se". 


Da weiter nach (4.3) lim H(y,%) existiert, existiert auch 
o>0, 


lim f H(y,%) Q (7) 40(n) %10(%) 47, 


8 Hier allerdings kénnen diese Terme sehr groBe Werte annehmen. Da dieses 
aber nur in einem sehr kleinen Intervall, etwa 0<7< 0,2 der Fall ist, und gerade 
hier die Greenschen Funktionen klein sind, ist dennoch der EinfluB auf die Korrektur 


des Eigenwertes klein. 
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“= 


dh. der erste Summand in der Klammer von (8.7) strebt wie o? gegen 0. Im 
Gegensatz zu den Verhaltnissen beim Gértler-Modell verhalten sich hier jedoch 
auch die restlichen Summanden wie o?, so daB fiir c->0, ein (lésbares) Eigenwert- 
problem mit jim (uo 0) als Eigenwert vorliegt: 


M 
ay (4) = a) Go (7,7) a9 (7) 21 9(%) 4%, 


M Ny 
orn) ee (14002) | f Hon, Q (7) 40) %10(%) 47 + Ane + Bap} far SM 
mit i 
A 1 1 —aM 
A SKE(-+aa)¢ i 


Ebenso bestatigt man dann auch, da8 die durch (7.16) definierte Stérungs- 
korrektur p, sich fiir o+0, verhalt wie o?. 


III. Sekunddre Wirbelstorungen einer durch Tollmien-Schlichting-Weilen 
primar gestorten Schichtenstr6émung 
89. Einfiihrung neuer Koordinaten und Herleitung der Storungsgleichungen fiir 
den Fall B>0 
Wir betrachten nun den Fall einer durch stromabwarts wandernde Tollmien- 

Schlichting-Wellen primar gestérten Grundstrémung. In einem mit der Phasen- 
geschwindigkeit B/x bewegten System — die Koordinaten seien wieder mit 
x,y, z,t und die Geschwindigkeiten mit u, v, w bezeichnet — haben die Grund- 
gleichungen dieselbe Gestalt (5.1), (5.2); dagegen lautet die Stromfunktion der 
wellenformig gestérten Str6mung bis auf Glieder héherer Ordnung in ¢ 

(9.1) p(x, ¥) =yo(v) —2y + e9(y) cos (a * +7(y)). 

Die primar gestérte Grundstrémung hat also bis auf Glieder héherer Ordnung 
in ¢ die Geschwindigkeitskomponenten 


Uy = Uo(¥) — as + eq’(y) cos (x +7(y)) — eg(y)7"(y) sin (ax +7 (y)), 


(9.2 
7 dy =eug(y) sin (ax +7(y)). 


Zur Beschreibung der Stromlinien benutzen wir den (nicht notwendig reellen) 


y-Wert, der im Wellental, genauer fiir x = — f 7(v), angenommen wird. Dem- 
a 

gemaB werden Stromlinien durch konstante 7-Werte?! in 

(9.3) p(*, ¥) = po(n) — S n+ €q(n) = F(y) 


1 Diese Definition von 7 weicht von der des § 6 ab. Wir haben jetzt in dem uns 
besonders interessierenden Wellental 7 = y und erreichen zugleich die Auszeichnung 
des Stromungsbereichs au8erhalb der Katzenaugen durch reelle Werte von 7». Dagegen 
erhielt durch die friithere Definition die gewellte Wand gerade die Gleichung 7 =0 
oe nicht 7 = — e, so daB die Greenschen Funktionen der §§ 3 und 4 benutzt werden 
konnten. 
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gekennzeichnet. In einem gewissen Abstand von der kritischen Schicht gilt 
dann, wie man durch Taylor-Entwicklung nach e verifiziert, 


(9.4) n=y+e— 10) (cos(a x +7(y)) — 1); 
. Ug(y) — ie 


man hat dort also das erwartete Bild wellenférmig gestérter Stromlinien, die 
oberhalb der kritischen Schicht in positiver, unterhalb in negativer x-Richtung 
durchstrémt werden. In der Nahe der kritischen Schicht dagegen ergibt sich auf 
Grund einer Entwicklung nach Potenzen von y—y, bis zu den quadratischen 
Gliedern 


' (e TV 204 (4 —cos(a* + 7)) + (y —y,)?4 


; 
+ F(6q' cos(ax +1) 297 sin(ax +7)) (yy); 


dies fiihrt zu dem in Abb.6 etwas 

vereinfacht wiedergegebenen ,,Katzen- SES 

augenbild‘’. Beim Durchgang durch die 

kritische.Schicht-géhen, wie. schon der -  =——<——~"__ —__.._—""__ 
Vorzeichenwechsel des Nenners von  Abb.6. Stromlinienbild in der Nahe der kritischen Schicht, 
(9.4) anzeigt, Wellentaler in Wellen- bezogen auf ein mit der Phasengeschwindigkeit bewegtes 


Koordinatensystem 
berge itber und umgekehrt. Da aber 
gleichzeitig %) einen Nulldurchgang hat, haben wir z.B. langs der Kurve 


is 1 iy) sowohl oberhalb als unterhalb der kritischen Schicht gekriimmte 
a 


Stromlinien, die konkav sind zur Richtung wachsender Werte der Geschwindigkeit. 
Wir haben also sowohl oberhalb als unterhalb der kritischen Schicht das zur 
Wirbelentstehung notwendige instabile Gleichgewicht. Wir wollen deshalb auch 


die Umgebung z.B. der Kurve * = — -j (vy) durchweg als Wellental bezeichnen. 


In der Nahe der kritischen Schicht kénnen die in (9.1) vernachlassigten Glieder 
hdherer Ordnung in ¢ einen gewissen EinfluB haben. Ebenso werden hier bei der 
uns im Moment allein méglichen Zugrundelegung einer linearisierten Tollmien- 
Schlichting-Welle die durch die analytischen Rechnungen bedingten Glieder 
hdherer Ordnung in ¢ von Bedeutung sein. Wir fiihren deshalb diese Rechnungen 
bei Zugrundelegung von (9.1) vorerst in Strenge durch und definieren somit die 
Vektoren ¢,, €y, ¢; des auf die wellenférmig gestérte Strémung bezogenen Drei- 
beins durch 


By Pay 2 TE ERI 
Vero a Yerer?’ 
2 Waa + 8 +ViteR’ 
Cale 


Dabei ist eR aac! nach (9.2) eine GréBe erster Ordnung in ¢ und + \1 + @2 R? 
U 


9 . . 
ist positiv (negativ) zu nehmen fiir %)>0 (#%)<0). Mit Ausnahme der Schnitt- 
punkte der singuldéren, die Katzenaugen berandenden Stromlinien, in denen 
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tip =d)=0 ist, sind die Vektoren e,, ey, es eindeutig definiert. Beim Durchgang 
durch die kritische Schicht auBerhalb der Wellentaler drehen sich e, und e, stetig 
mit der Neigung der Stromlinien, wechseln also hierbei im wesentlichen ihre 
Orientierung. 

Aus (9.6) folgen fiir die Umrechnung der Komponenten des Geschwindigkeits- 
vektors auf Grund von (6.7) die Beziehungen 


= —____——__ , = —_____—_. c=—wW. 
Pe i +)1i+ eR? +)i+eR? 
Fiir die Grundstrémung (9.2) gilt damit 
(9.8) dy = bi V1 + PR? =VR+05, bo =0, cg =0. 


Da ¢,, €3, ¢; orthogonale Einheitsvektoren sind, haben die Ableitungsglei- 
chungen eine sehr einfache Gestalt. Fiir die Differentiation nach ¢ etwa gilt 


dR : dR 
5 ee 
(9.9) dey dt me dey _ dt e, des __ 
dt 4-+ 6? R? dt 1+ 6? R? dt 


Damit lassen sich die Navier-Stokesschen Gleichungen wie in § 6 beztiglich des 
Dreibeins (9.6) zerlegen und die Stérungsgleichungen zur Grundstrémung (9.8) 
ableiten. 

Die Beschleunigungs- und Druckanteile kénnen wir noch durch Einfiihrung 
einer auf den wellenformigen Stromlinien konstanten Koordinate 7 vereinfachen. 
Wir benutzen hierzu die durch (9.3) definierte GréBe 7 =F 4 (p(x, y)), wobei F4 
die Umkehrfunktion der durch (9.3) definierten Funktion F(7) bezeichnet. Weiter 
ist wegen der lokalen Betrachtung der Wellentaler wie in § 6 die Einfiihrung der 
GroBe € als zweite Koordinate zweckmaBig. Demgema8 benutzen wir als neue 
Koordinaten 


§=x+—4(9), 

H= Ft (vo (9) = By + €q(¥) Cos (% x +i()), 
(9.10) 1 

cee, 

(oie 


Die Umrechnung der Ableitungen erfolgt dann gemaB 


n= amne a) 

On. VOR ee 

1d: Dyes te tae Re 
(6.11) aah og Soe ea ae 

oe 

yee. 

Oye 

“Oten Ove 


Dabei ist F’ =u,— ze +e q'=0 offenbar gerade auf den die Katzenaugen beran- 


denden Stromlinien, da auf diesen 7 konstant und in deren Kreuzungspunkten 
F"=%)=0 ist. Nach (9.11) werden also diese singularen Stromlinien der ebenen 
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wellenférmigen Strémung zu Singularitaéten (z.B. Trennflachen) der sekundaren 
Langswirbel werden. Insbesondere wird auch deren Schnittpunkt bei lokaler 
Betrachtung der Wellentaler eine singulare Stelle sein, was schon nach (9.6) zu 
vermuten war. 

Wir wollen die Umrechnungen und die linearisierten Stérungsgleichungen im 
einzelnen hier nicht angeben, da keine prinzipiellen Anderungen gegeniiber den 
Ausfiithrungen des § 6 vorliegen; wir wollen nur bemerken, da es auch bei voll- 
standiger Beriicksichtigung der -Abhangigkeit die Reibungsterme sind, die durch 
das Auftreten singularer Anteile die wesentlichen Schwierigkeiten bedingen. 
Vernachlassigt man wie in §7 in den Wellentadlern die Abhangigkeit von der 
Koordinate &, so liefert der Ansatz sekundarer Wirbelstérungen 


a= 4a,(h) cose ec”, 
b =b,(n) cosale”, 
(9.12) ‘ 1 (7) C 
C= (7) Sino €*, 


p = br (m) cosa e”" 
entsprechend (7.2) die Stérungsgleichungen 


ya, + tid, + 209d, = 2 [Lila] + L,[6,]], 


Oe Bea gue: 640) b—= —p, Pe |— Lala.) + Lyell, 


Re 
(9.13) , 
SERA Ma ree Male 
2 if 
ye Ppa bl 4 oe =O 
Ug 
mit 
e2 g2 2 (2+ 172 
I, [a] = Aa Ee Z ) a, 
0 
ab Ree 32d EG eee 
14) a . . : ? 4 
ve 2eP a9" q? | 2eag’y’teaqs’ , (woteg’—eqs*)eagy |p 
‘ia a tho a 2 
2 172 a 
Aa, = ay — e* FF ga; — oa. 


0 


Diese sind unter den Randbedingungen 


(9.15) ee Ge) EU by) —— On Ind ie 100 


zu losen. 


§ 10. Erste Lésung des Eigenwertproblems der Sekunddrinstabilitat 


Der Ansatz (9.12) sekundarer Wirbelst6rungen in den Talern der Tollmien- 
Schlichting-Wellen fiihrt gemaB (9.13), (9.15) wieder auf ein Eigenwertproblem 
mit der Wellenamplitude «¢ als Eigenwert. Dieses unterscheidet sich von dem- 
jenigen des § 7 hauptsdchlich dadurch, daB einige Koeffizientenfunktionen der 
linearen Differentialgleichungen Pole an der Schicht #)=0 haben und deshalb 
eine Singularitat der Lésung an dieser Schicht bewirken. 


19* 
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Vergleicht man zundchst die Gleichungen (7.2) und (9.13) miteinander, so 
sieht man, daB auch hier bei den Beschleunigungsgliedern die Terme auftreten?, 
die frither die Glieder nullter Ordnung lieferten. Alle anderen Glieder, insbesondere 
auch die mit den singularen Koeffizientenfunktionen sind mit dem Faktor ¢ be- 
haftet und werden somit nur in der Nahe der Schicht # = 0 zu beriicksichtigen sein. 

Eliminiert man deshalb wie in §7 die Funktionen #, und c, und ersetzt 
Re }, durch b,, so sind die Gleichungen fiir neutrale Wirbelstérungen wieder 


von der Form 
Al ae 
Gy ON 8 alee Ole 


10.4 E 
ie) bY — 207 by + 04d, = — po? qa, + eo? L,[(%, by], 


(10.2) Gy = 0, = 0; =O) 1 ey Oat eo: 


Dabei ist analog zu (7.5) zur Abkiirzung « =2¢ Re? gmyax bei durch qyax=1 
normiert gedachten g(7) gesetzt worden; L, und L, bezeichnen wie friiher lineare, 
jetzt aber in der Nahe der Schicht #)=0 singulare Differentialausdriicke in a, 
und b,. Die in diesen Restgliedern zusammengefaBten Anteile brauchen nur in 
der Nahe der kritischen Schicht beriicksichtigt zu werden. Wir entwickeln des- 
halb hier die Lésung (in Abhangigkeit der 6 Integrationskonstanten) in eine Reihe 


nach Potenzen von (n —y, +e S ). Fiir kleine bzw. groBe Werte von 7 dagegen 


/ 


Uo 
legen wie die wie friither vereinfachten Gleichungen zugrunde und erhalten hier 


durch Reihenentwicklung um 0 bzw. durch strenge Lésung fiir 7>M und an- 
schlieBende numerische Integration unter Beriicksichtigung der Randbedingungen 
(10.2) fir 7=0 bzw. 7 co Naherungslésungen, die jeweils noch von 3 Integra- 
tionskonstanten abhangen. Schlie8t man diese 3 Losungen an geeigneten Zwischen- 
stellen zusammen, so erhalt man 2 Systeme von je 6 homogenen Gleichungen 
fiir die insgesamt 12 Integrationskonstanten und damit eine Eigenwertbedingung 
fiir die noch zu bestimmende GroBe w bzw. «. In Anbetracht des Arbeitsauf- 
wandes konnte dieser Lésungsweg numerisch noch nicht durchgefiihrt werden?. 

Wir miissen uns deshalb hier mit der Angabe einer Naherungslésung begniigen, 
bei der die Anteile eL, und cL, vernachlassigt und die Singularitaten des Strom- 
linienbildes somit nur tiber das benutzte System (9.6) der Basisvektoren beriick- 
sichtigt werden. Dieses erscheint bei unseren numerischen Erfahrungen bis zu 
einem gewissen Grade gerechtfertigt. Das somit zugrunde liegende Eigenwert- 
problem ist mit (7.10), (7.142) nach Ersetzen von Qa, durch q identisch und kann 
entsprechend den Ausfiihrungen des §8 gelést werden. 

Da auch hier q(7)>0 ist® fiir 0<7< o, hat das so vereinfachte Eigenwert- 
problem nach R. JENrzscu [15] wieder einen betragsmaBig kleinsten positiven 


‘Man hat nur Qay=q zu setzen, um auch formale Ubereinstimmung zu erhalten. 

* In Strenge erfordert dieses Problem sogar ‘die Beriicksichtigung der nicht- 
linearisierten Storungsgleichungen, da sich nach (10.1) eine Unendlichkeitsstelle der 
Losung der linearisierten Stérungsgleichungen ergibt, sowie die Beriicksichtigung der 
Nachbarschaft der Stelle £=0 wegen der Unbestimmtheit des Koordinatensystems 
in der kritischen Schicht fiir = 0. 

, 3 Wie erwahnt, wechseln in der kritischen Schicht die Stromlinienkriimmung und 
die (relative) Geschwindigkeit ihr Vorzeichen, so da8 sowohl oberhalb als unterhalb 
dieser Schicht eine instabile Situation vorliegt, was sich in ¢(y)>0 und a(n) >o 
ausdriickt. : 


Instabilitat laminarer Stromungen 281 


Eigenwert, fiir den die zugehérigen Eigenfunktionen a, (7) und 0,9(7) ihr Vor- 
zeichen nicht wechseln. Da aber die Basisvektoren beim Durchgang durch die 
kritische Schicht gemaB Abb. 6 im wesentlichen ihre Orientierung dndern, besagt 
dies, daB die entsprechenden Eigenfunktionen u,9(7) und v,9(y) hierbei ihr Vor- 
zeichen wechseln. Dieses wiirde dem Auftreten zweier durch die kritische Schicht 
getrennter Langswirbel mit entgegengesetztem Drehungssinn entsprechen. 


Abb. 7 enthalt die so erhaltenen Kurven des kleinsten Eigenwerts fiir die 
beiden von H. SCHLICHTING [19] sowie zwei von J. A. ZAAT [23] berechneten 
neutralen Wellenst6rungen. Ihre Daten sind Re = 2620, « =0,466, B/x =0,350; 
Ke =6080, « =0,737, B/a=0,350° bzw. 
Re =1465; "a= 0,047, Bla=0;372 und 
Re =3095, «=1,072, B/«=0,372. Die erst- 
malig zur Sekundarinstabilitat fithren- 
den Amplituden sind in der GréBen- 
ordnung von 104 Grenzschichtdicken. 
Die erste dieser Stérungen hat auch i 
G. B. SCHUBAUER [24] seinen Messungen Jo 
zur Klarung des Mechanismus’ des lami- 
nar turbulenten Umschlags zugrunde ge- 
legt. Wie bei allen theoretischen Stabili- 
tatsuntersuchungen ist der errechnete 
kritische Wert des Stabilitatsparameters W00 
(hier um einen Faktor 10) kleiner als der ae ae yap id 
experimentell beobachtete ; dagegen be- Abb. 7. Eigenwertkurve des vereinfachten Systems in 
steht eine gute Ubereinstimmung (nam- Form des Gértler-Parameters Mo= 26 Re? o? gy fiir die 
lich von etwa 10%) des errechneten kriti- ee ee ee rae 
schen Achsenabstandes unserer sekun- 
daren Langswirbel mit dem Abstand der StraBen hdherer Energie, die sich auf 
Grund der Messungen von G. B. SCHUBAUER nach dem Auftreten der Tollmien- 
Schlichting-Wellen in Stromungsrichtung bilden. 


J00 


250 


750 


Zusammenfassung 


In dieser Arbeit wurde gezeigt, daB eine laminare Strémung mit wellen- 
formigen Stromlinien in den Talern dieser Wellen instabil ist gegentiber St6rungen 
vom Typ der Taylor-Gértlerschen Langswirbel, falls nur die Stromlinienkriim- 
mung geniigend groB ist. Die dieser kritischen Kriimmung entsprechende Wellen- 
amplitude lag bei der wellenférmig gestérten Blasiusschen Plattenstrémung in 
der GréBenordnung von 10~4 Grenzschichtdicken. 

Besonderes Interesse hat diese Fragestellung im Hinblick auf den laminar- 
turbulenten Umschlag. Nachdem W. ToLLMIEN [7] die Instabilitat der laminaren 
Plattenstrémung gegeniiber wellenformigen Stérungen oberhalb einer gewissen 
Reynoldsschen Zahl gezeigt und dann G. B. SCHUBAUER und H. K. SkRAmstap [8] 
diese Wellen als Anfangsstufe des laminar-turbulenten Umschlags erkannt 
hatten, ist die Frage des weiteren Ubergangs in die turbulente Strémung 
viel diskutiert worden, ohne jedoch bisher in befriedigender Weise beantwortet 
zu sein. 
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Nach den Resultaten dieser Arbeit liegt die Vermutung nahe, da8 die wellen- 
formig gestérte, aber immer noch ebene laminare Stromung auf Grund ihrer 
Instabilitat gegentiber den wirbelartigen Sekundarstérungen in eine dreidimen- 
sionale laminare Strémung tibergeht und der Umschlag in die dreidimensionale 
turbulente Strémung schlieBlich als Folge von Instabilitaten dritter und héherer 
Ordnung erfolgt. Eine letzte Entscheidung hieriiber kann nur das Experiment 
liefern. Es wird deshalb besonders auf die erst kiirzlich veréffentlichten neuen 
experimentellen Untersuchungen von G. B. SCHUBAUER [24] hingewiesen, die mit 
unserer Vermutung in einem fiir den gegenwdrtigen Stand von Theorie und 
Experiment durchaus befriedigenden Einklang stehen. 
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Expansion and Boundedness Theorems for Solutions 
of Linear Differential Systems 
with Periodic or Almost Periodic Coefficients 
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1. Introduction 
In this paper systems of the form 


aw _ 
ia 


(1.4) —K Aw +AB(z)w 

are studied, where the independent variable is real or complex, and the variable 
coefficients which are the components of the matrix B(z) are periodic, almost 
periodic, or have general non-harmonic Fourier expansions 


(1.2) B(z) = By +> Be 

kE1 
with arbitrary complex exponents 70,. Solutions are obtained which are gener- 
alizations of those for periodic systems established by the FLogueEt theory [3]. 
They are of the form 
(1.3) (Zee ee 

Meek 
where the w, are linear combinations of the @, with non-negative integral coef- 
ficients. Explicit expansions are obtained both for the coefficients c, and for 
the determining equation (in general transcendental) whose roots are the charac- 
teristic exponents T. 


The method by which these expansions are obtained and by which their 
convergence is proved is essentially that developed by Wasow [13] and Go- 
LOMB [6] in connection with non-linear systems having small forcing terms similar 
to the components of B(z). The formal procedure involved had been suggested 
previously by E. WEBER (see Wasow [73]). These expansions (an appropriate 
name for them would be “‘free parametric expansions’’) are quite different from 
the expansions in powers of A and, at least in form, also from various versions 
of successive approximation expansions with elimination of secular terms that 
have been used extensively for systems of form (4.1) with periodic B(z). The 
literature on this subject is so vast that the references in this paper are necessarily 
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limited to a few specific items most closely related to it. For the general back- 
ground the reader is referred to CESARI’S monograph [2]* with its unsurpassed 
bibliography. 

One advantage of the method used here is that it leads in many cases to less 
formidable calculations than other existing methods. Even in the case of the 
classical Mathieu and the general Hill equation, for which very elaborate theories 
exist [see, e.g., 10, 11, 12, 14|, its application seems to yield some unknown 
expansions (Sections 4,5). The method is also applied to systems of second 
order of the form 


2 m 
(1.4) Ft 4 GB tty +24 YB y (2) My =O, eg Wee) 
v—-1 

which (with the restriction to periodic b,»(2)) were studied by Cxsart [J], 
HAACKE [7], HALE [8, 9], GamBILt [4, 5], and others. CEsARI observed in 1940 
that for systems (1.4) with m>1 various conditions of symmetry of the coef- 
ficients 6, assure stability while, when no condition of symmetry of the coef- 
ficients b,,, prevails, instability may well occur. These studies were carried on 
by Crsari himself and by his students HALE and GAMBILL, who also applied 
CESARI’S method of successive approximations to nonlinear systems. 

As with the other general methods, convergence of the expansions (1.3) is 
proved only for sufficiently small 4 (Section 6). As with the method of CEsart, 
it is not essential that the parameter 4 appear as it does in equation (1.1). Since 
the expansions are not in powers of A, the system 


(1.5) 1K 2A) +B2)u, 


where A(A) tends to A and B(A, z) tends to 0 as A approaches 0, can be treated 
by the same method with little change. However, this study is limited to the 
model system (1.1). It is also assumed that the expansions (1.2) for B(z) are 
absolutely convergent. This condition, too, can be relaxed for special classes of 
B(z), as will be shown in another paper. In the case of periodic B(z), CESARI’S 
method of successive approximations [7] was shown by HALE [8] to converge 
for L-integrable B(z) as well as for those with absolutely convergent Fourier 
series. The one essential condition in the convergence proof is that the difference 
of any two characteristic values o, (which are the values of o that make the 
matrix ok —A singular) have a positive distance from the “compound spec- 
trum” of B(z) which is the set of linear combinations with non-negative integral 
coefficients of the exponents 9,. Because of this condition there are no “‘small 
divisors” in the expansion of the solution. For periodic B(x) it is implied by the 
seemingly weaker condition that none of the differences o,—o, belong to the 
compound spectrum of B(x), which is also a necessary condition in CESARI’S 
theory. 

The characteristic exponents 7t,, lie, for small A, near the numbers zo,. In 
Section 7 the determining equation whose roots are the t,, is examined, and it 
is proved that there are at least as many linearly independent solutions of (1.1) 
of form (1.3) as there are distinct characteristic values o,. If the @, in (1.2) and 


* The manuscript of this monograph was available to the author before its 
publication. 


286 MiIcHAEL GOLOMB: 


the t, are real, then the corresponding solutions are bounded on the real line 
and system (1.1) is (at least restrictively) stable. Conditions which derive the 
reality of the roots t, from various reality and symmetry conditions on the 
matrices K, A, B(x) are given in Section 8 (Theorems 8.1 through 8.5). Several 
of these were found among others by CESARI, HAACKE, GAMBILL and HALE in 
the papers cited above. Finally, a general criterion for boundedness and un- 
boundedness of all or some of the solutions on the positive real line is established 
(Theorem 8.6). 
2. Formal Expansions 


The systems considered are of the form 
Ne Oey 
(2.1) ~ Ka =AwtsB)y. 


The independent variable z is complex, the unknown w is a vector whose 
components are functions of z. The matrices K and A are constant, 4 is a com- 
plex parameter, and B(z) is a matrix with the expansion 

(2.2) B(z) = By + Dd, Be, 

kel 

where the B,, B, are constant matrices and the o, are complex numbers different 
from 0. It is assumed that there is a part 8 of the complex z plane including 
z=0 and z=1 in which the series (2.2) is absolutely and uniformly convergent. 
Since 8 may be assumed to be convex, 8 either contains interior points and the 
components of (2.2) are analytic in %, or B is an interval of the real axis including 
0 and 1. In the first case the components of the vector w are sought as func- 
tions analytic in 8; in the second case z is considered as a real variable and the 
components of w are to be differentiable functions. 


Let P denote the sequence 


(2.3) P= (015 2, ++) 
and 3(P, 8) the class of functions f(z) possessing an expansion 
(2.4) f(2) = ay + Dag ee 

kE1 


absolutely convergent uniformly in 8. We look for solutions w(z) of system (2.1) 
whose components are of the form f(z) exp(itz) where t is a complex number, 
fEB(Q, B), and the elements of Q are linear combinations of 0,, @),... with 
non-negative coefficients. 

Let the sequences of non-negative integers only a finite number of which 
are +0 be arranged as a sequence N,, N,, Ny,... in such way that N+N=N, 
implies tr +s (see [6]). Here N+ N, denotes the usual sum of the sequences 
N,, N,, and Ny denotes (0,0, ...). Also let 


Na) = (ay OF Oar) 

(2.5) Nosy = (0, 1, 0, £2.) 

and put 

(2.6) N,- P=) m, On, |N,| = 2) 
kE1 kel 
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for N=(m,2,.-.). Then 
7) 0, = Na: P, [PVN ass 


Some of the complex numbers N, - P may be equal to 0. The set of indices 7 
(including 7 =0) for which N,- P =O will be denoted as 2. We set up 


(2.8) (2) = (2, + >} a, en 4 gue, 
rgEm 


where the constant vectors a), a, a,, ... and the number 7 are to be determined 
so that (2.8) satisfies (2.1). For this purpose we put 


(2.9) (NP +2)K—A)a,=AD Bay, Not Ny=N, 7M, 
(k, k’) 


where the sum is to be extended over all those k, k’ for which NaytNy=N,. 
In particular, for y=(k) formula (2.9) gives 


(2.10) (G+) K —A)agy=AB, a, k=4,2,.... 


Since k'<yv, system (2.9) is a recursive system for the a,,7¢M. If we now make 
the assumption that 


(2.11) (N,-P +1t)K —A_ is non-singular, FON, 
then (2.9) determines the a, uniquely as linear functions of dp: 
(2.12) GP TAA tay vEN. 

Next, let us assume that the sum 
(2.13) ZB AMLD, (Myth) P=0,  r€M, 


extended over all those r¢9 and 7’ for which (Ni,+N,) - P =0, is convergent. 
Then to complete (2.8) as a solution of (2.1) we set 

(7,7) 
This is a linear homogeneous equation for a). It has a non-trivial solution if and 
only if 
(2.15) F(A, 1) = det [cK Aa Be YBa 7)) =. 

(1,7) 

This is the determining equation* for the characteristic exponent t. If this equation 
is satisfied, if a) is determined from (2.14) and a, (r¢%) from (2.12), and if the 
two series 


(2.16) 2A, (A, ) eiNe Pe AE t) et (Ne+N(s))* P2 
s2l1 


converge uniformly and absolutely in %, then (2.8) is easily verified to be a solu- 
tion of system (2.1) of the desired form. 


x The expressions ‘“‘determining equation’”’ for the equation whose roots are the 
characteristic exponents was used before by GAMBILL & HALE in their paper [15], 
and probably by others. 
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For later application we determine the structure of the 4,(A, 7) (v€M) in 
some detail. Put 


(Q47) C=C) SNe ae) =: rN. 

By (2.9) for 7¢ 

(2.18) js AC ete Noy + Ny =N,. 
(k, k’) 


For any k’-+-0 the same formula may be applied again, and repeating this until 
no a, with a subscript k=+0 is left, one obtains 


(2.19) a, = (2 CD Bs pie oped eae Oot B,) Qo, 
where 
Noy +My to tM =N,, ner, 
(2.20) Nos) a viel a Nis) =a vee 7 EN, 
Nosy = N,,> r,EN, 


and the sum in (2.19) is extended over all 7,,...,7;, 5;,.--, 8; (f21) for which 
(2.20) holds. From the first of the equations (2.20) it is seen that 


(2.21) #=N,|- 

By (2.12) and (2.19) we have 

(2°22) Vie Wiel hike RUIN Cie g SO) 2h Dre Re 
(7,8) 


and (2.15) becomes 
(223) FAs) dete A A Sr Ca Bo: 
(7,8) : 
In (2.23) the sum is extended over all indices 7, ..., 7;, 8, ..., 8; for which (2.20) 
holds, and sp is restricted to those values for which 


(2.24) (Nis) | Nis) | avers i Nis,)) . P — O. 
This sum will be denoted as 
(2.25) SoA, tae BOR deta Asp Dy. 


(7,8) 
The results are summarized in 


Theorem 2.1. Assume the matrix B(z) of the system 


1 7 dw 
has an expansion 
(2.2) B(z) = By +> B, ei Nw Ps 
kE1 


absolutely convergent uniformly in B. Assume C,=(N,-P+1)K—A its non- 
singular for N,. P==0 and the two series 


(2.26) R (A, T, 2) a pS A Ox BS tee Ce Be er Nee $2 
5) 
(2.27) S(A, 0; 2) = D1 B,C" B,... Cz? B,, ent Ne) Pe 


(r,8) 
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are absolutely convergent uniformly in 8. (Here 1,,...,%,, S1,-.., 8; extend over 
all solutions of system (2.20) and s ranges over all positive integers.) If ay is a 
solution of the linear system 


(2.28) . [tK —A—AB,— # Sy (A, t)] a, =0, 

where So(A, t) is the sum (2.27) with s=Sp restricted by (2.24) then, 
(2.29) w(z) =[1+AR(A,t; 2)] e** a, 

1s a solution of (2.1). 


For the important special case when B(x) is periodic of period 2z/w (w>0), 
equations (2.20), (2.24) take a simpler form. We then may put 


(2.30) P=(0, =o, 20;,—20=.>.). 
Each N,,,)- P is a positive or negative multiple of w; we set 


Ney: P= = 41, Uae pret A 


(2.31) 
By, = B (ep +1): 
Then by (2.24) 
(2.32) My + Xe o> + Hy, = 0, 
and by (2.20) 
(2.33) N,,° P = (oy ++ Hg + +°: + %,) @ = 0, [Shae WP od fe 
We also set 
GAIN, ol aT el ed 
(2.34) ={(q tx+--+%,)o+}K—A 


= C(x, +H +--+ +), — 1 BO NS i 
Then (2.25) becomes 


(2.35) So(A, 2) =>, A B (v4) C3 (54) B (249) C2 qt). CA qt +) Bg 1), 
() 

and the sum is extended over all systems of positive and negative integers ~,, 

Hg, +++, %444 (f21) for which, by (2.32) and (2.33), 

(2336) 2 Me oe ppp HO, Hy Me HO, [SSH ie 


The series (2.26) becomes 
(2.37) R(A, 32) =>, Al? CU xe) B(oqy C2 (xe +04). CA pag +++ p21) Bre,) 4, 
(*) 


where the sum is extended over all systems of positive and negative integers 
Ht, %y,...,%, (f21) for which 


(2.38) w#wtmte-ty=0, *+%4+-°:°°4+%,+0, {Si1sf—1. 


3. Systems of Second Order 


The system of second order equations 


a? m 
(3.1) at + of, u, 4 21 >" Buy (2) 4 = 9; Hig A Onin Me. — 0, == 0, 


v=1 
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where b,,€(P, 8) (u,v=1, 2,.--,m), special cases of which were studied by 
Crsart [1], HAackeE [7], Haze [8, 9], GamBiLt [4, 5], is easily reduced to form 
(2:4) Set 


3.2) Oy Wey a= — Fy MypttMy, Fy Wey = — Fy Myp—tU,, = 1,2,...,M, 
and let w be the vector whose components are W,, W2,..-,Wg,- Then (3.1) 
becomes 
1 Ci = 
(3.3) KK =Aw+ABQ)w, 
where 
(3.4a) K = diag (01, — 01, +++) ms — Om): 
(3.4b) A = diag (Gf oR e.0,,.05)) 
by eb 
‘ pees) Wr ed he a ee 
(3.4c) ee i ! 


CESARI et al. consider the special case where the 0,,(z) are periodic of the 
same period, say 22/w. In this case P=(w, —@, 2m, — 2m, ...), and the numbers 
N,- P are integral multiples of w. The matrix C, of Theorem 2.1 is a diagonal 
matrix with the elements o,(xw +t—0,), —0,(xw +1T+0j),...,6,,(x@ +T—O,,), 
—o,,(%@ +t+0,,), Where x is an integer. The condition that C, be non-singular 
for 7¢% requires in this case that o,--0 and tto,+xw (u=1, 2,...,m; 
~=+1, +2,...). The latter conditions will be satisfied for any t near oj, 


Opa On tt 


(3.5) O,+=0, 6,+0,—- 0, FU a’ ere |e eae ee 
These are among the conditions used by CEsart. The condition o,--o, for w=» 
is not required here. 
4. Mathieu’s Equation 
A special case of system (3.1) is the Mathieu equation 


(4.1) eo +o7u + 2A(cos 2z2)u=0, 

which on substitution of 

(4.2) M=—ou+iw, we=—au—in' 

becomes the system 

Sy * on = 07 w, + A(cos 22) (w, + wp), 
= ae = 07 w, + A(cos 22) (w, + ws). 


In this case P =(2, — 2), and condition (3.5) becomes 
(4.4) QP See. MO dl oom 3 


To obtain the expansions (2.35), (2.37) for the Mathieu equation, we first 
infer from (4.3) that there are only two matrices B, different from 0, 


ae i 
(4.5) B= By =(4 )=2 
BOF 
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which correspond to Ny)- P=2 and Nay: P=—2. The matrix C, of (2.17) is 
diagonal; its elements are o(2% +1—o), —o(2x+1-4 a), where x is some integer. 
Thus 
f A j res i res ie a mab 
(4.6) ong ee (2% +1 — 0) (2% + t— @) 
2¢\— (2x +1+6)1 —(2x%+7-+0)1 


)=Fe, ay Oe 


Since all the B, are 0 except B, and B,,, we conclude that the only integers x; 
used in (2.35) are +4 and —1. Thus, using (4.5) and (4.6) we obtain 


(4.7) SoA, 2) = DW Fe) Pg 2) OF a eee Le 
(4) 
where each xz, in the sum is +1 or —1 and all systems %, ..., x; (f21) are used 
for which by (2.36) 
(4.8) Hy Het +20, tet J 


It follows that # must be odd, and we write (4.7) as 


(4.9) Slt, te © A a) Fog 4-2) & FO ee, eo Ho j—1)- 


( 
By (4.5) and (4.6) 3 
(4.10) E F(x) = [(2% 4+ 1)? —o?|*E. 
Hence, (4.9) may be written as 
(444) SAL) = 2 RI (2%, + 1)? — 07] 4 [ (2%, + 2%, + 1)? 0?]4 
wee [(22ey + 2% >>> + 2%,_, + 1)? — 0]. 
If (4.11) is substituted in (2.28), the determining equation is obtained 


eee 8 ah Ce et as eer tot 
7. (4) 

we+ [(2% + 2% + +++ +2%5;-1 + 1)? — 67], 
or in developed form 


1 1 
(4.13) BAP treat (“24a a8|T 
+4] “lee 2 - ae + (2°). 
((2-- z)\e—— 0)" ((44-2)*— 07). ((—2-42)*= 07)? ((— 4-7)? 6") 


From this equation one derives after some calculation 


=e 1 —8+350%—1504 4 Z 
(4-14) oe +|o+ 4a(1— 0?) ae 64 03 (1 — a?) (4 —o?) Beg: | 


as the two characteristic exponents, exact mod 8. 


If t is equated to 0,1, 2,... in the determining equation, one obtains the 
relation between o and 4 for which the resulting Mathieu equation has the 


periodic solutions cey(z, A), ce, (z, A), ... and se,(z, A),..... For example, for t=0, 
equation (4.14) yields readily 

B=. 21 90 744 6 
(4.15) o s+ ae + (a8) 


(compare [10, p. 15]). 
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a 


Equation (2.26) for a) gives immediately 
(4.16) ESM UG Boos a ed) 
With the two values of t from (4.12) equation (4.16) represents two linearly 


independent vectors except when t=O (Mathieu equation of order 0). To 
calculate the complete solution (2.29), we use (2.37) and (4.6). We find 


(4.17) R(A, 132) =D MF (e) Fe +m)..-Fe+44+---+%)P™, 
(*) 


where x is an integer +0, each of the ~,...,%; (721) is +1 or —4, and the 
sum is extended over all such systems for which, by (2.38), 


(4.18) eb ye ge ed, ey a 0, 1 Ss ee 


By (4.6), for any integers %,, x, 
(4.19) F (1%) F(x) = [(2%2. + 1)? — o?} F(x). 
Hence, (4.17) may be written as 
wey RA, v3 2) =D 7 Wee Et) 074 
Sete [(2% + 2% +--+» + 2x,+ 1)? —o | F(x). 
Using (4.16) we have 
(4.21) F(x) ay = ((2% +14 —o)4, —(2x+4+0)}). 


On substitution of (4.20) and (4.21) in (2.29), the following expansion of the 
solution to MATHIEU’S system (4.3) is obtained: 


4 im or 
W) (2) iene (2 +") +a3/ Q2e7—o) Jax 
(4.22) Ws (z) Or (x) \— (2 +17 +0)7 
: x [(2%+ 2% + 1)?—o7] 4... [(2x+2%,4-- +2%,+ t)?—o?]4 a ; 
Forming (#,-++W,)/20 yields the following solution to MATHIEU’s equation (4.1) 
u(z) = én {t +A >) Ai [(2% + 1)? — 0?) [(Q% + 20 +1)? —o?]2 
(*) 


(4.23) | 
eer ea fs | 24 Tae at aact: 


7 


Replacing t by —r in (4.23) results in another solution of (4.1). From these 
two solutions the following two are easily derived as linear combinations (denoted 
here as cu(z, a), su(z, d)): 


cules) eos te 


+4 2) Ai [(2% + 1)? — o?]4 [(2% + 2%, + 7)? 0?]2 
(x) 5212 08 (2% + 1) 2 
sin(2x +1) z 


(4.24) su(z,4) =sinrz 


we. [(2% + 224 bss + 2,4 7)? 


For t= (n=O, 1, 2, ...) this results in the following novel expansion of the 
Mathieu functions ce,,(z, A), se, (z, A) 


CO;(25A) == COS 


+A DA [(2% +n)? — 07] [(2% + 2%, +n)? o?] 1 
% cos (2% +n) z 
sin (2% +) 2° 


(4.25) se (2, A) == sin n z 


wo. [(2% + 2a, +. 4 2x, +n)? — 0%] 
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The following result is noteworthy. Put 
(4.26) cu(z,A) =costz +4) c,cos(2x +1) z. 
|x|=1 


Then by (4.18) and (4.24) 
(4.27) 3 (G4 +e4) =AD A [(2% +1)? 08]... [(2e +2%,+---4 2%; + t)2=67|74 
(*) 


with x ranging over the values +1, —1 and the x, restricted by (4.18). With 
these choices for x, equations (4.18) are the same as equations (4.8) and the sum 
in (4.27) is the same as that in (4.12). Thus 


(4.28) 


This is a known result for ceg(z, 4) where t=0, 1 =c_, (see, e.g., [14, p. 410)]). 


5. Hill’s Equation 
The general equation of HILr 


(5.1) ae + ou + 2A >'a,,(cos 2m 2) u =0 


nZa1 


becomes on substitution of (4.2) the system 


noe ay AG Os daze (o,, cos 2 z) (1 + Ws) xi 


1 az n> 
(5. ) o dw W. w =0 
pny 
: Z== 02 WD, ae Bl pabes 2 2) ( 1 2) , 


to which we apply the developments of Sections 2, 3 in much the same way as 
was done in the preceding section for MATHIEU’s equation. In the present case 
P=(2, —2,4, —4,...), and condition (3.5) is again 


(5.3) ; CU We, aot == OF. 

The matrices B,, B, which correspond to N,)- P=2s, Nv); P=—2s are 
dy ae 

(5.4) B= By=o.( 3 i )=ok 
2 2 


with E as in (4.5). Thus if the diagonal matrix C, has the elements o (2% +7 —o), 
—o(2%-+1t-+0o) (x a positive or negative integer) we have 


| (2x +t—o)t peli a ee) 


f —1 —C- Boa 
De ee es rae (2% +7—o)2 —(2%+7+0)2 


20 


with F(x) as in (4.6). 
The matrix S)(A, t) of (2.35) becomes 
(2G) Wn glAr i= Le Ani O pd is Doey nal Lea) (oq + og)... F(tq + a+ +++ +), 
(#) 
where each x; is a positive or negative integer and the sum is extended over all 
systems (%,, %:,---,%;41) (7=1, 2,.-.) subject to the conditions 


(527) My he Hy = 0; ty +H +++ +4, 0, (SUES 
Arch. Rational Mech. Anal., Vol. 2 20 
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By use of ie equation (5.6) can be written as 
(5.8) So(A, 1) =E 2 DI Oy| pal «© Fpggal [(241+ 1)? — 07] [ (22+ 2%—-+ 1)? — 07] * 
J [(2ty oe 420 +2)? 07]. 
f (5.8) is substituted in (2.28), the determining equation is obtained: 


— o?+ > M* Gy) On| 06.9 lees 441 [(2%,+ ial = Fa We [(2%,+ 2Ho + T) = vl 


5.9 x 
( ) (*) : [ (2 Ar | 2%; 4 a2 il 
or in developed form 

2 
Py espera, 90) lz| 
To =, Ms Oar (2%-+1)2?—o? af 
(5.10) des 43 ; wel pal Obes — %2| / ‘3 (A4) . 
ae ree 
Ser [(2%, + a7] [ r( (2%,-+7)?—o?] 
Hy + Hy 
I'rom this equation one derives 
A o2 23 O\, | 91 ..,| oh, — xy) 4 
5.41) t=4+ Ie ED ate Il bel “Pe L(A 
( ) 5 402 2—o? | 326 Fis ala Hy Yo (24 AL Gy (x,+o) ae 
Ky Hy 
as the characteristic exponents exact mod (A*). 
If t is equated to 0,1, 2,... in the determining equation, one obtains the 


relation between o and A for which the resulting HILL’s equation has periodic 
solutions. 
Equation (2.26) for ay) yields, as in the case of MATHIEU’s equation, 


(5.12) Go \0 >t. Ot). 


With the two values of t from (5.9) this represents two linearly independent 
vectors except when t=0. To complete the solution (2.29), we use (5.5) again. 
Then (2.37) becomes 


(5.13) RA, 752) = QM Oy, Ope) «++ Seq | F(%) F(x +4) «.. Pte ie Heal teeny Fey ee 
(#) 


where x, x, ..., %;41 (7 =0, 4, ...) are positive or negative integers and the sum 
is extended over all systems (x, %,..., #;14) Subject to 
(S044)5 mecpepah eee 05, eee ieee sb ee Os Ae ie 


Using (4.19), we may write (5.13) as 
DAU, RS cers » A Oh Opaline= Sgn [(2% + 2%, + t)? —o?]4x 
(5.15) X [(2% + 2 + 2%. + 1)? — o?]4 
. ((2% + 2% ees 2%, + t)2 = 07 | 1 F(x) ett 


If (5.12) and (5.15) are substituted in (2. 29) and (4.21) is used, the following 
expansion of the solution to HILt’s system (5.2) is obtained: 


W,(2 jog ZEA eis 
(ey) (eure Siege Fyn) Bal al Sagal X 


\ Ws (2) —% () \— (2%-+4 +0)7 


X [(2% + 2% bit)? o8] 1, (Dron, ee. 226; 0)? otjiaie, 


(5.16) 
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Forming (w,-+-W»)/20 yields the following solution to Hi1tr’s equation (5.1): 


u(z)=e'* Mt ae Ax M Oe Pegler poy 4 [(2%-+1)?—o?] 4 [(2x+2x,+1)?—0?]4 


(5.17) : 
eee my. 2h, meio 


Using the two exponents t, — 7, one obtains from this the two solutions (denoted 
again as cu(z, A), su(z, A) 


CU 2,4) == COS aa? 4A DB op) Opes) «++ Ogg] (2% + 7)? — 7] 2 X 
(*) 


(5.18) Suz, Ay — sin t2 
x [(2%+2%,-+1)?—07] 4... [(2%+2%,+---+2x;+1)?—0? | 


= 08(2%- T) 2 


sin(2x+1)z 
There is a result which generalizes formula (4.28) to the solution (5.18) of 
HILv’s equation. Put 


(5.19) cu(z,A) =costz +4 >) c,cos(2% + 1) z. 
es 


Then by (5.14) and (5.19) 
3 DI Op, (Cy ae C_y) = > Ma Oo}! lx, | +++ Dlr 4 1 [(2% =F t)* na o? | s 


(5.20) “rE (%) 
.»» [(2% 4+ 224+ -+--+ 2%, + 1)? — o?]44 

with the systems (x, x, .--, %;,) subject to condition (5.14). Comparing (5.20) 
with (5.9) we find 

1 Co Oe 
(5.24) z Di P(r +4) = a 

kE1 

miichereuuces tO (4.20) ito; — 1, o5=—Gg=~ = ==0: 


6. Convergence of Expansions 


Theorem 6.1. Assume the matrix B(z) of the system 


(6.1) __K — =Aw-+A/B(z)w 
has the expansion 
(6.2) Biz) = By +) By eo 
ke 
absolutely convergent uniformly in 8. Let 61,09, .--, Gy, (m<n) be the values for 


which the matrix oK —A 1s singular and assume 


(6.3) inf |N-P+7-—g,|>0. 
1 m 


Then the series of Theorem 2.1, 


(6.4) RA, T; 2) rT. 2 Ne (Oe BS ate Cre Be eee. 
(7,8) 
(6.5) Si(Apos 2) a Dome Be On) Pacha Naas 2 


converge absolutely uniformly in % for sufficiently small |A|. 
20* 
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Proof. Hypothesis (6.3) implies that each of the matrices 


(6.6) C,=(N) P +7) K—A, rEM 
is non-singular and moreover that 
(6.7) y =sup||C;*|| <0. 

rEn 
Put 
(6.8) b(2) = TB — 2 Bill lene 

kel k=l 
and let 
(6.9) |o(2)| = bla) S Bet FCB. 
The series 
(6.10) @ (2) =a Dy" \2"| 8" (Zz) 
n=0 

converges absolutely uniformly in ¥ if 
(6.11) [Al <(By)>, 
and w(z) satisfies the equation 
(6.12) @=at+ Aly d(z)o. 


Let the expansion (6.8) be substituted in (6.10); let the multiplications of 
exponentials be performed according to 


(6.13) ee a) il ee Saeed ee eh N,+N=N,. 
Let the terms with | e’*’?*|, 7 fixed €9, be combined to form the term 4,| e'%”"P*| 


and the remaining terms with |e'’?*|, y€%, be combined to form one con- 
stant term a). The resulting series 


(6.14) @ (z) elk oa 
rgEM 


is, like the series (6.10), absolutely convergent uniformly in %. 

The coefficients a, can also be obtained in a recursive way by substituting 
(6.14) and (6.8) in (6.12), rearranging as above and comparing coefficients. One 
obtains 


(6.15) 4, =|Aly 2 I Ball ae, Noy + Ny =N,, PEs 


where the sum is extended over all the k, k’ for which Ny) + N,=N,. Comparing 
(6.15) with equation (2.9) according to which the coefficients a, in the expansion 
(2.8) of the solution w(z) are determined, one concludes that 


(6.16) \|a, || 4, rer 

provided |! a || <a). Thus, the series in 

(6.17) w(z) = (40 vas ay gite 
rgm 


which is identical with that of (6.4), is dominated by the series in (6.14), hence 
is absolutely convergent uniformly in 8. 
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From (6.12) it follows that if the product D(z) @(z) is expanded with the use 
of (6.8) and (6.14) in the way described above, the resulting series is absolutely 
convergent uniformly in 8. This series dominates that of (6.5). Thus, the asser- 
tions of the theorem are proved. The only restriction on A is the inequality (6.11). 

In the case of periodic B(z), where the 0, (k =1, 2, ...) are integral multiples 
of some number w, condition (6.3) will be satisfied for all t sufficiently close to 
one of the o,,...,6,, provided 


(6.18) O,— 0,4 %0, TD mek On ee, Amana es Op ae 


This is the minimum condition, used also in the preceding sections, to make the 
C, (ry €M) nonsingular. Condition (6.3) is satisfied also for all 7 sufficiently close 
to one of the o,, ...,0,, if the o, (k =1, 2, ...) are real, but none of the 0, ..., 6,» 
is real. It is not satisfied if the o, are ea but not integral multiples of some w 
(case of almost periodic nonperiodic B (z )) and at least one of the o,,..., 0, is 
real. Whether convergence of the expansions (6.4), (6.5) can also be proved under 
these conditions is doubtful. 


7. Characteristic Exponents 


Whether there are m or fewer linearly independent solutions of (2.1) which 
have the form (2.8) depends on how many solutions the determining equation 
F(A, t)==0 has. For A=0 the roots of this equation are .¢,,09,....;O, (¢w=n) 
which are the zeros of det(aK —A), and for sufficiently small A each of the 
zeros of F(A, t) will lie near one of the o,. Let U, be a neighborhood of oa, 
(uw=1,2,...,m), U the union of the U,, chosen so that the condition 


(7.1) epee en 
vEN 
implies 
(72) C= nf |N,.P+7—a,|>0. 
ses 1 
rEN 


We look for zeros of F(A, t) in U, so that the condition (6.3) of Theorem 6.1 can 
be replaced by condition (7.1). 

Theorem 7.1. Assume the matrix tKk —A 1s singular for the m distinct num- 
bers T=01,0p,---,6,, and U is an open set of the t-plane containing the points 
G15 Oz) 2. », O),. Assume the series 
(7.3) Seca AC nD mn Ba 

(7, 8) 


[the sum is extended over the systems (11, ..-, 7, So» +++, 84) restricted by (2.20), 
(2.24)] is absolutely convergent for some A+-0. Then the determining equation 


(7.4) F(A, t) = det (1K — A —AB,— 2?S,(4,2)) =0 


has, for all sufficiently small A, m roots t,, Tz, +++) Tm, 1% U and there is a null 
vector ay, of the matrix (t,K —A —ABy—#So(A, t,)) (uw =1,2, -.., 1m) such that 


Agi» 4o2>+++» 4m ae linearly independent. 
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Proof. The function /(0, t) has the zeros tT =0j, 02, ..-, Gp- If these are of 
order 1 then 


é $ 
(75) 2F(0, theo, = detleK = Aly, 0, 9p SH} 


and, by the implicit function theorem, the equation F(A, t) =0 has, for sufficiently 
small 4, a root t=t,, in U which is an analytic function of A and is of the form 


(7.6) T, =o, (mod A), = 1, ee 

Since the system 

C7) (cK—A —ABy—# Sold, t)) ay = 0 

has for 4 =O the solutions t =o, 05, ..., 6,, With m linearly independent vectors 
Gg = Cty Cay nem tt Das, fOr Asurucicntny, small, 7 = T,, 2 Solution dy)= 4p , which 


is an analytic function of A and is of the form 


(7.8) ay, =e, (mod 4), eE=Ne2nek we 
Thus the Gramian determinant of a,, ..., 4, iS an analytic function of A and 
+0 for A=0. It follows that the vectors ao,,..., 4,, are linearly independent 
for sufficiently small A. 

Assume next thatthe zero o, (#=1, 2,....,m) of F(O,t) is of order g,, 


1<q9,<n. By a well known extension of the implicit function theorem there 
is, for sufficiently small 4, a zero t, of F(A, t) in U which is an analytic function 
of A” and is of the form 


(7.9) T, = 6, (mod AM%) , Ue ANS? a aie 

System (7.7) has, for 4 sufficiently small, t=1,, a solution ay=a, (u=1, 
2,...,m) which is an analytic function of 4!“ and is of the form 

(7.10) Aon =e, (mod AMZ) , [bMS ors Oe, 

where @,,é@),...,é,, are constant linearly independent vectors. It follows as before 
that the vectors a ,, a2, ..., 4», are linearly independent for A sufficiently small. 


The following observation for the common case where By=0 is worth noting. 
In this case both equation (7.4) for t and equation (7.7) for a) contain no terms 
which are of power one in A. Therefore, (7.9) and (7.10) become 


(7.44) %,==o,, (mod) 1744) sag e,, (mod 7/1) , WO, 
The linear independence of the vectors a1, Qo, -.., 4m implies the linear 
independence of the solutions w,(z), w2(z), ..., w,,(z) of system (2.1) which are 


obtained by putting t=1, and ay=ay, (u=1, 2,...,m) in formula (2.29). 
Combining results of Theorems 2.1, 6.1, and 7.4, we obtain the following: 


Theorem 7.2. Assume the matrix B(z) of the system 


1 7 dw ; 
(7:42) ; K Gp Aw +A B(2) w 


has the expansion B(z) = By +)  B,e'%’?* absolutely convergent uniformly in 8. 
Assume these 


(7.1) inf |N,-P+o,—«,|>0, 


y reey 
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where 01,02, ..-, Om (m<n) are the values of o for which the matrix oK —A is 
singular. Then, for A sufficiently small, there exist m linearly independent solutions 
of (7.12) with the expansions 
(2) (40. +>) a, iN Ps) ad ae HN aD teins HE 
. rEM 
absolutely convergent unifornily in 8. The characteristic exponents t, and the vectors 
ay, are solutions of the system 


(7.7) (cK — A —ABy— 2S, (A, 2)) a =0, 
and the vectors a,,, are given by 
(7.13) Ay, =AA,(A, Ty) Bou, rEM. 


Explicit expressions for A,(A, t) and Sy(A, t) are given in (2.22) and (2.25). The 


exponents t, and the coefficients ay, 4,, ave analytic functions of AM%u, and 


T, =o, (mod Aven) , 


(7.14) Ay, = &, (mod AV%) , 


Mu Ht 


tin == © (mod A!™l/4x) , AA == AMO) ey mus) Oe TUE 


where (o,.K —A)e,=0 and q,, 1s the order of the characteristic value o,. If By=0 
then 1/q,, 1s replaced by 2/q,, 1m (7.14). 

To illustrate the occurrence of characteristic values o, of order >1 consider 
the equation of HILL’s type 
au 
az 


(7.45) +(09 +29, cos 2n2)u =0. 


n=a1 


The developments of Section 5 do not apply to this equation since o =0 in this 
case. In place of substitution (4.2) we put 


(7.16) W,=U, Wye=—1H’ 

and obtain the system 

eh) 1 i= Oy, —iwh=A(o,+2Yo,cos2n2) 1, 
n=1 


which is of form (7.12). The only characteristic value is o,=1 and its order is 
g,=2. Condition (7.1) is satisfied. Proceeding as in Section 5, one finds the 
determining equation 
t=hAoyt ~ Ae Freg| Ocal +» Flog al (2% + Tb) (2%q  2%9 Aes 

ged (eg erieee ae 2 te) Pa 
with the x; as restricted by (5.7). Thus, the characteristic exponents are analytic 
functions of A? if j= 0 and of 4 if oy=0, with the expansions 


(7.18) 


Re A 2 5 
(7.19) T= tOR/? (1 + Gan >, | (mod /?), oy) 0, 
#21 
— 1 2 : A 1 Oy, 5 Fhe, Ore, Olsey— tee 
r=ti(tyaiiii4t(iys ae) 
2 % BENZ x ; KY HE 
(7.20) “21 | /21, [eco] 21 


Hy Hy 


(mod A); 0, = 0. 
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With these values of t system (7.7) has the solution 
(7.21) “ga ihs ee 


The solution of system (7.17) is readily completed. Its first component is the 
solution of equation (7.15) 


u“(z) = veld +A DS MGiaiiaes Cine (2% + 1)? (2% + 2%, + 1)? 
(7.22) (x) | 
ss (AE ae a oa zy atin 


with the x, x; subject to condition (5.7). Using the two exponents t, —t this 
yields the two solutions 


CU(24/) COS: 72 


SA Dy A Open AO aap Ier Ek aa pe ee 
S _» COS(2% + T)z 
sin (2% + T) 2. 


(7.23) su (2x Ayr=-sitea 2 
vin (2s 1 te D3 St) 


Finally we observe that if the set 9¢ is vacuous, that is, no linear combination 
of the 0,, @2,... with non-negative integral coefficients is zero except the trivial 
one, then the determining equation is very simple. In this case equation (2.24) 
has no solution, and therefore S)(A, t) =0. The determining equation is 


(7.24) det (¢K =A —AB,) =0; 


and if B)=0, then the characteristic exponents t,, coincide with the characteristic 
values o, for all A. The expression for R(A, t; z) [see (2.26)] is also simplified 


in this case: there is no restriction at all on the non-negative integers s,, Sy, ..., S; 
in the sum 
(7225) Rina) ae eee ee 


(7,5) 


and the 7; are determined from the s; according to 


(7.26) Ni Map Ne cy reas t=41,2,...,f. 


8. Stability and Boundedness 


In this section the independent variable z is assumed to be real and is denoted 
by x. It ranges either over the whole real axis denoted by R or over the positive 
real axis R,. It is assumed that 


(8.1) B(x) = By+ >) Bye %* = By + > B,etNw Px 
k=1 k=1 

is absolutely convergent uniformly on R,, hence 

(8.2) Im 0, => 0, Reeels, 

We discuss the stability of solutions of the system 


(8.3) a e =Aw+AB(x)o 
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which are of the form 
(8.4) Ww (x) = (2, =e a a, NP) ete 

rE 
where it is assumed that the series expansion converges absolutely on R,. We 
also assume that the expansion (2.25) for S)(A, t) converges absolutely so that + 
can be determined from equation (2.28). As proved in Theorem 6.1, these hypo- 


theses are satisfied for sufficiently small A in the important special case of periodic 
B(x), provided 


(8.5) CO, — 0, >= 4W TRE eas) Wy ee ear Ae ta OPO 
where 0;,...,6,, (mn) are the distinct characteristic values which make the 


matrix ok —A singular, and 2z/w is the period of B(x). In this case there are 
at least m linearly independent solutions of (8.3) which are of the form (8.4). 


For linear systems, stability and boundedness are equivalent. Clearly solution 
(8.4) is bounded for sufficiently small 4 on St, if and only if 


(8.6) Ini 0: 
The characteristic exponent 7 satisfies the determining equation 
P(A, t) = det (« K — A. — ABy— 27 S5(A, 7)) =0, 


where S,(A, t) is given in (2.25). In this expression only those coefficients B,, 
of (8.1) are used for which condition (2.24) is satisfied. Thus, B, does not affect 
the determining equation if Im @,>0. Therefore, in deciding whether solution 
(8.4) is bounded or not on #,, the terms B, exp(io,x) of (8.1) with Im o,>0 
may be ignored. For this reason we assume from now on that the o, (k =1) are real. 

For sufficiently small A the characteristic exponents Tt lie in the neighborhood 
of the characteristic values o which make the matrix of —A singular. In fact, 
by (7.9), if a, is a characteristic value of order g,, then there is a characteristic 
exponent t,, which is an analytic function of 4"/% and for which t,=<, (mod 4%) 
holds. It follows that the boundedness condition (8.6) can be satisfied for suf- 
ficiently small 4 only if 


(8.7) Rin Gpes Oe 


we 
If Imo,>0, then condition (8.6) is satisfied; hence solution (8.4) is bounded 
on %, for A sufficiently small. In the following only the more difficult case of 
veal o,, is discussed. We prove that the characteristic exponent 7 is real, hence 
soution (8.4) is bounded on , in several special cases. 


Theorem 8.1. Jf K,A,B,; (7=0,1,2,...) are real matrices, then, for all 
sufficiently small real i, there are at least as many linearly independent solutions 
of system (8.3) bounded on NR as there are real characteristic values o of order 1 which 
make ok — A singular. 


Proof. Let o,,05,...,0, be distinct real characteristic values of order 1. 
Then there is a unique root t=t, of the determining equation F(A, t) =0 for 
which 


(8.8) lim a 
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The determining equation is easily seen to have only real coefficients. Therefore, 
T,, is also a root, and by (8.8) lim t,=0,. Thus t,=T,, T, is real (4 21 2c 
It will be noticed that B(x) = By+ > B,e'%* need not be real in Theorem 8.1. 
If B(x) is real, then the condition that the B; (7 =0, 4, 2, ...) be real is equi- 
valent to B(x) being even, B(x) =B(— *). 
In the following theorem the matrix B(x) is supposed to be composed of 
blocks which are alternately real and purely imaginary. More precisely, assume 


(8.9) B; ee (Bayete paleheys Mm? i = 0, A; 2, ae 


juny is ann, Xn, matrix (m,21; m+ M21 +n,,=%) which is real 
for «-+y even and purely imaginary for “ +¥ odd. 


where each 6, 


Theorem 8.2. If K, A are diagonal real matrices and B(x) 1s as in (8.4), 
(8.9) with 


(8.10) B 


then there are at least as many linearly independent solutions of system (8.3) bounded 
on St as there are characteristic values o of order 1 which make oK —A singular. 


ip» veal for w+y even, B,,, tmaginary for u+yv odd, 


Proof. For convenience let us say that an” x m matrix H is “‘essentially real”’ 
if H =(H,,)u.=1,2,...,.m» Where each H,, is an m, Xm, matrix (m,21; m+M%+-:: 
+n,,=%) which is real for ~+y even and purely imaginary for w~+~y» odd. 
It is readily verified that if D is a real diagonal matrix and H, H’ are essentially 
real matrices then D+H, DH, HD, H +H’, HH’ are essentially real matrices. 
Since by hypothesis the matrices B; are essentially real, one concludes by applying 
the preceding remark several times in succession that the determining matrix 
tk —A—AB,—A?S,(A, t) is essentially real for A, t real. The determinant of 
an essentially real matrix H is real. To see this let the number 7 be factored 
from each of the rows numbered 2,+1,...,%+%; m+m”.+73+1,..., 
My +Ny+N3+MN4;... and also from each of the corresponding columns. The 
remaining elements of the determinant are all real, hence det H is real. It follows 
that f(A, t) is real for real 4 and Tt. 


The remainder of the proof is the same as that of Theorem 8.1. 


Remark. In the important case of real B(x) condition (8.10) is equivalent 
to B(x) being “essentially even’; this is to mean B(x) == Balx) es pk ae 
where each B,,,(x) is an , x, matrix whose elements are even functions of x 
for «+» even and odd functions of « for «+» odd. 

In the following theorem some of the B ; are real and some are purely imaginary. 
We rewrite B(x) as 
(8.41) B(x) = By + D Bei’ 47> Bl eieb’s 


kE1 kE1 


’ 


where By, BY), BY), ol), o®) (k=41) are real. 


Theorem 8.3. If K, A, By, BY, B®) (k=1) are real matrices and if > ny oD + 
Ding of” =0 (each n\), n®)=0 or 1) implies Yn” =0 (mod 2), then there are at 
least as many linearly independent solutions of system (8.3) bounded on KR as there 
are veal characteristic values o of order 1 which make oK —A singular. 
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Proof. The only possible non-real coefficients in the determining equation 
F(A, t) =0 arise from S (A, t). The typical term in this expression [see (2.25) ] 
is for the present case 


(8.12) Geb Okay Bide Cle) BOE 
where each ¢; (7 =0,1,..., f) is either 1 or 2 and / is the number of 2’s among 
the ¢;. By (2.24) 
(to) (Ga) i oie Ge= 
(8.13) Qs, ae Qs, a a 0; Sw 


and, by the hypothesis of the theorem, this implies 4 =0 (mod 2). Therefore 
(8.10) is a real matrix for real T. 


The remainder of the proof is the same as that for Theorem 8.1. 


Remark. In the important case of real periodic B(x) let 


(8.14) B(x) = By+ >) BY coskwx + >) B® sinka x 
kE1 k=l 

with w>0, By, BY, B® (k=1) real. Then if 

(8.15) BY) = 0, BY) =0, (Mak i Oxo de 


the hypothesis of Theorem 8.3 concerning B(x) is satisfied. 


Theorem 8.4. Jf K, A, B(x) are (complex) Hermitian matrices then, for all 
sufficiently small real i, there are at least as many linearly independent solutions 
of system (8.3) bounded on N as there are real characteristic values o of order 1 
which make oK —A singular. 


Proof. In (8.1) we may assume that 0, @,, 02, ... are distinct numbers. Then 
since B(x) is Hermitian [that is, B(x) is equal to its conjugate transpose B*(x) |, 
we conclude that B,=B¢ and that if B,exp(io,x) is a term of (8.1) so is 
By exp(—io,x). We set 


PE Or. 01; 0a. — O20 +)» 0, > 0, 
(8.16) Or = Ng): JE (Qs = Ney: 12, k = 4h 
ING Peet IN all, rE. 
Then 
(8.17) Bota 1d, eee ee eee 


If t is a characteristic exponent, then for real A 


(8.18) det (cK—A—A By—2? Sy (A, 2)) =O=det (7 K*— A*—A BE —7? SF (A, 7)). 


(8.19) Sade) Te ath Siete ain? 
where the So, ..., 5;,%,-.-, 7 are restricted by conditions (2.20) and (2.24). By 
(2.17) and (8.16) since K = K*, A =A* 


(8.20) C*=(N,-P+7)K—A, rN. 
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Using (8.17) and (8.20), one may write for (8.19) 
(8.21) Sti t =D Aa Be Cae By CxB ss 


(7, 8) 
Now by (2.20) and (8.16) 
(Noy + Mga) + ++: + Noy) P, 


— (Nea bo + Ney) P, 


8.22 
(8.22) Ze (RiBy yet aN a) igs 

me N «2 Pe ea ee Oy tee ef ade 
This shows that af the systeni<sp, sj, -+-ss77-!15 - >> 77 satisfies conditions (2.20), 
(2.24) so does the system sj, Sa) (23)'S05 73%) Me and! vice-versa) Thas the 
sum in (8.21) is the same as in (2.25), and we have 
(8.23) So (A,2)c lA a 


It now follows from (8.18) since K*=K, A*=A, B*=B, SoA, 2) Sela) 
that 7 is a characteristic exponent if T is. 

The remainder of the proof is the same as that of Theorem 8.1. 

If Theorems 7.2, 8.2, 8.3, 8.4 are applied to the system considered in Section 3 


a2 “ 
He + pt, +2A > dyy(x)%=0, fe =4,2,....m, 


ad x 


(8.24) 


y=1 
where the real b,,(x) have the absolutely convergent Fourier expansions 


(8.25) buv(*) = Ouy.0-+ >, 00) ,cosk@ x >) OM) sin kava; 
kE1 kD1 


the following theorem is obtained. 


Theorem 8.5. Given o,>0, Cyto, xO, 206,--xW (U,y=1,..., mM, MEY; 
~=0, +1, +2,...), then for sufficiently small real i every solution of system (8.24) 
1s bounded tf one of the following conditions 1s satisfied: 

(a) The matrix (On neers ae m US Symmetric. 

(b) The matrix (b,4(%))u»=1,...,m 2S essentially even. 

(c) Oy 2}—1 =O Oy. 2¢= Os oak 

The sufficiency of condition (a) and of the special case of condition (b) where 
essentially even is replaced by even was proved first by CESARI ([1]; see also 
[7], [8]). Condition (b) is equivalent to one given by HALE [9] (it can be derived 
from HALE’s by a convenient change in the numeration of the unknowns and of 
the equations). If x is replaced by x—a/2@, condition (c) is reduced to con- 
dition (b). 

The following theorem is of a different nature. It deals with boundedness 
on ht, rather than on jt. For this reason the conditions involved are inequalities 
rather than equalities and symmetries. These inequalities may be satisfied for 
all sufficiently small positive (or negative) 4, but not necessarily for all suf- 


ficiently small real 2. If the characteristic exponent near the real characteristic 
value o is 


(8.26) CH= OLD 


TZ 
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and if 
(8.27) Inyo, lms === 1m 0, f= 0; Im 6,>0, 


then Im t=0 for all sufficiently small positive 4° and the corresponding solution 
w(t) is bounded on }t,. Hf, on the other hand, Im 6,<0, then, for all sufficiently 
small positive 4°, w(t) is not bounded on %,. 

In the following K and 4 are real matrices while the B, (j =0, 1, 2,...) are 
arbitrary complex matrices. The characteristic value o which makes oK —A 
singular is assumed to be real and of order 1. Let e be a real null vector of 
ok —A and e* a real null vector of oK*— A*, that is 


(8.28) (ck —A)e=0, (CK 3A \icr == 0: 
We assume 
(8.29) (Ke e*) == 0; 


where the bracket symbol denotes the scalar product. If K and A are symmetric 
matrices and A is moreover positive-definite, as they are for the systems con- 
sidered in Section 3, then e* may be taken as e and (8.29) is satisfied since 
(ice we)--— Oviimplies (4-¢-¢) —O; hence e —0. 

Let all the terms in S)(A, t) [see (2.25)] containing the same power of / be 
combined and put 
(8.30) ABy HA So(A, 4 ea MDa t) +1D;(z)), 

Yr 

where D,(t), D/(t) (y=1) are real matrices. These can be calculated from the 
matrices A, K, By, B,,..., and they are polynomials in the matrices By, B,,... 
if only a finite number of these are +0. 


The characteristic exponent t is a root of the determining equation 


(8.31) F(A, t) = det {7K — A — AE. (t) +4D,(x))} = 

Let the root t=t, which lies near o for small A be expanded in powers of A: 
(8.32) T= o + DA (5, +45), 

where 6,, 6, (y=1) are real es RT now assume that 

(8.33) Dye Do D2 5 (t) 0 

for all veal + near o and some s=1. Then it follows readily that 

(8.34) eee il 0, 

Also if we expand 2By)+A?So(A, t) in powers of A, 

(8.35) A Boy + 2? So (A, To) epee ,ttE,) 


where E,, E; (r=1) are real constant matrices, it follows that 


(8.36) ee ge ep) a =1).(0). 


Thus, EZ) can be calculated without the knowledge of the root 7,. 
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The matrix in (8.31) has for 7 =1, a nullvector /,, which may be expanded as 
(8.37) fo =Co t+ DM (h, +ih,), 
Peak 


where h,, hi, (y=1) are real vectors. It follows again from (8.33) that 


(8.38) hice hpertnicigile jee, 


Expressing the fact that /, is a null vector and using (8.34), (8.36), (8.38) one has 


(+ 044, 4146) K —A— SNE, —i4E,)- 


(8.39) i. ie . 
(6, Ei path) =O) med 
id ‘ 
Equating the imaginary part of the coefficient of A° in (8.39) to 0, one finds 
(8.40) (0.K —E)e,+(oK — A)h, =0. 


If the scalar product of this vector by eZ is formed, one obtains 


1 _ (Egeg,e*) _ (Ds (2) eo, 26) 
a) @, (Keg, et) ee eter 
Thus by (8.34) 
(8.42) Im 7, = (Ps) 0x08) 9s" mod AS+4). 


(Ke, , es) 


which can be calculated without knowledge of the exponent t,. If (8.42) is 
positive, then the solution w=w,(x) of (8.3) is bounded on %, for sufficiently 
small A°, and it is not bounded on Nt, if (8.42) is negative. Thus we have proved 


Theorem 8.6. Suppose K, A are real matrices, the characteristic value o of 
order 11s veal and o K —A,o K*— A* have the real null vectors e, e* with (Ke, e*)=- 0. 
Let the matrix whose determinant is the determining equation of system (8.3) be 
expanded as im (8.31) and suppose that 


Di (x) = Di(t) =--- = Di a(x) =0 


for all veal t near o. Then for all sufficiently small A with 28>0 the solution 
w=w,(t) of (8.3) belonging to the characteristic exponent t, is bounded on MN, if 


Di(0) eg, e#) 

8.4 ( s o» &o 

see (re 

is positive, and unbounded on Ji, if this quotient is negative. 


To illustrate this criterion consider the system 


du, 2 : 
| 
fo 0] % = 2A(sin@ X) uy, 


2 
aus 
dx? 


(8.44) 


+ 03 U, = 2A(cosw x) uy, 
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where 0,>0, 0,>0, w>0, A real. By the substitution 


y= —O,% +i, Wy= — 0,4, —in,, 
(8.45) 


Da 7 
W3 = — OnU2 + 1s, Ws = — OnUy + 1s, 


system (8.44) becomes 


2 w, =o7u,+4 a (sin w %) (ws + w,), 
2 

= 1 ws = 01 Wy aa (sin w x) (ws + w,), 
(8.46) ; 


“2 ig = G5, 1 A i (cos w x) (w, + Ws), 


- Op == O50, A “a (cos w x) (w, + wy), 
1 


which is of form (8.3). Condition (7.1) becomes 


(8.47) 0, to,+xo, = se 1, ns 


and is assumed to be satisfied. The matrix B(x) in this case is B(x) =B,e'®?*+ 
B_e~*®* where 


ee ( Os — 1 (04/02) i ue ( Os t (0,/09) it 
(62/0,) E Oz, (6/0) E Oz 


where 0, is the 2 x2 zero matrix and E is defined in (4.5). Let the matrices C;* 
belonging to N,- P=w, —w be denoted as C+, Cz=' respectively. They are 
diagonal matrices whose elements are +(-+-w+1)o,—o} (i =1, 2) denoted as 


vy; ,¥; (j=1, 2,3, 4) respectively. Thé matrices D(z) of (8.30) are Dj(t) =0, 


; CORRE Bea ewe. Oy 
8.48 DA — Lh 
ae aM | Oz (—yt —ywe +yi +2) £ 


The characteristic values making ok —A singular are o,, —0o,, 063, —@_ (not 
necessarily distinct) and corresponding characteristic vectors are e; = (1, 0, 0, 0), 
é, =(0, 1, 0, 0); €5 =(0, 0, 1,0); €2 =(0, 0,0, 4)= The quotients (8.43) become 


(D3 (0) er, Cue ee (D3 (= 0;) eT, ei) ney eee fy Wiens 
(Key, ef) (A ez, er) [w? — (0, +6,)?] [w? — (a, — o4)?] 
8. 
Be ie a COME 6) ge) EA were ed os 
(Ex es, e3) (K ez, e2) [@7 =~ (0,77 0)" |" [(@ — (0, — 05)" )" 


We conclude that the two linearly independent solutions of system (8.46) cor- 
responding to ey, ey are bounded on Wt, for all sufficiently small / if 


(8.50) |o,— o,|<w<0,+ 0, 
and unbounded on %, if 


(8.54) wo<|o,—o,| or w>0,+ 65. 
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For the two linearly independent solutions of (8.46) corresponding ‘to €3 , €3 the 
situation is exactly reversed. Thus, for no w are all the solutions bounded on 
jt, no matter how small 4 is. The same conclusions hold for the solutions of 
system (8.44), which may be chosen as real. That system (8.44) has unbounded 
solutions was discovered by CEsARI [J] for particular values of o,,0, and by 
GamBILL [6] for the general case. It is to be noted that the transition from 
boundedness to unboundedness occurs at the highest “‘resonance frequencies” 
@ =0,-+0, and w =|0,—0,|. No comparable changes occur at the other resonance 
frequencies excluded by condition (8.47). It is also to be observed that it was 
not necessary to assume o,=- 02 since the linear independence of the characteristic 
vectors e},e],e}, 2 assures the linear independence of the corresponding 
solutions. 
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